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CRUISE NARRATIVE 


Highlights 


Area Designation: 


Expedition Designation: 


Chief Scientist: 


Ship: 


Ports of Call: 


Cruise Dates: 


Mission Summary 


Overview 


NAFO Regions: 5Y, 5Ze, 4X, 4W, 4Vs, AVn, 3Pn, 3Ps 
Extent: 41° 51'N - 47° 35'N; 056° 08'W - 066° 08'W 


EN2017606 or 32EV17606 (ISDM format) 


Dave Hebert 

Ocean Ecosystem Science Division 
Marine Ecosystem Section 
Department of Fisheries and Oceans 
Bedford Institute of Oceanography 
PO Box 1006 

Dartmouth, NS, Canada B2Y 4A2 
David.Hebert@dfo-mpo.gc.ca 


R/V Endeavor (call sign - WCE5063) 
Oceanographic research vessel out of the University of 
Rhode Island. 


Nov 24", 2017 — Depart BIO, Dartmouth, NS 
Nov 26", 2017 — Arrive BIO, Dartmouth, NS 
Nov 27", 2017 — Depart BIO, Dartmouth, NS 
Dec 5". 2017 — Arrive Sydney, NS 

Dec 5%, 2017 — Depart Sydney, NS 

Dec 9*. 2017 — Arrive BIO, Dartmouth, NS 
Dec 11%, 2017 — Depart BIO, Dartmouth, NS 
Dec 16", 2017 — Arrive BIO, Dartmouth, NS 


Leg 1: Nov 24^ — Dec 5" 
Leg 2: Dec 5% _ Dec 16^ 


The planned departure of the R/V Endeavor from BIO was planned to be at 1000 LT on 
November 24". An issue with steerage delayed departure until 1045 LT. The start of the 
recovery of the Nova Scotia Current Mooring (M1996) started at 1500 LT and the new 
mooring (M2024) was deployed at 16:30 LT. Then, we headed to HL 01 to start the 
Halifax Line throughout the night, completing HL 03.3 at 0630. The AMAR mooring 
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(M1949) in Emerald Basin was recovered on the 24" on the 25th. We headed back to 
BIO to wait out an approaching storm, docking at 1830LT. We departed BIO at 1830 
LT on the 26" and were running several hours late out to HL 04 due to sea state and 
weather. The weather was rough throughout the night so we waited until daylight before 
heading to HL 05. A release test was conducted just prior to the station occupation at 
HL 06. Stations occupations were then completed in order out to HL 06.3. At 0420 LT 
on the 28". the conditions became too poor to continue in the southeast heading so it was 
decided to change headings to begin mooring work in Dawson Canyon. The weather 
deteriorated further and a decision was made to hold position near HL 06.7 and heave-to. 
At O815LT, it decided to head to HL 06.7 and complete that station in addition to 
HL 07. Due to time constraints imposed by the planned port call in Sydney and 
impending weather in the area, the rest of the extended Halifax Line stations were 
dropped from the schedule. 


On November 29", a series of CTD casts were conducted in Dawson Canyon to await the 
deployment of the AMAR mooring at first light. The mooring deployment (M2027) was 
completed at 0930LT and we headed to Logan Canyon. The AMAR mooring at Logan 
Canyon (M2028) was deployed at 1600 LT and a CTD cast was conducted nearby before 
heading to the Gully. 


While occupying SG. 28 on November 30%, conditions deteriorated enough that we could 
not conduct a vertical net tow. After the CTD cast, we hoved-to until the weather 
improved. At the same location, a release test was conducted at 1630 LT. We deployed 
an AMAR mooring (M2026) at 1920 LT at the offshore Gully location and conducted a 
MCAL survey. The remainder of the Gully station occupations were followed by a 
recovery and deployment of an AMAR mooring (M1948 and M2025 consecutively) 
during the afternoon of December 1“. 


Following the Gully work, the Louisbourg Line was occupied. After LL 09, the AMAR 
mooring at Stone Fence (M1950) was recovered at 0930 LT on the 2". The Louisbourg 
Line was completed and the St. Anns Bank Line was started. After completion of the 
line, the AMAR at the end of the line was deployed (M2029) and another recovered 
(M1947) on the morning of December 4%. We headed to M1999, a St Anns Bank 
mooring that could not be recovered on an earlier mission. There was no communication 
with the release so it 1s likely that it was released during a previous mission in November. 
Later in December, the ADCP and SUB were found in Newfoundland. On December 5", 
we headed in to Sydney Harbour to disembark some of the science party (Jay Barthelotte, 
Adam Hartling, Jenni Tolman, lan Luddington and Jennine Winkel) and the chief 
engineer. After the change, the Cabot Strait Line was completed on December 6" before 
heading to STAB 06. At 1600 LT, winds increased and our heading limited our ship 
speed to 6 kts and later, to 4 kts. 


The line across the mouth of Laurentian Channel began in the afternoon of December pe 
The line was completed on the 8% and plans were made to head to BBL 01 to avoid an 
offshore storm. On December 8% at 1600 LT, the Captain decided the storm on 
Saturday/Sunday was too large to stay out given our close location to Halifax. We 
returned to BIO at 1400LT on the 9%. The ship departed BIO at 2130 LT on December 
10”. 
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The Browns Bank Line began at around noon on December 11%. The section across the 
Northeast Channel was occupied, but BBL 07 was dropped due to impending weather. It 
was decided to head to the western end of the Yarmouth Line. On the way to YL 10,a 
CTD cast was undertaken at PL 08 on December 13". The Yarmouth Line was started 
at 1430 LT on December 14%, At YL 06, communications to the CTD was lost as it 
started the upward portion of the cast. A decision was made to switch the CTD/net 
winches due to a shorted termination on Winch +1. Due to weather and timing, YL 03 
was the last station occupied before heading back to BIO on December 15", HL 02 was 
occupied before heading into Halifax. The RV Endeavor arrived at BIO at 1545 LT on 
December 16". 


Over the 23 day mission, the R/V Endeavor logged ~2861 nm and AZMP science staff 
conducted 175 operations at 87 stations (Figure 1). Table 1 breaks down the operations 
by sampling gear for each leg of the trip. The table also points to figures that display the 
deployment locations for each gear type. Each of these figures is accompanied by a table 
of coordinates detailing each deployment of that gear type. Table 2 contains the break 
down in time allocated to each gear type. 


*Note that approvals for work in the Gully and St. Anns Bank MPA are included in 
Appendix | of this report. 
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Figure 1. EN2017606 stations. Overlapping event labels may not visible. 
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Table 1. Station operation summary. 


Operation # of Operations Figure 
CTD 79 9 
Vertical Ring Net Tows 76 16 
ARGO Float Deployments 6 20 
Mooring Recoveries 5 21 
Mooring Deployments 6 2] 


Table 2. Operational time by gear type. 


Gear Operation Duration (hrs) 
CTD ~50 

Vertical Net Tows ~24 

Mooring Recoveries ~3 

Mooring Deployments ~2 

Argo Float Deployments ~] 


* Surface water parameters were recorded throughout the mission. Refer to the Underway Seawater 


system Section of this report for more information. 


Mission Participants 


A complete ship’s crew list for this mission can be found in Appe 


Table 3. EN2017606 Science Staff. 


Name Affiliation Duty Leg(s) Shift 

l Barthelotte, Jay DFO — OESD Mooring Ops l Day 
CTD 

2 Belzile, Mélany DFO — OESD Operator\Elog\Deck Both Day 
Ops 

3 Benjamin, Robert — DFO—PCSD Data Manager Both Day 

4 Caverhill, Carla DFO — OESD Lab Tech\Deck Ops Both Day 
CTD 

5 Cogswell, Andrew DFO — OESD Operator Elog Deck Both Night 
Ops 

6 Hartling, Adam DFO - OESD Mooring Ops l Day 

7 Hebert, Dave** DFO — OESD zd id Both Day 

8 Luddington, lan DAL — Erin Bertrand Lab Tech l Day 

9 MacIsaac, Kevin DFO — OESD Deck Ops\Biologist Both Night 

10 Perry, Timothy DFO — OESD Lab Tech Deck Ops Both Night 

11 Spry, Jeffrey DFO - OESD S E ina Both Day 

12 Tolman, Jenni DAL — Julie LaRoche Lab Tech l Day 

13 Winkel, Jeannine ECCC — CWS Piro ane Mammal I Day 


**Chief Scientist 
DFO: Department of Fisheries and Oceans Canada 


Observer 
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OESD: Maritimes - Ocean Ecosystem Science Division 

PCSD: Maritimes - Program Coordination and Support Division 

ECCC-CWS: Environment Canada and Climate Change - Canadian Wildlife Service 
DAL: Dalhousie University 


Objectives 


There were 15 defined objectives for EN2017606. Table 4 describes whether each of 
these objectives was met along with any relevant supporting commentary. 


Primary 


ll; 


Obtain observations of the hydrography and distribution of nutrients, 
phytoplankton and zooplankton at standard sampling stations along “core” 
Atlantic Zone Monitoring Program sections within the Maritimes Region 
ue Nir: ero M ada - http//www.bi0.gc.ca/science/monitoring- 


Additional 


2; 


Occupy stations in support of the extended Halifax Line (XHL) (HL 08 and 
greater) (Contact Dr. Igor Yashayaev) 

Carry out hydrographic, chemical and biological sampling at stations in the Gully 
in support of Gully MPA monitoring initiatives by Oceans and Coastal 
Management Division (Contact Dr. Dave Hebert - http:/Anter-w02.dfo- 
mpo.gc.ca/Maritimes/Oceans/OCM D/Gullv/Guiiv-MPA). 

Nutrients and hydrography across the Northeast Channel and Gulf of Maine as 
part of NERACOOS Cooperative Agreement, (Contact Dr. Dave Hebert - 
http//ww.neracoos.org/). 

Deploy 6 ARGO floats in support of the International Argo Float Program 
(Contact Dr. Blair Greenan - híiip//vwww.meds-sdmm.dío-mpo.ge.ca/1sdm- 
edsi/argo/index-eng. html). 

Collect underway and CTD water samples at specified locations and depths to 
fulfil the regional component of an Aquatic Climate Change Adaptation Services 
Program (ACCASP) initiative investigating the delineation of ocean acidification 
and calcium carbonate saturation Le ol n ooo zone — Dr. Kumiko 
Azetsu-Scott - htip://www.dfo-mpo.gc.ca/science/oceanography- 
oceanographie/accasp-psaccma/index- an 
Collect water samples for the Bertrand lab at Dalhousie University to evaluate 
whether and how organic and organometallic micronutrients influence primary 
productivity and phytoplankton community structure on the Scotian Shelf 
(Contact Erin Bertrand — Erin. Bertranc@dal.ca). 
Collect water samples from strategic locations and depths to support a microbial 
community analysis via DNA, RNA and flow cytometry, as well as the isolation 
of novel diazotrophs (Contact Dr. Julie Laroche - 
http //www.dal.ca/facuitv/science/biology/faculty-staff/our-facultv/1iulie- 
laroche.html) 

Bird and mammal observations as part of EC-CWS sea-bird observation program 
and in fulfilment of Gully MPA occupation requirements (Contact Carina 
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Gjerdrum - carina.gjerdrumíocanada.ca). 


10. Carry out hydrographic, chemical and biological sampling at stations in the St. 
Anns Bank MPA as a continued monitoring effort in support of Oceans and 

Coastal Management Division (Contact Dr. Dave Hebert - hitp//veww.dío- 

mpo.gc.ca/oceans/mpa-zpm/stanns-sainteanne-eng. htm] ). 

11. Attempt to recover a single mooring (M1999) deployed during the fall 2016 
AZMP shelf survey (HUD2016027). An unsuccessful attempt to communicate 
with the acoustic release was made prior to EN2017606 during Dr. Ed Horne's 
mission aboard the CCGS Perley (Contact Dr. Dave Hebert). 

12. Conduct hydrographic, chemical and biological sampling across the mouth of the 
Laurentian Channel (BP and BANQ stations). This transect has been proposed to 
enhance our understanding of hydrographic phenomenon in these areas in support 
of current modelling efforts (Contact Dr. Dave Brickman). 

13. Collect 200 um ring net zooplankton samples at 8 predefined stations across the 
Scotian Shelf to supplement the Canada C3 program sample collection (Contact 
Dr. Claudio Dibacco - bttps://canadac3 .ca/en/expedition/the-research/) 

14. Recover and deploy the Nova Scotia Current Mooring. This work, funded by 
AZMP, supports the operation of a mooring that continually monitors the Nova 
Scotia Current. These data are used to validate shelf circulation models. (Contact 
Dr. Dave Hebert). 

15. Recover 4 Autonomous Multichannel Acoustic Recorders (AMAR) from Emerald 
Basin, the Gully MPA, the Stone Fence Lophelia Conservation Area and the St. Anns 
Bank MPA. In addition, deploy a total of 5 AMAR moorings; 4 across the eastern 
Scotian Shelf break at Dawson Canyon, Logan Canyon, and the Gully MPA and 1 
deployed within the bounds of the St. Anns Bank MPA. (Contact Dr. Hilary 
Moors-Murphy - http://www.dfo-mpo.ge.ca/science/publications/article/2016/1 1- 

15-16-eng.html) 


Table 4. EN2017606 objectives status. 


Objective Status Comments 


2 Cancelled Due to early delays all XHL stations were dropped. 


Modified ` Dal could only participate for the 1” leg and thus was unable to make 
Complete collections from the western Scotian Shelf as originally planned. 


10 Complete 


Modified An additional station BP 00 was added to the NL shelf and all other stations 


12 Complete were successfully occupied. 
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SUMMARY OF ACTIVITIES 


CTD Summary 


Narrative 


As summarized in Table 1, there were a total of 79 CTD casts during the mission (Figure 
2 and Table 5). The configuration file used for the mission is provided in Appendix 3. 


At the beginning of the mission, DFO staff members were given a tutorial by the Ship's 
Tech on deploying and recovering the CTD off the starboard side of the vessel. 
Deployments and recoveries required 1 crane operator and 3 science staff. One science 
member was responsible for providing hand signals to the crane operator and controlling 
the swing of the CTD before it reached the rail. The other 2 science staff operated the tag 
lines on both deployment and recovery. On deployment, tag lines attached to the inboard 
rail of the ship, would loop the line around the vertical post of the CTD frame and then 
back to cleat mounted on the deck. On recovery, tag line operators used a long pole to 
secure a clip to metal extensions radiating from the CTD frame. Once clipped to the 
frame, they would put the free end of the line around a cleat and pull the line taught, 
which prevented the CTD from swinging as it was guided over the rail and into position. 
Once in position, the CTD was secured to the deck with ratchet straps and eyes screwed 
into the decks bolt pattern. 


The ship was able to deploy and recover gear in winds and waves comparable to the 
CCGS Hudson. Nonetheless, because of the low freeboard and dynamics of the ship, 
science staff were regularly exposed to wash on the deck both during recovery and 
deployment and also during water collection. 


Water sampling went smoothly but the ship was often required to hold position or steam 
slowly between stations. This impacted our program efficiency compared to our typical 
platform. The sampling area around the rosette was somewhat cramped and those 
sampling the starboard side of the CTD were often exposed to the wind and waves and 
were precariously close to a very low rail. Initially, water sampling was tricky, because 
sample bottle racks were brought out one at a time because of the risk of them being 
washed away or broken deck wash. Half way through the mission, a shelf was created on 
deck to accommodate racks while samples were being collected. The distance to the lab 
was minimal and made sampling more efficient. The proximity of lab space to CTD 
controls also made it simple for staff to stay in communication throughout the mission. 
Well positioned cameras around the decks of the ship allowed staff to gauge the state of 
operations. 


The CTD performed very well. Only 1 CTD cast (event 165 at YL 06) was aborted 
when the deck unit through an error. After some quick diagnostic work by the Endeavor 
Technician, it was determined that the best course of action was to switch the other EM 
cable and move the net to the cable which required re-termination. This meant a delay of 


10 


000010 


just less than 2 hours between the recovery of the aborted CTD and its redeployment. For 
the remainder of the mission the CTD and net were deployed from these new positions. 


Overall, science staff were pleased with the experience, competence and helpfulness of 
the Ship's Tech, Crane Operators, Engineers and Bridge Staff that made these CTD 
operations a success. 


Preliminary section plots of temperature (°C), salinity (p.s.u.) and sigma-t (kg/m?) can be 
viewed in Appendix 4. 


TOW SS W ge 


58"W SE W 


48°N 


44°N 


ENGOE CTE 
ENDGOO Mission Track 


42°N 


coo Protected Areas (includes Proposed) 


Figure 2. Locations for CTD casts during EN2017606. Each cast is labelled with the 
consecutive mission event. 
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Oxygen 


The oxygen data collected by the CTD sensors and Winkler titration method will be used 
to create new calibration coefficients before the final run of the CTD processing. It will 
be necessary to extract these corrected oxygen values when they are produced so they can 
be accurately reflected in our data archives. 


The adjusted Soc values are calculated by a 2 step process. First, a “threshold field” 1s 
produced that subtracts the mean difference between the sensor and the average Winkler 
value for all samples, from the individual sample difference between the sensor and 
Winkler: 


(SBE O2 - Winkler O2) - mean(SBE O2 — Winkler O2) 
The next step calculates a new slope term by using the following equation: 
NewSoc = mean(previousSoc*(| Winkler O2]/[SBE 02])) 


Before the Soc can be calculated however, comparisons between the primary (71230, 
calibrated August 2, 2017) and secondary (#0345, calibrated August 2, 2017) sensors 
were completed to remove outliers (Figure 3). The 1.5 * inter quartile range (IQR) was 
used to determine “outlier” data that could bias the results. The first 15 events (444601- 
444660) showed an average sensor difference greater than the rest of the mission (Figure 
3A). The difference during the mission also seemed to be changing slightly. There were 
some other minor removals, but another sequence of bad data was noticed during 
CSL 04 (Figure 3B). For oxygen sensors to be this close throughout the mission was 
actually quite good despite the removal of these outliers before proceeding to the next 
step. 


Comparisons were also made between Winkler replicates (Figure 4). There were a total 
of 7 Winkler replicates removed from further Soc analysis (events 43, 81, 95, 122, 124, 
163 and 174 which correspond to sample ID numbers 444795, 444999, 445062, 445241, 
445243, 445492, and 445547). The average difference between the Winkler replicates 
before outlier removal was 0.002 ml/l. The “threshold field” was then calculated and 
remaining outliers were removed (Figures 5 and 6). Values beyond the IOR of the 
difference between the sensor and the Winkler minus the mean difference between the 
sensor and the Winkler, were removed before calculation of the revised Soc values. For 
the primary sensor, 17 outliers were removed before calculation of the revised Soc 
(Events 24, 26, 28, 30,43, 48, 52, 56, 61, 71, 116, 118, 149, and 172 which 
correspond to sample ID's numbers 444701, 444726, 444744, 444752, 444771, 444783, 
444786, 444797, 444798, 444833, 444847, 444871, 444929, 445213, 445225, 445397, 
and 445529). Only one more threshold outlier for the secondary sensor was removed 
(Event 18, 444677) prior to calculating the new secondary sensor Soc value (Figure 6) 
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Table 6 shows the previous and revised Soc values and ratio for both the primary and 
secondary oxygen sensors (#1230 and #0345). 


The sensor values were then multiplied by their new corresponding Soc ratios to produce 
corrected primary and secondary sensor values. This correction brought both sensors 
closer to at 1:1 relationship with their respective Winkler replicate values (Figure 7). 
With the corrections applied the mean difference between the average difference between 
the primary and secondary sensor went from -0.0180 ml/l before correction to -0.0016 
ml/l after correction (Figure 8). 
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Figure 3. The difference between primary oxygen sensor #1230 and secondary oxygen 
sensor 70345. Outliers in red were removed prior to proceeding with Soc calculation: A) 
outliers from Events 1-15 (444601-444660), and B) Event 99 (CSL 04: 445093 - 
445102). The mean difference between sensors before outlier removal (solid blue line) is 
0.0168 ml/l. The upper and lower dotted blue lines are 0.0625 and -0.0277 ml/l 


respectively. 
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Figure 4. The mean difference (solid blue line) between 1* and 2" Winkler replicates (- 
0.002 ml/I). The lower and upper dotted blue lines are -0.17 and 0.16 ml/l respectively. 
Note the 7 outliers in red that were removed prior to proceeding with Soc calculation 
(sample ID numbers 444795, 444999, 444062, 445241, 445243, 445492, and 445547). 
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Figure 5. Outlier “threshold” values for the primary sensor were removed. The solid 
blue line is the mean value of the primary sensor threshold (~0.001 ml/l) and the lower 
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and upper dotted blue lines are -0.13 and 0.15 ml/l respectively. These outlier data points 
were removed and the remaining data were used to calculate the primary Soc values. 
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Figure 6. There was just a single "threshold" field value removed for the secondary 
sensor after the “bad” primary sensor threshold data had been removed. 
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Clo Secondary Oxygen Correction 
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Figure 7. The comparison between the corrected (blue) and uncorrected (black) A) 
primary (#1230) and B) secondary (#0345) sensors to the mean Winkler value. 


Table 6. Old and new Soc values for the primary and secondary SBE Oxygen sensors. 


po Old Soe Ratio (New:Old) 


Primary Sensor #1230 5.0347e-1 5.2597e-1 1.044693 
Secondary Sensor #0345 | 3.8281e-1 4.0107e-1 1.047702 
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Figure 8. Black dots — non-outlier differences between primary (71230) and secondary 
(#0133) sensor values before correction (black line is the mean = 0.0180 ml/l). Blue 


squares — Soc corrected difference between the primary and secondary sensor (blue line is 
the mean = 0.0016 ml/l). 
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Salinity 
(With portions extracted from HUD2014017 Cruise Report) 


Conductivity Calibration 


The salinometer outputs the conductivity as a ratio with the standard; therefore, some 
conversions are done to get the conductivity of the bottle. The standard has a given K15 
value: 


K15 = conductivity of standard seawater at 15°C and 1 atm/conductivity of KCI solution 
(32.4356g/kg) at 15°C and 1 atm. 


Where K15 = 0.99984 for this particular standard and the conductivity of KCI standard = 
4.29140 S/m and can be found in the seawater Matlab package (gsw C3515 function). 
Knowing K15 and the conductivity of the KCI solution, the conductivity of the standard 
seawater can be determined. Then, by multiplying by the conductivity ratio from the 
salinometer, the conductivity of the sample can be determined. 


It should be noted that these samples were analyzed with a bath temperature of 24°C 
rather than the 15°C that the standard conductivity was defined. The salinometer program 
accounted for this temperature difference so that the output sample conductivity ratios 
with the standard are at 15°C. 


Now we have the conductivity of the sample at 15°C and at the pressure of the bath in the 
salinometer; however, this needs to be converted to conductivity at the temperature and 
pressure of the CTD. This can be done using some functions from the same Matlab 
package (adopted for R using the Dan Kelley's oce package). 


First calculate the salinity of the bottle using the conductivity and pressure from the 
salinometer and a temperature of 15?C. 


Salinity bottle = gsw SP from C(Conductivity salinometer[mS/em[, TIC P bath) 


Then re-calculate the conductivity from this salinity value using temperature and pressure 
from the CTD. 


Conductivity bottle = gsw C from SP(Salinity bottle, T CTD,P CTD) %[mS/cm] 


This now gives conductivity values that can be compared to the CTD values. To correct 
the CTD conductivity a linear regression is done on this equation: 


Bottle conductivity — bI + b2*CTD conductivity 


to find an intercept, bl, and slope, b2, that will make the CTD conductivity better match 
the bottle conductivity. 


Figure 9 shows the difference between the primary (#3220 calibrated Dec 16, 2016) and 
secondary (#0864 calibrated Dec 15, 2016) sensors throughout the mission, filtered by 
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IQR to identify outliers. In Figure 10, the difference between the primary and the 
salinometer is examined and outliers are identified and removed using the IQR method. 
The mean difference between the primary sensor and the salinometer 1s -0.007 P.S.U. 
before outliers are removed, with an upper IQR threshold of -0.0098 and a lower 
threshold of -030306. All data points highlighted in red were removed before 
proceeding. 


Figure 11 compares the difference between the secondary sensor and the salinometer and 
identifies 3 additional outliers that were removed before coefficients were calculated. 
The mean difference between the secondary and salinometer was -3.45e-03 and the IQR 
upper limit is 1.05e-02 and the lower limit is -2.13e-02 (dotted blue lines). . 


The slope and intercept coefficients for both the primary and secondary sensors after are 
shown in Table 7. Figure 12 shows the difference between the 2 sensors both before and 
after correction. Before correction with new coefficients the average difference between 
filtered primary and secondary conductivity values was —3.36e-02 mS/cm. After 
correction, the average difference between sensors improved to 6.16e-05. 
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Figure 9. A) The mean sensor difference throughout the mission was was -0.0033 P.S.U 
(blue line). The lower and upper dotted blue lines are -0.007 and -0 ml/l respectively. 
Erroneous data (red dots) were removed before proceeding to the next step. 
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Outliers Outside 1,540QR 


Primary - Salinometer (P.S LL) 


0 200 400 600 Sub 


Ordered by Event and Increasing Sample D 


Figure 10. The difference between the primary sensor (73220) and the salinometer after 
the removal of erroneous sensor data. Erroneous values (red dots) were removed before 
proceeding. The mean difference (solid blue line) between the primary and salinometer 
was 7.00e-03 and the IQR upper limit 1s 9.80e-03 and the lower limit 1s -3.06e-02 (dotted 


blue lines). 
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Figure 11. The difference between the secondary sensor (70864) and the salinometer 
after removal of erroneous primary sensor and salinometer data. Only three additional 
erroneous values were removed prior to proceeding The mean difference (solid blue line) 
between the secondary and salinometer was -3.45e-03 and the IQR upper limit is 1.05e- 
02 and the lower limit 1s -2.13e-02 (dotted blue lines). 
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Figure 12. Before correction with new coefficients (black dots), the average difference 
between primary (#3220) and secondary (#0864) conductivity was —3.36e-03 mS/cm 
(solid black line). After correction (blue squares), the average difference between sensors 
was 6.16e-05 (solid blue line). 


Table 7. The revised intercept (b1) and slope (b2) terms calculated for both the primary 
(#3220) and secondary (#0864) conductivity sensors from EN2017606. 


Conductivity Sensor bl b2 
Primary (#3220) -1.0684e-02 . 1.000553 
Secondary (#0864) -7.7168e-03 . 1.000367 


Chlorophyll a 


Throughout the mission, ChlA was measured in-situ via a Wet Labs Eco-AFL/FL (SN: 
492 — calibrated Dec 15, 2016) attached to the CTD rosette (Appendix 3). Duplicate 
samples were regularly taken for ChlA analysis with a Turner Fluorometer from Niskin 
bottles fired in the upper 100 m. A comparison of the replicates showed that while the 
mean difference between replicates was -0.0061 ug/L, there were a total of 98 out of 577 
replicates that would be considered outliers (Figure 13). Outliers were selected via the 
1.5 * interquartile range (1.5 IQR) method discussed in the previous oxygen and salinity 
sections of this report. These outliers were removed before making the comparison 
between the WetLabs sensor values and the mean Turner replicate values (Figure 14). 
The relationship 1s confused and appears to be broken into 2 parts. There seems that 
there could be two separate relationships throughout the mission but the reason for this 1s 
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not entirely clear. The relationship is mostly above the 1:1 line from 0 to 0.4, above and 
below the line at 0.4 to 0.5 and mostly below the line after 0.5. Figure 15 shows the 
standardized percent difference between the sensor values and the Turner replicate mean 
throughout the mission. For roughly the first half of the mission (Halifax Line, Gully, 
Louisbourg Line, St. Anns Bank and Cabot Strait Line), the WetLabs fluorometer 
registered relative concentrations -40% greater than the Turner fluorometer. Over the 
second half of the mission, relative concentrations for both the Turner and WetLabs 
Fluorometer were roughly equivalent with a mean value of ~ -0.02 96. This suggests that 
the WetLab fluorometer was more in line with Turner readings for the second half of the 
mission and because the ship was in warmer nutrient rich waters to the southwest, the 
ChlA concentrations observed were generally greater. These two factors likely account 
for the unusual relationship observed in figure 14. 
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Figure 13. A total of 98/577 Turner fluorometer replicates were considered outliers 
using the IQR method. 
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Figure 14. The relationship between the WetLabs Fluorometer and the mean of the 
corresponding turner replicates. 
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Figure 15. The standardized percent difference between the fluorometer and the mean 
Turner fluorometer throughout the mission. 
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Water Samples for Chemical Analyses 


Station specific rosette bottle firing depths and water collections for chemical analysis 
can be found by referring to the CT? deck sheet binder and/or water chemistry sampling 
document prepared upon the conclusion of the mission and provided to ODIS. Table 5 


highlights CTD casts where water collections were made. 
pH Sensor 


The pH sensor (#1307, calibrated February 3, 2017) was deployed on the rosette only 
when the maximum depth was less than or equal to -1200 m. The CTD casts for which it 
was deployed are noted in Table 5. The sensor was included to support an ACCASP 
initiative investigating the delineation of ocean acidification and calcium carbonate 
saturation state of the Atlantic zone. 


Biological Program 


Narrative 


The “core” biological program conducted as part of cruise EN2017606, with some 
modifications, was a continuation of studies began in pre-AZMP years to describe the 
large-scale (spatial and temporal) variability in plankton biomass, productivity and 
biogenic carbon inventories on the Scotian Shelf. 


The program currently consists of essentially 2 elements: 


1. mesozooplankton community structure, population growth and biomass, and 
2. dissolved organic carbon measurements 


Table 5 provides a review of the stations where water samples were taken from rosette 
bottles for element 2 above. The mesoplankton sampling program is described below in 
more detail. This is followed by descriptions of “non-core” or ancillary biological 
sampling that includes text describing water sampling efforts in support of projects 
investigating: organic and organometallic micronutrients and their influence on primary 
productivity and phytoplankton community structure on the Scotian Shelf (Erin Bertrand 
— Dalhousie University), and water samples from strategic locations and depths to support 
a microbial community analysis via DNA, RNA and flow cytometry. The Biological 
Program section is concluded with a summary of pelagic seabird and marine mammal 
observations during EN606, provided by Carina Gjerdrum of the Canadian Wildlife 
Service. 


The ultimate aim of “core” studies is twofold: 


1. to provide a description of the inventories of biogenic carbon, their turnover rates and 
variability in space and time as part of Ocean Ecosystem Science Division's (OESD) 
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continuing climate studies, and 

2. to provide a description of plankton life-cycles and productivity on the Scotian Shelf 
and its influence or contribution to ecosystems in support of OESD's ecosystem- 
related research. 


Mesozooplankton Sampling 


Remarks/Comments 


In order to estimate the mesozooplankton community abundance and biomass, a conical 
ring net of 202 um mesh size with an aperture of 75 cm in diameter (filtering ratio 1:5) 
equipped with a KC Denmark flow-meter was towed vertically from the bottom to the 
surface at most stations (or from a maximum depth of 1000 m — AZMP standard). In 
total, there were 76 vertical ring net tows during the mission (Table 8, Figure 16). Of 
these, 1 was a 76 um mesh (30 cm diameter and 1:5 filtering ratio) at HL 02 (event 6). 
The 76 um net tow at HL 02 serves the same purpose of quantifying the community but 
targets a smaller fraction of the mesozooplankton community (1.e. smaller developmental 
stages, eggs and nauplii). Regardless of the mesh size, contents of the cod end were 
preserved in 4% buffered formaldehyde. 35 of the 202 um mesh tows were conducted at 
stations along core AZMP sections (HL, BBL, CSL and LL) (Table 8). The remaining 40 
200 um casts were conducted at ancillary stations throughout the mission (Table $, 
Figure 16). 


Six out of 76 casts were aborted for various reasons throughout the mission (HL 03.3, 
HL 06, HL 06.7, SG 28, STAB 05 and event 173 at HL 02). Of these, HL 06, HL 3.3 
and HL 02 were successfully reattempted. It should also be noted that BBL 07 was not 
occupied because of forecasted inclement weather. 


There were 6 - C3 genetics samples taken throughout the mission (of the 8 proposed 
locations) in support of Objective 13 (“Collect 200 um ring net zooplankton samples at 
predefined stations across the Scotian Shelf to supplement the Canada C3 program 
sample collection") (Table 8). 


Overall, net operations were successful during the mission. As with the CTD, net 
deployments occurred on the starboard side of the vessel, exposed to the wind and waves. 
This occasionally made operations technically difficult during inclement weather. 
Despite the challenges, the rate of unsuccessful tows was no more or less than typically 
experienced on our primary oceanographic platform (CCGS Hudson). 
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Figure 16. Locations for vertical ring net tows during EN2017606. Each tow is labelled 
with the consecutive mission event. 
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Microbial Protein and Organic Micronutrient Sampling 


Principle Investigator: Dr. Erin Bertrand (Dalhousie University, Department of 
Biology) 
Sampling by: Jenni Tolman and Ian Luddington (Dalhousie University) 


Objective 

To collect underway and rosette samples for protein and vitamin analyses in order to 
determine whether and how organic and organometallic micronutrients influence primary 
productivity and phytoplankton community structure on the Scotian Shelf. Sampling 
locations were coordinated with the LaRoche lab since our data types are synergistically 
informative. 


Microbial Protein Sampling 


Purpose 
Proteins are key to microbial activity: the type and amount of proteins present 


determines, in large part, the contributions microbes make to the ecosystems they occupy. 
Proteins can also be used as indices for nutritional status: elevated expression of specific 
proteins can be diagnostic for different nutritional states, such as nitrogen starvation, iron 
starvation, or vitamin starvation. Protein sequences also contain taxonomic information 
and can be used to assess contributions of different organisms to specific functions. 


samples were collected for targeted, mass spectrometry- based proteomic analyses of 
microbial communities in order to characterize the role of organic micronutrients in 
structuring phytoplankton communities on the Scotian Shelf. Primary objectives include 
measuring phytoplankton nutritional status indicator proteins (nitrogen, vitamin Bj, 
vitamin B, starvation) and vitamin- production biomarker proteins. Development and 
application of peptides for primary producer community composition analyses is a 
secondary focus. 


sampling Methods 
A total of 31 size- fractionated microbial protein samples (10L of water each) were taken 


from the CTD rosette at depths ranging from the surface to 250 m depth (Table 9) along 
the Halifax and Louisburg Lines, and in the Gully. In each case, water was pre-filtered 
(330 um) while dispensing from the Niskin bottle into 10L carboys. Water was then 
filtered through 3 and 0.2 um polycarbonate filters via peristaltic pumping. Filters were 
then frozen immediately at -80°C. 


Vitamin Sampling 

Purpose 

To determine the particulate and dissolved concentrations of organic and organometallic 
micronutrients on the Scotian Shelf. Organic and organometallic micronutrients are 
required by many phytoplankton groups and only produced by a select few microbes, 
setting up a series of interactive dependencies between microbial groups. The importance 
of these dependencies are not well known, as they have not yet been studied on the 
Scotian Shelf. Measuring the concentrations of these micronutrients in the particulate 
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and dissolved phases is one step towards understanding the role of microbial interactions 
in driving primary productivity and phytoplankton community structure. 


sampling Methods 
A total of 31 particulate and 23 dissolved vitamin samples (1L each) were taken from the 


CTD rosette at depths ranging from the surface to 250 m depth along the Halifax, Gully, 
and Louisburg lines (Table 9). Samples were protected from light and gently vacuum 
filtered through 0.2 um nylon filters. Filters were frozen at -80?C and dissolved samples 
were frozen in amber HDPE bottles at -20°C. 


Table 9. Protein and vitamin sampling, Bertrand Lab EN2017606. 


Station Event Depth ID# Protein RUE MEE 
(m) Vitamin Vitamin 
] 444635 - - - 
HL 02 7 20 444629 - - - 
7 40 444624 - - - 
80 444618 
l 444675 l I l 
HL 04 15 20 444670 l l l 
B 40 444666 l l l 
60 444662 1 l 1 
l 444721 l l l 
HL 06 24 20 444716 I I ] 
7 50 444712 l ] ] 
80 444707 l l ] 
l 444781 l l l 
HL 07 30 20 444777 l l l 
50 444772 I l ] 
l 444832 l l - 
GULD 04 E 20 444827 I I - 
_ 40 444823 I I 
60 444820 I ] - 
I 444870 I l l 
LL 09 T 20 444866 l l l 
B 80 444858 l l ] 
250 444854 l l ] 
I 444908 l l - 
LL 07 64 20 444903 1 ] - 
7 80 444896 l ] - 
250 444891 l 1 - 
LL 04 71 l 444946 l l l 
m 20 444940 l 1 1 
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40 444936 I l I 

80 444932 l l l 

l 444986 l l I 

LL 01 77 20 444981 l l l 
40 444977 l l l 

60 444973 l l l 


Microbial Community Analysis 


Principle Investigator: Dr. Julie LaRoche (Dalhousie University) 
Sampling by: Jenni Tolman and lan Luddington (Dalhousie University) 


Purpose 


Microbial communities and their associated processes are the foundation of marine life. 
Of particular interest to our group is the marine nitrogen cycle, comprising complex 
microbially-driven reactions whereby atmospheric nitrogen is fixed into a biologically- 
available form and cycled through the ecosystem. Though nitrogen is an essential 
element for life, the availability of fixed nitrogen can be a limiting factor for primary 
production and thus diazotrophs — organisms capable of biological nitrogen fixation — can 
be key to the productivity of an ecosystem. 


samples were collected for genomic and fluorescence-based analyses of the microbial 
communities on the Scotian shelf. Community composition will be assessed via 16S tag 
sequencing (bacteria and chloroplasts), and the naturally-fluorescent population will be 
characterized via flow cytometry. The latter method can also be used to quantify the 
bacterial community via nucleic acid stain SYBR green. Community function will be 
assessed via metagenomic sequencing, and qPCR assays for selected functional genes. 
Further samples were taken for manipulation in the lab, including targeted metagenomics 
and single cell isolation via fluorescence-associated cell sorting (FACS), and enrichment 
culturing of putative diazotrophs. 


sampling Methods 


Genomics: 


At 12 select stations along core AZMP lines (Halifax, Louisbourg, St Ann's Bank, and 
the Gully), duplicate 4L water samples were collected from the CTD rosette each of 4 
depths ranging from the surface to 300 m (Table 10). During collection, water was pre- 
filtered through a 330 um mesh to remove zooplankton. Each water sample was then 
sequentially filtered through 3 and 0.2 um polycarbonate filters by peristaltic pump until 
the water was depleted or the filters clogged. Filters were immediately frozen at -80 °C. 
samples have been collected at selected stations to provide time-series continuity with 
previous years (2014 and 2016). 
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Flow Cytometry: 


At each station and depth where genomic samples were collected, duplicate 2mL water 
samples (330um filtered) were fixed with 2% paraformaldehyde (PFA) for 10 minutes at 
room temperature, then frozen at -80°C for later enumeration of bacteria and 
characterization of the naturally fluorescent microbial community via the Accuri C6 flow 
cytometer. 


At select stations (Table 10), 45 ml of 330 um-filtered water were mixed with 5 ml of 


gly-TE buffer and frozen at -80 °C for later cell sorting on the BD Influx FACS 
instrument. 


Enrichment Cultures: 


At select stations (Table 10), large (1L) 330 um-filtered water samples were collected. 
These samples were spiked with phosphate (200 nM) and iron (2 nM) and stored in 
conditions approximating natural light/dark cycles and ambient temperature until return 
to the lab. 


Table 10. Microbial community samples, LaRoche lab EN2017606. 


Sorting 


auon Exent pow IDE c cam c Do ERO M: MU 
y y Cytometry 
] 444613 2 2 - - 
20 444609 2 2 - - 
HL 01 4 
B 40 444605 2 2 - - 
80 444603 2 2 - - 
] 440574 2 2 - - 
20 444630 2 2 - - 
HL 02 7 
a 40 444625 2 2 - - 
80 444619 2 2 - - 
] 444676 2 2 - - 
20 444671 2 2 - - 
HL 04 15 
m 40 444667 2 2 - - 
60 444663 2 2 - - 
] 444720 2 2 - - 
20 444717 2 2 - - 
HL 06 24 
B 50 444711 2 2 - - 
80 444708 2 2 - - 
HL 07 30 ] 444782 2 2 - - 


37 


000037 


I 444779 - - I I 
20 444776 2 2 - - 
50 444773 2 2 - - 
80 444769 2 2 - - 
l 444831 2 2 - - 
20 444826 2 2 - - 

GULD 04 50 
i 60 444819 2 2 = I 
250 444814 2 2 - - 
I 444869 2 2 I I 
20 444865 2 2 - - 

LL 09 56 
i 80 444859 2 2 š I 
250 444855 2 2 - - 
l 444909 2 2 - - 
20 444904 2 2 - - 

LL 07 64 
u 80 444897 2 2 - - 
250 444892 2 2 2 E 
I 444945 2 2 - - 
20 44494] 2 2 - - 

LL 04 71 
B 40 444937 2 2 - - 
80 444931 2 2 - - 
I 444985 2 2 - - 
20 444980 2 2 - - 

LL OI 77 
E 40 444976 2 2 - - 
60 444972 2 2 - - 
I 444998 2 2 - - 
10 444995 2 2 - - 

STAB 01 79 
E 20 444993 2 2 - - 
40 444990 2 2 - - 
l 445045 2 2 - - 
20 445041 2 2 à - 

STAB 05 88 
7 80 445035 2 2 - - 
300 445030 2 2 - - 
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Pelagic Seabird and Marine Mammal Observations 


Seabird Survey Report 
Leg1: 24 Nov — 4 Dec, 2017 
Canadian Wildlife ie Environment Canada 


Observer: Jeannine SES 


Background 


The east coast of Canada supports millions of breeding marine birds as well as migrants 
from the southern hemisphere and northeastern Atlantic. In 2005, the Canadian Wildlife 
Service (CWS) of Environment Canada initiated the Eastern Canada Seabirds at Sea 
(ECSAS) program with the goal of identifying and minimizing the impacts of human 
activities on birds in the marine environment. Since that time, a scientifically rigorous 
protocol for collecting data at sea and a sophisticated geodatabase have been developed, 
relationships with industry and DFO to support offshore seabird observers have been 
established, and over 100,000 km of ocean track have been surveyed by CWS-trained 
observers. These data are now being used to identify and address threats to birds in their 
marine environment. In addition, data are collected on marine mammals, sea turtles, 
sharks, and other marine organisms when they are encountered. 


Methods 


seabird surveys were conducted from the port side of the bridge of the Endeavor during 
the Scotian Shelf AZMP from 24 Nov to 4 Dec, 2017 (Leg 1). Surveys were conducted 
while the ship was moving at speeds greater than 4 knots, looking forward and scanning a 
90° arc to one side of the ship. All birds observed on the water within a 300m-wide 
transect were recorded, and we used the snapshot approach for flying birds (intermittent 
sampling based on the speed of the ship) to avoid overestimating abundance of birds 
flying in and out of transect. Distance sampling methods were incorporated to address 
the variation in bird detectability. Marine mammal observations were also recorded, 
although surveys were not specifically designed to detect marine mammals. Details of 
the methods used can be found in the CWS standardized protocol for pelagic seabird 
surveys from moving platforms’. 


Results 


Seabird sightings 
We surveyed 628 km of ocean from 24 Nov to 4 Dec, 2017. A total of 878 birds were 


observed in transect (1259 birds in total) from 7 families (Table 11). Burd densities 
averaged 4.5 birds/km/ (ranging from 0 — 88.5 birds/km^). The highest densities of birds 
(> 50 birds/km^) were observed on the Canso Bank and Western Banks (Figure 17A). 
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Dovekie accounted for 32% of the sightings (Table 11) and were scattered throughout the 
survey area (Figure 17B). The Scotian Shelf (and Grand Banks of NL) is an important 
wintering ground for Dovekie breeding in Greenland. Other Alcids observed in lower 
numbers include the Atlantic Puffin and Thick-billed Murre (Table 11). Northern fulmar 
and Black-legged Kittiwake were also relatively common (Table 11), especially in the 
deeper water (Northern Fulmar) and on the Canso Bank (Black-legged Kittiwake; Figure 
17C). A complete list of all species observed can be found in Table 11. 


Marine Mammal sightings 
A total of 36 marine mammals were recorded during the surveys (Table 12), 86% of 


which were long-finned pilot whales, observed in the eastern sections of the survey 
(Figure 17D). A single Grey Seal was also identified. 


Gully MPA 
Surveys were conducted within the Gully MPA in the afternoon of 30 Nov and the 


following morning on 1 Dec. A total of 46 birds were observed and 16 marine mammals 
in this area (Table 13; Figure 18). 


St. Anns Bank MPA 


Surveys were conducted within the St. Anns Bank MPA in the morning and early 
afternoon of 4 Dec before steaming to the mouth of Sydney Harbor to end Leg 1. A total 
of 74 birds and 11 marine mammals (all long-finned pilot whales) were sighted here 
(Table 14 and Figure 19). 
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Table 11. List of bird species observed during surveys on the Scotian Shelf AZMP, from 
24 Nov to 4 Dec, 2017. 


Number Total 
Family English Latin observed number 
in transect observed 


Phalacrocoracidae Unidentified Cormorant Phalacrocorax 16 16 


White-winged Scoter Melanitta fusca 0 2 
Anatidae Black Scoter Melanitta nigra l I 
Unidentified Duck All duck genera 0 7 


Dovekie Alle alle 278 552 
Atlantic Puffin Fratercula arctica 12 13 
Alcidae Thick-billed Murre Uria lomvia 10 10 
Unidentified Murre Uria 0 10 
Unidentified Alcid Alcidae 15 30 


Total 878 1259 
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Table 12. List of marine mammals observed during surveys on the Scotian Shelf AZMP, 
from 24 Nov to 4 Dec, 2017. 


English Latin Total number observed 


Unidentified Cetaceans Cetacea ] 


Unidentified Seals Phocidae 3 


Table 13. List of species observed in the Gully Marine Protected Area during surveys on 
the Scotian Shelf AZMP, from 24 Nov to 4 Dec, 2017. 


Species Latin Number observed in transect 


Herring Gull Larus argentatus 13 


Northern Fulmar Fulmarus glacialis 6 


Unidentified Skua Stercorarius 2 


Unidentified Cetaceans Cetacea 1 
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Table 14. List of species observed in the St. Anns Bank Marine Protected Area during 


surveys on the Scotian Shelf AZMP, from 24 Nov to 4 Dec, 2017. 


Number 
Species Latin observed in 
transect 


Dovekie Alle alle 23 


Black-legged Kittiwake Rissa tridactyla 9 


Great Black-backed Gull Larus marinus 3 


Total sightings 85 
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Figure 17. Density of A) all bird species combined, B) Dovekie, C) Northern Fulmar 
and Black-legged Kittiwake, and D) marine mammals observed during the seabird survey 
on the Scotian Shelf AZMP, from 24 Nov to 4 Dec, 2017. 
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Long-finned Bird Density 
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Figure 18. Density of birds and counts of marine mammals observed in the Gully 
Marine Protected Area on 30 Nov and 1 Dec, 2017. 
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Figure 19. Density of birds and counts of marine mammals observed in the St. Anns 
Bank Marine Protected Area on 4 Dec, 2017. 
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ARGO Float Deployments 


Contributions by: Ingrid Peterson 


Narrative 


There were a total of 6 APEX ARGO floats deployed during the mission (Figure 20 and 
Table 15). These floats continue to acquire data and their latest temperature profiles can 
be accessed on the following site by searching for their WMO numbers, 3901637- 
3901642 (Table 15). As of January 29^. 2018 the float profiles are not on the website but 
should be soon. 


hito://www.argodatamet.ore/Access-to-data/Description-ot-all-floats? 


AUN 


Q ARGO Foai 


“we ENGOS Mission Track 


> Protected Areas (includes Proposed 


AZUN 


Figure 20. The locations for each Argo float deployment during EN2017606. Refer to 
Table 15 for more details. 
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Mooring Operations 
Contributions by: Jay Barthelotte 


Narrative 


Over the duration of the mission there were 5 moorings recovered and 6 deployed (Figure 
21; Table 16). Please refer to Appendix 5 for the mooring diagrams. The Nova Scotia 
Current Mooring (M1996) was recovered on November 24^. and M2024 was deployed in 
its place on the same date and with the same sensor configuration. The first acoustic 
mooring (M1949) was recovered on November 25" from Emerald Basin and was not 
replaced. 


After completing the occupation of HL 07 late on November 28%, we began the steam 
towards Dawson Canyon to deploy acoustic mooring (M2027). We arrived on site just 
after midnight on the 29" and spent the next ~ 7 hrs doing 5 CTD casts (DC 01 — 
DC 04) in close proximity to the planned mooring deployment location. At day break, 
M2027 was deployed and we began the steam towards Logan Canyon to deploy M2028 
later in the afternoon of the 29%. This was followed by a CTD in close proximity 
(LC 01) before proceeding to the Gully MPA. After station SG 28 was occupied early 
in the morning of the 30%, acoustic release tests were conducted until early evening when 
M2026 was deployed in the Gully MPA. For the remainder of the 30% and overnight on 
December 1%, the remaining stations in the Gully were occupied before M1948 was 
recovered in the morning and then later replaced by M2025 at ~ the same location. 


After completing Gully operations, the ship sailed towards the deep end of the 
Lousibourg Line (LL 09). After occupying LL 09 in the early morning of the 2"4 the 
ship began the steam towards the Lophelia Conservation Area to recover M1950 by mid- 
morning of the 2". After this recovery, the rest of the Louisbourg Line and all of the St. 
Anns Bank Line was occupied before deploying the final acoustic mooring (M2029) on 
the morning of December 4^ Later on the same day, the final acoustic mooring 
(M1947) was recovered. 


It should be noted that the ADCP mooring (M1999) planned for recovery via dragging 
could not be contacted on the afternoon of December 4” and plans for dragging for the 
mooring were cancelled. After the mission we received information that parts of the 
mooring were discovered in Newfoundland on January 252018. 
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Figure 21. Mooring recovery and deployment locations during EN2017606. 


Table 16. List of mooring operations during EN2017606. The coordinates provided 
below are in decimal degrees and represents the ship's position at the time of the 
operation. 


Slat SLong 
Date Event Operation Station (DD) (DD) Program Comments 


24/11/2017 2 Deployment M2024 44.2455 -63.1631 NSCM Hebert 


27/11/2017 21 Release Test HL 06 


29/11/2017 40 Deployment M2028 43.4424 -59.9428 Acoustic 00" 


30/11/2017 45 Release Test SG 28 43.7567 -58.9669 


01/12/2017 53 Recovery M1948 43.8620 -589135 Acoustic Moors- 
Murphy 
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Moors- 


02/12/2017 59 Recovery M1950 44.4647 -57.1838 Acoustic 
Murphy 


04/12/2017 90 Recovery M1947 . 46.3562 -58.7306 Acoustic 


Murphy 


Underway Sampling 


Contributions by: Robert Benjamin’, Bill Fanning” 
! Program Coordination and Support Division, DFO 
^ Marine Technician V, Endeavor, University of Rhode Island 


Navigation 

Positional data and Date/time (SGPGGA and $GPZDA) from the ship's GPS was logged 
throughout the mission along with sounding data from the ships EK60 scientific echo 
sounder ($SDDBT). Heading data (SHEHDT) was also logged. These data were logged 
at 1 Hz throughout the mission using NavNet, a data logging and distribution system 
designed by NRCAN. Prior to the ship's return to BIO, navigation data was converted 
into daily coordinate logs at 1 second intervals in both .csv and .shp formats. 


The Endeavor's data logging systems where employed by the ships technician during the 
mission. This allowed logging of their TSG and ADCP systems internally and display of 
all systems where available during the mission including: ADCP, TSG, Wind direction 
and speed, Winch line-out and pressure, current position maps, and many on-board 
camera displays. A complete list of available sensors on the Endeavor can be found here: 
R:\Science\BIODataSvc\SRC\2010s\2017\EN606\Ship 

Deliverables\EN606 Hebert\scs\docs\sensor.html 


Underway Seawater System 


The Endeavor’s underway seawater system was used throughout the mission. The 
configuration file for the Thermosalinograph (TSG) on EN2017606 can be found in 
Appendix 6. 


Twenty-five gallons per minute is available from an intake located in the starboard sea 
chest, 48 feet from the bow. Seawater passes through a steel shut-off valve to a non- 
metallic pump. 1" PVC pipe to 1" PVC valves located in the Wet lab, 01 lab and on the 
01 deck supply a constant flow for devices such as incubators. The water in the Wet lab 
flows through a debubbler to a low-pressure manifold suitable for supplying flow through 
instruments. The flow through instrumentation included a SBE 21 SEACAT 
Thermosaliograph (TSG), and SBE3S remote thermistor located near the water intake, a 
WetLabs WetStart Fluorometer and WetLabs Eco-AFL/FL. 


50 


000050 


Prior to sailing the underway system was also plummed to include a water bath housing a 
ProOceanus CO;-Pro Atmosphere system to measure the partial pressure of CO». 


Every day, a single PCO2 and TIC sample, along with 2 ChlA samples were acquired and 
provided with a unique sample ID. The scanned paper log for these samples will 


eventually be located here: 
R:\Science\BIODataSvc\SRC\2010s\201 7\EN2017606\SCANNED LOGS and the 
digital e-logs can be found here: 


R:\Science\BIODataSvc\SRC\2010s\201 7\EN2017606\ELOG\Flow-Through Log. In 
total there were 21 PCO», 21 TIC and 42 Chla samples taken over this period. 


All underway sea-water system data was submitted to ODIS upon conclusion of the 
mission. Dr. Dave Hebert Dave. Fiebert{@dto-mpo.gc.ca) is the point of contact for these 
data. 


Other Underway Data 


The vessel also acquired a suite of underway measurements that are detailed on the 
Endeavor website. These data include vessel mounted ADCP, air temperature, humidity, 
wind speed and direction, barometric pressure, precipitation, short and long wave solar 
radiation and dual frequency (3.5 and 12 KHz) bathymetry. These data, as with all other 
data collected by ship provided equipment, were distributed to the DFO Data Manager 
and Chief Scientist upon the conclusion of the mission. They have been submitted to 
ODIS and can be found in the mission folder as specified in Appendix 7. 


Data Management 


Prepared by: Robert Benjamin 
Division: Program Coordination and Support Division, DFO 


Please refer to Appendix 7 for a table detailing the data collected during EN2017606. 


Data Collection 


In addition to standard AZMP manual data collection methods (i.e. various equipment 
specific deck sheets) ELOG, an electronic logbook system for collecting event metadata 
including position and sounding was used during EN2017606. This electronic logbook 
was accessible via computers connected to the RV Endeavor's network, including ship's 
data displays. Two locations in the main lab were used for data entry and one location in 
the Upper Lab for data Management. Metadata related to each piece of equipment was 
collected in the electronic log including position/time deployed, on bottom and recovered. 
Additional logbooks were employed to act as an itinerary, a daily operational log and a 
logbook to monitor the flow through. All digital logbooks were backed up daily and at 
the end of the mission were sent to ODIS for storage. After each event, the logbooks 
were entered into the Mission database. 
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CTD data was collected using the Endeavor's CTD system, setup and managed by the 
Endeavor data technician and backed up on the science server. After each CTD cast, the 
data was processed using CTDDAP and entered into the Mission database. 


Nav-Net, an on board ship's data collection system was used to send data to Elog. In 
addition, Regulus was also used to record ship's data sent to the science team during the 
entire mission. These data will be located in the archive here: 


R: Science BlODataSve SRC 201082201 1EN606 Nav 


At the end of the mission, the Endeavors’ data technician supplied a drive which 
contained all data collected by the vessel during the mission including TSG data. These 
data can be found here: R:\Science\BIODataSvc\SRC\2010s\2017\EN606\Ship 
Deliverables 


NOTE: pC02 data was collected from the TSG system during the mission but NOT stored 
in SCS. These data will be stored here: 


R:\Science\BIODataSvc\SRC\2010s\2017\EN606\pCo2 


Salinity, Winkler Oxygen and Chlorophyll was analyzed while at sea. Data from the 
Analysis was routinely backed up and entered into the mission database. 


Data Input Template 


Reports were generated from shipboard input data in the AZMP Template Database to 
compare with corresponding CTD sensor data and conduct preliminary analyses included 
in this report. 


GIS 


Daily navigation and operations were maintained in a graphical information system 
(QGIS). Final line and point shapefile were generated from these data for the cruise 
report. 


Hardware 


One laptop was used to run the NavNet software. GPS data and Sounding data was sent 
to this computer via serial RS232 and logged. Data was transferred to our other 
computers via the ships network. The Endeavor’s TSG system was used during the 
mission and data was saved in the ships SCS repository. pCo2 was collected in the Wet 
Lab along with the TSG. 
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APPENDICES 


Appendix I. Gully and St. Anns Bank MPA Activity Approvals 


Gully Approval 
Signed Letter from RI 


StAnnsMPA EN2017 
606 approval. pdf 
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Document Released Under the Access to 
Information Act / Document divulgué en vertu 


ef Fisheries and Oceans Pêches et OcéÁ £ la Loi sur l'accès à l'information. 


Canada Canada 

PO Box 1006 File / Référence 
Dartmouth, NS GMP A-201 5-05 
B2Y 4A2 


Mr. Andrew Cogswell 

Ocean and Ecosystem Sciences Division 
Bedford Institute of Oceanography 

P.O. Box 1006 

Dartmouth, NS B2Y 4A2 


Dear Mr. Cogswell: 


RE: Gully Marine Protected Area (MPA) Activity Approval 


I am pleased to inform you that your request to collect seawater chemistry and zooplankton 
samples in the Gully MPA during Fall and Spring cruises of the Atlantic Zone Monitoring 
Program through October 2018 has been approved under Section 6(1) of the Gully MPA 
Regulations. Information provided in your 30 June application demonstrates compliance with 
the regulatory conditions that must be satisfied for Ministerial Approval. Any changes to the 
approved activities that have not been described in the submitted plan must be discussed with 
the Oceans and Coastal Management Division (OCMD) prior to commencement. 


While in the MPA, you will be expected to comply with all applicable federal legislation. In 
particular, we'd like to emphasize that holding a Ministerial Approval issued pursuant to the 
Gully MPA Regulations and the Oceans Act does not satisfy the requirements of the Species at 
Risk Act or the Fisheries Act. Neither does the MPA Approval given here substitute for any 
permits or licences required under those statutes. It is your responsibility to ensure any 
necessary authorizations are acquired prior to undertaking the approved MPA activities. 
Further, you are expected to comply fully with any conditions stipulated therein. 


To support conservation and protection of the MPA ecosystem, you are asked to adhere to the 
following requests while undertaking the approved research: 


l. Maintain a watch during daylight hours for marine mammals as described in the submitted 
Activity Plan. Record those observations and any sightings of turtles and marine debris 


(e.g., abandoned fishing gear, plastics, other garbage or pollutants). Provide sightings 
information to Oceans and Coastal Management Division (OCMD). 
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Document Released Under the Access to 
Information Act / Document divulgué en vertu 
de la Loi sur l'accès à l'information. 


2. Report environmental emergencies or other incidents, including unintentional discharges 
and mammal collisions, to the Canadian Coast Guard within two hours of the occurrence. 
Notify OCMD as soon as possible and file an incident report. 


3. Provide an activity report to OCMD that details MPA arrival and departure dates & times, 
and outlines operations undertaken within the MPA. 


We have enclosed a set of templates and instructions to assist with the documentation being 
sought in the requests attached to this Approval. The activity report, sightings data, and 
incident notifications should be submitted to Tanya Koropatnick in OCMD. Please feel free to 
follow-up with Tanya should you have any questions or need further clarification: 


Tanya Koropatnick 

Fisheries and Oceans Canada 
Bedford Institute of Oceanography 
1 Challenger Dr., B500 
Dartmouth, NS 

B2Y 2V9 


Phone: (902) 426-6048 
E-mail : Tanya.Koropatnick(g)dfo-mpo.gc.ca 


Thank you for your continued interest in the Gully MPA. The bi-annual AZMP sampling 
program contributes data to the ecological knowledge base and is directly linked to the site's 
long-term monitoring needs. Your AZMP trips also provide a consistent and predictable 
platform for the conduct of other important research and monitoring activities. OCMD remains 
grateful for this ongoing Science Branch support and welcomes the AZMP for another 3 years 
of sampling in the Gully MPA. 


Yours sincerely, 


Z, Moda 


Glen Herbert 

A/Regional Manager 

Oceans and Coastal Management Division 
Ecosystem Management Branch 
Maritimes Region 


Attachments: Approved Activity Plan 
Activity and Incident Report Template 
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Fisheries and Oceans Péches et Océans 


Canada Canada 

PO Box 1006 File / Référence 
Dartmouth, NS SABMPA-2017-10 
B2Y 4A2 


Andrew Cogswell 

Bedford Institute of Oceanography 
1 Challenger Dr 

Dartmouth, NS 

B2Y 4A2 


Dear Mr. Cogswell: 
RE: St. Anns Bank Marine Protected Area (MPA) Activity Approvals 


I am pleased to inform you that your request to conduct the Atlantic Zone Monitoring 
Program (AZMP) in the St. Anns Bank Marine Protected Area (MPA) has been approved 
under section 10(1) of the St. Anns Bank MPA Regulations. Based on the information you 
submitted on October 5, 2017, Fisheries and Oceans Canada (DFO) has determined that your 
activity meets the regulatory conditions required for issuance of Ministerial Approval. 


We will anticipate your advance notices before sailing under any vessel used. Any changes to 
the approved activities that have not been described in the information submitted must be 
discussed with the Oceans and Coastal Management Division (OCMD) prior to 


commencement. Approval renewals for future years are subject to reassessment. 


While in the MPA, you will be expected to comply with all applicable legislation. In 
particular, we would like to emphasize that holding a Ministerial Approval issued pursuant to 
the St. Anns Bank MPA Regulations under the Oceans Act does not satisfy the requirements 
of the Species at Risk Act or the Fisheries Act. Neither does the MPA Approval given here 
substitute for any permits or licences required under those statutes. It is your responsibility to 
ensure any necessary authorizations are acquired prior to undertaking the approved MPA 
activities. 


As per Section 11 of the St. Anns Bank MPA Regulations, a post-activity report on the 
activity conducted in the MPA is required 90 days following the completion of the activity. 
We will contact you shortly with more details on what is required for reporting. 


Vessels must comply with all relevant provisions of the Marine Mammal Regulations 
pursuant to the Fisheries Act. Further guidance for vessel operation can be found in Section 
A2 of the Annual Notices to Mariners (Marine Mammal Guidelines and Marine Protected 
Áreas) as published by the Canadian Coast Guard (https://www.notmar.gc.ca). 


To support conservation and protection of the MPA ecosystem, you are asked to adhere to 
the following requests while undertaking the approved research: 


1. Maintain a watch during daylight hours for turtles, marine mammals and marine debris 
(e.g., abandoned fishing gear, plastics, other garbage or pollutants). 


coll 
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2. Report any marine mammal collisions, entanglements, distressed or dead animals to the 
Marine Animal Response Society's emergency hotline (1-866-567-6277), or via VHF 
channel 16. Sightings of healthy marine mammals should be reported to 
XMARwhalesightings@dfo-mpo.ge.ca. The following information about the sightings 
should be included: date, time, location and species. Photos and videos should be 
submitted if available. 


3. Report all live and dead sea turtle sightings and incidents (e.g., entanglements, collisions) 
to the Canadian Sea Turtle Network's hotline (1-888-729-4667) or online at 
http://seaturtle.ca/turtle-sighting/. The following information about the sighting or 
incident should be included: date, time, location, species, and condition of the animal. 
Photos and videos should be submitted if available. 


4. Report environmental emergencies or other incidents, such as unintentional discharges, to 
the Canadian Coast Guard within two hours of the occurrence. Notify OCMD as soon as 
possible and file an incident report. 


Once completed, the activity report and incident notifications should be submitted to Melanie 
MacLean in OCMD at: 


Melanie MacLean 

Fisheries and Oceans Canada 

Bedford Institute of Oceanography 

1 Challenger Dr., B500 

Dartmouth, NS 

B2Y 2V9 

Phone: (902) 440-5727 

E-mail: melanie.maclean@dfo-mpo.ge.ca 


Please confirm with Melanie MacLean that you have received this letter, and feel free to 
follow-up anytime should you have any questions or need further clarification. 


Yours sincerely, — 


S. Ss 


Wendy Williams 

A/Director, Oceans Management 

Oceans and Coastal Management Division 
Ecosystem Management Branch 

Fisheries and Oceans Canada, Maritimes 
Region 
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Appendix 2. Crew List for the R/V Endeavor 606. 
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Appendix 3. CTD Configuration File — EN606.xmlcon 


Date: 01/17/2018 


Instrument configuration file: 
R:\Science\BIODataSvc\SRC\2010s\2017\EN606\CTD\CTD PROCESSING\EN606\EN 
606.XMLCON 


Configuration report for SBE 911plus/917plus CTD 


Frequency channels suppressed : 0 
Voltage words suppressed  :0 


Computer interface : RS-232C 

Deck unit : SBE1 plus Firmware Version >= 5.0 
Scans to average : 1 

NMEA position data adde : Yes 

NMEA depth data added : No 

NMEA time added : No 


NMEA device connected to : deck unit 
Surface PAR voltage added  : Yes 
Scan time added : NO 


1) Frequency 0, Temperature 


Serial number : 2902 
Calibrated on : 15-Dec-16 


G : 4.34451712e-003 
H : 6.447303 10e-004 
I : 2.28889365e-005 

J : 2.12526223e-006 

FO : 1000.000 

Slope : 1.00000000 


Offset : 0.0000 
2) Frequency 1, Conductivity 


Serial number : 3220 
Calibrated on : 16-Dec-16 


G :-0,77555876e+000 
H : 1.34416455e+000 

I : -2.86467321e-005 

J : 6.94809804e-005 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
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3) Frequency 2, Pressure, Digiquartz with TC 


Serial number : 0444 
Calibrated on : 20-Dec-16 


Cl : -5.378517e+004 
C2 : -3.498580e-001 

C3 : 1.648580e-002 

DI : 4.036100e-002 

D2 : 0.000000e+000 

TI : 2.984744e+001 

12 : -3.538190e-004 

T3 : 3.972770e-006 

T4 : 2.922330e-009 

T5 : 0.000000e+000 
Slope : 0.99989692 
Offset : -0.45761 
AD590M : 1.125800e-002 
AD590B : -8.763490e+000 


4) Frequency 3, Temperature, 2 


Serial number : 2034 
Calibrated on : 13-Dec-16 


G : 4.41249522e-003 
H : 6.41293978e-004 
I : 2.37750205e-005 

J : 2.28693904e-006 

FO : 1000.000 

Slope : 1.00000000 


Offset : 0.0000 
5) Frequency 4, Conductivity, 2 


Serial number : 0864 
Calibrated on : 15-Dec-16 


G : -3.93005749e+000 
H : 5.65787779e-001 

I :-6.14331081e-004 

J : 6.37838626e-005 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
6) A/D voltage 0, Transmissometer, WET Labs C-Star 


Serial number : 969DR 
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Calibrated on : 06-Dec-16/15-Feb-17field 
M : 19.5917 

B : -1.1363 

Path length : 0.250 


7) A/D voltage 1, Fluorometer, WET Labs ECO-AFL/FL 


Serial number : 492 

Calibrated on : 15-Dec-16 

Dark output : 0.0250 

Scale factor : 2.40000000e+001 


8) A/D voltage 2, Altimeter 
Serial number : 49899 
Calibrated on : 30-Mar-15 
Scale factor : 15.000 
Offset : 0.000 
9) A/D voltage 3, PAR/Irradiance, Biospherical/Licor 


Serial number : 70513 
Calibrated on : 21-Nov-16 


M : 1.00000000 

B : 0.00000000 

Calibration constant : 9900990099.00989910 
Multiplier : 1.00000000 

Offset : -0.10245222 


10) A/D voltage 4, Oxygen, SBE 43 


serial number : 1230 
Calibrated on : 02-Aug-17 
Equation : Sea-Bird 


Soc : 5.03470e-001 
Offset : -5.15300e-001 
A : -3.72370e-003 
B : 2.04260e-004 
C : -2.88240e-006 
E : 3.60000e-002 
Tau20 : 1.81000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
HI : -3.30000e-002 
H2 : 5.00000e+003 
H3 : 1.45000e+003 


11) A/D voltage 5, Oxygen, SBE 43, 2 
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serial number : 0345 
Calibrated on : 02-Aug-17 
Equation : Sea-Bird 


Soc : 3.82810e-001 
Offset : -7.22200e-001 
A : -4.10370e-003 
B : 1.65940e-004 
C : -2.39230e-006 
E : 3.60000e-002 
Tau20 : 1.24000e+000 
DI : 1.92634e-004 
D2 : -4.64803e-002 
HI : -3.30000e-002 
H2 : 5.00000e+003 
H3 : 1.45000e+003 


12) A/D voltage 6, pH 


Serial number : 1307 
Calibrated on : 03-Feb-17 
pH slope  : 4.6258 
pH offset : 2.5383 


13) A/D voltage 7, Fluorometer, WET Labs ECO CDOM 
Serial number : 3745 
Calibrated on : 23-Nov-2017 
Dark output : 0.000 
Scale factor : 3.000 
14) SPAR voltage, Unavailable 
15) SPAR voltage, SPAR/Surface Irradiance 
Serial number :20190 
Calibrated on : 21-Nov-16 
Conversion factor : 1565.10000000 
Ratio multiplier : 1.00000000 


Scan length : 40 
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Appendix 4. Preliminary Section Plots and Anomalies (T/S/Sigma-T) 


Section plots were produced for Temperature, Salinity and Sigma-T all sections from 
EN606 (See map below). It should be noted that no anomalies were produced because 
this mission was well outside of the typical sailing dates for most of the preceding fall 
AZMP missions. Finally, BBL 07 was not occupied during the mission as noted in the 
mission narrative. 


Endeavor EN606, November 24 - December 16, 2017, 79 Stations 
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St. Anns Bank Line 
Section 
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Cabot Strait Line 
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Browns Bank Line 
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Appendix 6. Endeavor TSG Configuration File — 27Nov201 7a.xmlcon 


Date: 01/29/2018 


Instrument configuration file: R:\Science\BIODataSvc\SRC\2010s\2017\EN606\Ship 
Deliverables\EN606 Heberttsg raw 27Nov2017a.X MLCON 


Configuration report for SBE 21 Seacat Thermosalinograph 


Remote temperature : SBE 3 
External voltage channels : 2 
sample interval : 6 seconds 


NMEA position data added : Yes 
NMEA depth data added : No 


NMEA time added : NO 
NMEA device connected to : deck unit 
Scan time added : No 


1) Frequency 0, Temperature 


Serial number : 1578 
Calibrated on : 14-Dec-16 


G : 4.19581328e-003 
H : 5.93731371e-004 
I : 3.79065958e-006 
J : -1.86524830e-006 
FO : 1000.000 

Slope : 1.00000000 


Offset : 0.0000 
2) Frequency 1, Conductivity 


Serial number : 1578 
Calibrated on : 14-Dec-16 


G : -4.01390242e+000 
H : 4.78841495e-001 

I : 1.20442267e-003 

J : -2.89460350e-005 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
3) Frequency 2, Temperature, 2 


Serial number : 0604 
Calibrated on : 15-Dec-16 
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G : 4.80105098e-003 
H : 7.12509078e-004 
I : 4.85485321 e-005 

J : 6.14112724e-006 

FO : 1000.000 

Slope : 1.00000000 


Offset : 0.0000 
4) A/D voltage 0, Fluorometer, WET Labs WETstar 
Serial number : 1177 
Calibrated on : 17-Mar-2017 
Blank output : 0.063 
Scale factor : 6.100 
5) A/D voltage 1, Fluorometer, WET Labs ECO-AFL/FL 
Serial number : 478 
Calibrated on : 12142016 
Dark output : 0.0170 
Scale factor : 2.40000000e+001 


Scan length : 34 
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Appendix 7. Data and Meta-data Collections 


The mission data 1s uniquely organized because the charter vessel provided us with data 
files upon our departure for all shipboard instrumentation. 
processed using the Endeavor's protocols but was also processed using CTD-Dap to meet 


AZMP Maritimes standards. 


The mission data and metadata 1s held here: 


R: Science BlODataSve SRC 201082201 7EN606 


This folder includes: 

Raw and processed CTD data, configuration files and plots 

Lists of stations and Navigation 

Logs as they are scanned 

Raw shipboard analysis (Winkler, Autosal, Turner Fluorometer) 
Operation metadata (Elog) 

The AZMP database template for the mission and summary reports 
Ship deliverables 


l. 


UN UA m 


HPRP 


ADCP 

CTD raw data and Endeavor processing 
Navigation 

SCS log 

TSG data and configuration files 
Winch logs 


The BioChem folder will contain land based laboratory analysis as it becomes 
available, and includes: 


ES 


HPLC/Absorption 

POC/PON 

Nutrients 

Zooplankton 

Flow cytometry (samples collected but may be late to process) 
PCO2 

TIC/TA 
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CRUISE NARRATIVE 


Highlights 


Ave Dec een NAFO Regions: 5Ze, 5Y, 4X, 4W, 4Vs, AVn, 3Pn, 3Ps 
š | Extent: 41° 46'N - 47° 35'N; 056° 02'W - 070? 16'W 


Expedition Designation: HUD2018004 or 18HU18004 CSDM format) 


Andrew Cogswell 

Ocean Ecosystem Science Division 
Marine Ecosystem Section 
Department of Fisheries and Oceans 
Bedford Institute of Oceanography 
PO Box 1006 

Dartmouth, NS, Canada B2Y 4A2 
Andrew.Cogswell@dfo-mpo.gc.ca 


Chief Scientist: 


Ship: CCGS Hudson (call sign - CGDG) 
i oceanographic research vessel 

April 6, 2018 — Depart BIO, Dartmouth, NS 

Apr 19%. 2018 — Emergency ship repair, Sydney, NS 

Apr 20°, 2018 — Depart, Sydney, NS 

Apr 24", 2018 — Arrive Mulgrave, Dartmouth, NS 


Ports of Call: 


Cruise Dates: Apr 6” — 24" 


Mission Summary 


Overview 


One day prior to departure on April 5", while conducting a pull test on the CTD winch, we 
realized the joy stick controls were malfunctioning. When the winch was set to the prescribed 
low speed required for producing enough torque to recover and deploy the CTD, the winch payed 
out regardless of the joy stick position. All testing on the winch to that point had been conducted 
in high speed and was functioning normally. Pennicon was called to service the winch and they 
worked with Hawboldt to diagnose the issue. Just after 1800 on the 5% the issue had been 
resolved and a re-termination of the CTD cable was completed by 2300 that evening. A special 
thanks to the ship's crew, Pennicon, Hawboldt, DFO Field Ops and Terry Cormier and Dave 
Levy for their extra effort to get the winch functioning properly the night before sailing. 
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Just after departure from the BIO Jetty at -9 am on April 6". the CCGS Hudson conducted a 
compass swing, boat and fire drills and testing of the starboard life boat prior to science 
operations beginning in the basin at ~1300. The underway system in the forward lab was started 
just after pulling away from the dock. The basin tests consisted of 3 CTD deployments to 
diagnose an issue with the altimeter and the altimeter cable. As well, the hydro winch was tested 
to make sure it was functioning properly. The ship began steaming towards BBL 01 at -1600, 
but off Chebucto Head the shaft temperature was increasing and an oil change was required that 
resulted in a short ~1-2 hour delay before we could begin steaming again. We arrived at BBL 01 
at ~0100 on April 8% to begin station occupations. At ~1830 on the 9^ we completed the BBL 
and PS lines and began the Plymouth Line at PL 09. At 1200 on the 9", the ship received word 
from Global Affairs Canada that the US State Department had provided our clearance to conduct 
work in US waters. Early in the morning of the 10 we began occupying stations along the PL 
line, crossing into US waters at ~1000 on April 10%. The ship continued down the PL line, but 
there was concern from the bridge that an evening arrival at the nearshore YL 10 station could 
make it difficult to navigate around fishing gear in the area. Upon completion of PL 04, the 
Captain suggested we occopy stations YL 06 and YL 07 during the night and early morning of 
April 11% and wait for day break to return to PL 03. We completed the PL line at ~1020 on the 
11% and began occupying the Yarmouth Line (YL 10) at -1200 on April the 11%. The ship 
worked its way down the line, crossing back into Canadian waters at -0230 on the 12%, and 
finishing the line ~8 hours later before proceeding to HL 01. 


Occupation of the Halifax Line stations began at -0100 on April 13%. A station occupation was 
added to the plan at HL 02.25 in close proximity (~2 km) to a DFO glider that was in the area. 
Upon completion, operations proceeded as normal until ~ 1800 on April 14 when the CTD deck 
unit threw an error during the bottom of the cast at HL 07 in ~2800 m of water. This led to a 
string of diagnostics and 7 coinciding CTD test deployments that finally culminated in a 
successful full profile at ~1030 on April 15". The resulting ~17 hrs of lost time meant that 
HL 12, 13 and 14 were cut from the mission plan. Approximately 24 hours later on April 16% at 
1000, the Halifax Line was finished at HL. 11 and the Hudson departed for the Gully MPA. The 
finer details associated with the CTD repairs, and all other aspects of the CTD operations 
throughout the mission, are discussed at length in the CTD summary of this report. 


The Hudson arrived in the Gully at station SG 28 at -2200 on the 16", followed by typical ring 
net and CTD occupations at GULD 03 and GULD 04 overnight and early morning of the 17 
At SG 23 the ring net was deployed, but the continued to deteriorate during deployment. The 
forecast suggested that the winds and associated sea state would be beyond safe working 
conditions for at least the next 12 hours, so the CTD at this location was cancelled and a course 
was set for station 9 on the South Eastern St. Pierre Bank Line (SESPBL - A line added to the 
plan to satisfy the core requirements of the NL AZMP program that could not occupy the line due 
to weather) During the early part of the transit to this location under heavy seas, the Senior 
Engineer's Cabin experienced flooding, due to a leaky starboard scupper that could not be 
accessed under a container. The decision was made to halt the transit and begin the steam back 
in to Sydney to make the necessary repairs. On the way back to Sydney, the Louisbourg Line 
was occupied starting at LL 07 at -2000 on April 175 LL 08 and LL 09 were dropped because 
of the sea state and wind direction (more on this later). The goal was to have the ship alongside 
in Sydney on the morning of April 19%. We were able to occupy LL 07 to LL 01 and STAB 01 
to STAB 04 before heading for Sydney at -2200 on April 18% 
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We arrived at the dock in Sydney on the morning of the 19% at 0830. The underway system was 
shut off upon docking and the container was removed from the ship by -0930. We were 
instructed by the Captain to await an update later in the day about the state of the repairs and a 
possible departure time. By mid-afternoon on the 19%, we received notice that the repairs were 
complete, the container would go back on board that evening and the ship would set sail at ~0800 
on the 20", bound for CSL 01. A cruise track was provided to the bridge based on the amount of 
remaining time and the likely speed of the vessel given the long range wind forecast. Prior to 
leaving the Harbour, the underway system was turned back on. 


Station occupations began at CSL 01 at ~1400 on April 20", and finished ~15 hours later at 0500 
on the 21* before proceeding to STAB 05. STAB 06 was dropped from the itinerary because it 
is not within the bounds of the MPA and was only recently added to the AZMP sampling suite. 
After completing the occupation of STAB 05 at 1330 on April 21*, the Hudson began the long 
steam towards LL 08, arriving at -0330 on the 22" Work at LL 08 and LL 09 was completed 
at ~1300 on April 22" before steaming for BP 00 on the eastern shelf at the mouth of the 
Laurentian Channel. Work began at -2200 on April Dral BP 00 and continued through the 
night and morning of the 23", finishing at BANQ B1 at ~1300 before setting sail for the crew 
change in Mulgrave. Just after departure from Sydney on the 20%, we were informed that the 
berth would not be available for the Hudson, so an alternate port was chosen for the exchange. 
We arrived at the Mulgrave dock early in the morning of the 24”. 


Over the 18 day mission, the CCGS Hudson logged ~2461 nm and AZMP science staff 
conducted 180 separate operations at 85 stations (Figure 1). Table 1 breaks down the operations 
by sampling gear for each leg of the trip. The table also points to figures that display the 
deployment locations for each gear type. Each of these figures 1s accompanied by a table of 
coordinates detailing each deployment of that gear type. Table 2 contains the break down in time 
allocated to each gear type. 
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5 — HUDSOT8004 Operations 


Figure 1. The locations for all 180 events during the HUD2018004 AZMP spring survey. Some 
overlapping event labels may not visible. 


Table 1. Summary of operations during the HUD2018004 AZMP spring survey. 


Operation H of Operations Figure 
CTD 95 4 
Vertical Ring Net Tows 85 29 


Table 2. Break down of operational time by gear type during HUD2018004. 


Gear Operation Duration (hrs) 
CTD ~66 
Vertical Net Tows ~27 


Mission Participants 


A complete crew list for this mission can be found in Appendix 2. Please refer to Table 3 for a 
list of Science Staff who participated during HUD201 8004. 


000092 


s.19(1) 


Table 3. List of science staff aboard the HUD2018004 fall AZMP mission. 


Name Affiliation Duty Shift 
1 | Benjamin, Robert DFO — PCSD Data Manager Day 
2 Cogswell, Andrew ** DFO — OESD Chief Scientist Day 
3 | Cormier, Terry DFO — OESD CTD Technician Night 
4 Waclawik, ee Dal — Bertrand/LaRoche Technician Split 
5 | Dal — Bertrand/LaRoche Technician Split 
6 Wilson, Colleen Dal — Hays CTD Operator/Elog Day 
7 Layton, Chantelle DFO — OESD CTD Operator/Elog Night 
8 Levy, Dave DFO — OESD CTD Technician Night 
9  Maclsaac, Kevin DFO — OESD CTD/Nets/Biologist Night 
10 Perry, Tim DFO — OESD Lab Technician Night 
11 Spry, Jeffrey DFO — OESD CTD/Nets/Biologist Day 
12  Thamer, Peter DFO — OESD Lab Technician Day 
13 Winkel, Jeannine ECCC — CWS Bird/Mammal Day 


**Chief Scientist 

DFO: Department of Fisheries and Oceans Canada 

OESD: Maritimes - Ocean Ecosystem Science Division 

PCSD: Maritimes - Program Coordination and Support Division 
EC-CWS: Environment Canada - Canadian Wildlife Service 
DAL: Dalhousie University 


Objectives 


There were 12 defined objectives in the final version of the Form B submitted to Coast Guard 
Headquarters on Aprilt 4% 2018 (below). Of these, 2 objectives were cancelled prior to sailing 
(5 and 12). It was decided that the AZOMP mission would deploy the ARGO floats intended for 
the AZMP mission and a small vessel (Sigma T) was used to deploy the ALSEAMAR glider. 
Table 4 describes whether each of these objectives was met along with any relevant supporting 
commentary. 


Primary 


1. Obtain observations of the hydrography and distribution of nutrients, phytoplankton and 
zooplankton at standard sampling stations along “core” Atlantic Zone Monitoring 
Program sections within the Maritimes Region (Contact Mr. Andrew Cogswell - 
http//www.bio.gc.ca/science/monitoring-monitorage/azmp-pmza-eng.php.). 


Additional 


2. Occupy stations in support of the extended Halifax Line (XHL) (HL 08 and greater) 
(Contact Dr. Igor Yashayaev). 

3. Carry out hydrographic, chemical and biological sampling in support of Gully MPA 
monitoring initiatives by Oceans and Coastal Management Division (Contact Dr. Dave 
Hebert - hitp:;//inter-w02.dto-mpo.gc.ca/Maritimes/Oceans/OCMD/Guliv/Guliv-MPA ). 

4. Nutrients and hydrography across the Northeast Channel and Gulf of Maine as part of 
NERACOOS Cooperative Agreement, (Contact Dr. Dave Hebert - 

http://www. neracoos.org/). 
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5. Deploy 3 ARGO floats in support of the International Argo Float Program (Contact Dr. 
Blair Greenan - hitp:///www.meds-sdmm.dfo-mpo.gc.cea/isdm-gdsVargo/index-eng.html). 

6. Collect underway and CTD water samples at specified locations and depths to fulfil the 

regional component of an Aquatic Climate Change Adaptation Services Program 

(ACCASP) initiative investigating the delineation of ocean acidification and calcium 

carponate sara non stale of nd — zone — Dr. Kumiko EU COH - 

http//www.dfo-mpo.gc.ca/science/oceanography-oceanographie/accasp-psaccma/index- 
IE 

7. Collect water samples for the Bertrand lab at Dalhousie University to evaluate whether 
and how organic and organometallic micronutrients influence primary productivity and 
phytoplankton community structure on the Scotian Shelf (Contact Erin Bertrand — 
Erin. Bertrandi dal.ca). 

8. Collect water samples from strategic locations and depths to support a mucrobial 
community analysis via DNA, RNA and flow cytometry, as well as the isolation of novel 
diazotrophs (Contact Dr. Julie Laroche - 

http //www.dal.ca/facultv/science/biologv/facultyv-staff/our-facultv/1ulie-Jaroche.html). 

9. Bird and mammal observations as part of EC-CWS sea-bird observation program and in 
fulfilment of Gully ang zt Pins Bank MPA occupation requirements (Contact Carina 
Gjerdrum — carina.gjerdr da.c 

10. Carry out Be mee cg fen cl and led sampling at stations in the St. Anns 
Bank MPA as continued support for the Oceans and Coastal Management Division 
monitoring strategy (Contact Dr. Dave Hebert - hito://www.cto-mpo.gc.ca/oceans/mpa- 


zpm/stanns-sainteanne-eng.htmi). 


11. Conduct hydrographic, chemical and biological sampling across the mouth of the 
Laurentian Channel and St. Pierre Bank. This transects have been implemented to 
enhance our understanding of hydrographic phenomenon in these areas in support of 
current modelling efforts (Contact Dr. Dave Brickman). 

12. Possibly deploy an ALSEAMAR glider in support of the DFO glider group (Contact Dr. 


Clark Richards). 


Table 4. Status of objectives upon completion of the HUD2018004 mission. 


Objective tatu Comments 


Partially Stations HL 12 to 14 were cancelled because of delays caused by CTD issues at HL 07 
Completed on April 14% and 15% 


Completed 


6 Completed Samples acquired at all stations identified prior to sailing. 
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SUMMARY OF ACTIVITIES 


CTD Summary 


Narrative 


As summarized in Table 1, there were a total of 95 CTD casts during the mission (Figure 4 and 
Table 6). Ten of these were aborted casts, which included: events 073 (YL 01), 098-104 
(HL 7),116 (SG 28), and 148 (CSL 02). 


Starting at event 001 in the Basin (HL, 00), the Valeport Altimeter (759017) was replaced with a 
Benthos altimeter (#49559) because the Valeport was not providing an accurate bottom distance 
(off by -50 m). A second test CTD cast in the Basin (event 002) threw a deck unit error, the cast 
was aborted and it was determined that the new cable for the recently swapped Benthos altimeter 
was the problem. From event 002 to event 035 the configuration file HUD2018004C was used 
(Table 5). Prior to event 037 at PL 09, the SeaPoint CDOM UV Fluorometer was removed. The 
data produced by this sensor has been jumping from what seems to be a reasonable value, to 0 
and back again during a single cast. The problem was observed as early as the spring 2017 
mission but the sensor had not been sent back for repair and there is no spare (a spare UV 
fluorometer should be purchased in 2018 if possible). The configuration file 
HUD2018004D.xmlcon was used between events 037 and 113 (HL 11) (Table 5). 


During event 073 (YL 01) the primary T/S/O sensors were returning unreasonable data values 
and the cast was aborted. Upon recovery, the primary plumbing was flushed with Triton and the 
pressure relief valve was replaced before the CTD was re-deployed at the same station for event 
074. The primary data continued to be incorrect upon descent, but upon ascent the blockage was 
removed and data values for the primary and secondary systems were in better agreement. For 
event 074 at YL 01 the secondary sensor values were entered into the AZMP template rather 
than the incorrect primary values. 


During event 098 at HL. 07 at ~2784 m on the upcast, the deck unit threw an error and an RS-232 
Communications Timeout. Prior to the redeployment at HL 07 during event 099, all sensors 
with the exception of the primary and secondary T/S/O and the Y (JT2) for the altimeter and the 
PAR were removed. The deck unit threw the same error at -2700 m of the upcast and the CTD 
was aborted. The secondary sensors and the Y to the altimeter and PAR were dummied prior to 
the next cast (event 100) so only the primary T/S/O sensors were connected to the 9+. This cast 
was successful, so the problem was narrowed down to the either one of the secondary sensors 
(T/S/O) or one of the sensors on the Y. When the secondary sensors were added on during event 
101, the CTD was also successful so the problem had something to do with the Y. Prior to event 
102, the Y was reattached to the 9+ with a new altimeter but the deck unit threw an error around 
550 m on descent, so the altimeter was not the problem. Now the issue had been narrowed down 
to the cables in the Y or the PAR sensor. Before the next cast (event 103) the Y cable had been 
replaced, but this cast also failed at 608 m on descent and was aborted. Finally, before event 
104, the altimeter cable leading to the Y was replaced. This failed as well upon descent, so prior 
to the next cast (event 105) the PAR was dummied off and the cast was successful. The trouble 
shooting at HL 07 took ~14 hours, leading to the cancellation of stations HL 12 — 14. 
Ultimately, the PAR sensor had failed and we did not have a spare so it was left dummied 
off for the remainder of the mission. 
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During the deep casts on the XHL, a periodic signal oscillation was noticed in the primary 
oxygen sensor that increased in amplitude as the casts grew deeper. As a result, after event 113 
(HL 11) the primary oxygen sensor (#0133) was replaced with #3030 for the remainder of the 
mission. From events 116 (SG 28) to event 159 (STAB 05) the configuration file 
HUD2018004e.xmlcon was used (Table 5). During event 116 the deck unit threw and error that 
was likely caused by a faulty PAR unit, which was removed from the CTD for the rest of the 
mission. We will need to service this PAR sensor and purchase another in 2018 as a spare. 
During event 121 (GULD 04), the primary oxygen sensor oscillations were observed again, so 
some caution should be used in interpreting deep water primary oxygen throughout the mission. 
At this point, the Y cable for the oxygen sensor and the primary oxygen had been replaced. After 
event 159 (STAB 05) on the long steam south to LL 08, the SBE 9 unit (21214) was replaced 
with unit 20370. The problem had been resolved during events 161 and 163 (LL 08 and LL 09) 
and configuration file HUD2018004f.xmlcon was used from event 161 to the end of the mission 
(event 180 at BANQ B1) (Table 5). 


Throughout the mission, 8 bottles on the CTD did not fire. As a result, as standard practice for 
both shifts after each cast, the CTD firing mechanism is now washed with Triton solution and 


rinsed with fresh water. 


Table 5. Configuration files used throughout the mission. 


Configuration File Event Range Problem Appendix 
HUD2018004b 001 Altimeter change la 
HUD2018004c 002 - 035 SeaPoint UV removed 1b 
HUD2018004d 037 - 113 Primary oxygen replaced 1c 
HUD2018004e 116 - 159 SBE 9 change Id 
HUD2018004f 161 — 180 le 
Conditions 


Preliminary section plots and anomalies (where available) of temperature (°C), salinity (P.S.U.) 
and sigma-t (kg/m^) in order of occupation (Browns Bank, Peter Smith, Portsmouth, Yarmouth, 
Halifax, Louisbourg, St. Anns Bank, Cabot Strait, and the mouth of the Laurentian Channel), can 
be viewed in Appendix 3. 


Water temperatures observed along the shelf break were much warmer than normal (Figure 2). 
In fact, water temperatures were in excess of 8 degrees warmer than normal at both BBL 07 and 
HL 07 (Appendix 3) and surface temperatures (—14 degrees) and salinities (—35.5 P.S.U.) were 
well mixed from the surface to approximately 200 m (Figure 3). While this is not entirely 
unprecedented, it is certainly rare within the 20 year time frame of the AZMP program. The 
conditions were unusual enough to incite some interest from news media who picked up the story 
during the cruise. This is the first time that temperature, salinity and density section plots had 
been created and disseminated during a mission and it is something that we will continue to do in 
the future. A special thanks to Roger Pettipas who created all of the plots and Dave Hebert who 
did the interviews with various news outlets: 
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http//nationalpost.com/pmn/news-pmn/canada-news-pmnu/researchers-find-abnormal 
atlantic-waters-off-nova-scotia 
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Discussions were undertaken post-cruise to determine the significance of these findings and 
whether further in-depth analyses are warranted. 


AVHRR Sea Surface Temperature Anomaly 
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Figure 2. Sea surface temperatures along the shelf break were much warmer than the 1981 — 
2010 climatology. 
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Figure 3. The temperature and salinity during the 
1997 to 2018. 
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Figure 4. Locations for the 95 CTD casts during HUD2018004 AZMP spring survey. Each cast 
is labelled with the consecutive mission event. 
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Oxygen 


The oxygen data collected by the CTD sensors and Winkler titration method will be used 
to create new calibration coefficients before the final run of the CTD processing. It will 
be necessary to extract these corrected oxygen values when they are produced so they can 
be accurately reflected in our data archives. 


The adjusted Soc values are calculated by a 2 step process. First, a “threshold field” 1s 
produced that subtracts the mean difference between the sensor and the average Winkler 
value for all samples, from the individual sample difference between the sensor and 
Winkler: 


(SBE O2 - Winkler O2) - mean(SBE O2 — Winkler O2) 
The next step calculates a new slope term by using the following equation: 
NewSoc = mean(previousSoc*([Winkler O2]/[SBE O2])) 


From events 001 to 113 (HL 11), the oxygen sensors 70133 (primary — calibrated on Dec 
1/17) and #0042 (secondary — calibrated on Dec 19/17) were deployed on the CTD. In 
deep water casts (below 1000 m), the primary oxygen sensor was displaying periodic 
oscillations that were not evident in the secondary. The decision was made to remove the 
primary after event 113 and replace it with #3030, but the oscillations continued. Finally, 
after event 159 (STAB 05), the SBE 9 (71214) was replaced with #0370 and the problem 
observed in Figure 5 was resolved. Note that the mean oxygen concentration during 
these oscillations seemed in line with conditions above and below. While it is 
recommended that the secondary oxygen should be utilized to interpret water column 
conditions, the primary data should not be thrown out. 
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Figure 5. The periodic deep water voltage oscillations observed for the primary oxygen 
sensors from events 001 to 159. 


The first step was to compare the Winkler replicates throughout the mission and remove 
any outlier data rows from further analysis (Figure 6.). In total, 12 of the total 80 (15%) 
rows where Winkler replicates were taken during the mission were removed prior to 
proceeding 


The analysis below is broken in to 2 parts for each of the primary sensors (#0133 and 
#3030). The oxygen report dataset was split into 2 subsets, <= event 113 and >= event 
113. 


Events 005 to 113 


Before the Soc can be calculated, some basic comparisons between the primary (#0133) 
and secondary (#0042) sensors for the first subset were completed to remove outliers and 
bad data (Figure 7). The 1.5 * inter quartile range was used to determine “outlier” data 
that could bias the results. From event 5 (BBL 01) to event 31 (BBL 05), there were 
many outliers between the primary and secondary sensor. This was driven by a 
malfunctioning secondary oxygen sensor that self-corrected after event 31 (Figure 8). 


The “threshold fields” were calculated for the both the primary (#0133) and secondary 
(#0042) sensors (Figures 9 and 10). Twelve more data rows were removed during the 
primary threshold calculation (Figure 9) before proceeding to secondary threshold 
calculations (Figure 10). The primary and secondary threshold calculations remove 13 
additional rows of data. 
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Table 7 shows the previous and revised Soc values for the primary and secondary oxygen 
sensors (#0133 and #0042) between events 005 and 113. The ratio of the new and old 


Soc values was calculated for each sensor. 


The original outlier free primary and secondary sensor values were then multiplied by 
their new corresponding Soc ratios to produce corrected primary and secondary sensor 
values (Figures 11 and 12) that more closely matched corresponding Winkler values. 
The difference between the primary and secondary sensors after correction (Figure 13) 
was slightly improved by the revised Soc values. 
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Figure 6. In total, 12 of 80 (15%) Winkler replicates were considered outliers and 
removed from further analysis. The mean difference between replicates was nearly zero 
at 4.4e-16. 
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Figure 7. The difference between primary oxygen sensor #0133 and secondary oxygen 
sensor #0042 for events 005 - 113. The mean difference before outlier removal (solid 
blue line) is ~ -0.04 ml/l. ~96% of the outlier difference between the primary and 
secondary were observed between events 5 and 31. 
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Figure 8. The primary oxygen sensor (A — blue line +0133) has a matching up and down 
trace. The secondary sensor (B — blue line #0042) overestimated on the down cast and 
came back to values similar to the primary on the up cast. This occurred between events 
005 and 030 but then self-corrected for the remainder of the mission. 
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Figure 9. Seven outlier “threshold” field values were removed for the primary (#0133) 
sensor after other outlier data (1.e., Winkler and sensor differences) were removed. 
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Figure 10. An additional 6 outlier “threshold” field values for the secondary sensor were 
removed and all of the remaining data were used to calculate the Soc value. 


Table 7. Previous and New Soc values for the primary and secondary SBE Oxygen 
sensors from events 005 to 113. 


| Old Soe Ratio (New:Old) 


Primary Sensor #0133 4.0488e-1 4.213026e-1 | 1.040562 
Secondary Sensor #0042 | 4.3400e-1 4.475889e-1 | 1.030598 
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Figure 11. Black dots — uncorrected outlier free sprimary sensor values (#0133) and 
Blue squares — Soc corrected primary sensor values. 
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Figure 12. Black dots — uncorrected outlier free secondary sensor values (#0042) and 
Blue squares — Soc corrected secondary sensor values. 
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Figure 13. Black dots — uncorrected difference between primary (#0133) and secondary 
sensors (#0042) between events 005 and 113 (mean = -3.8866e-002 ml/l). Blue squares — 
corrected difference between the primary and secondary sensors (mean = 2.0948e-002). 


Events 117 to 180 


Comparisons between the primary (#3030) and secondary (#0042) sensors for the second 
subset were completed to remove outliers and bad data (Figure 14). From event 117 
(SG 28) to event 180 (BANQ BI). there were 17 outliers between the primary and 
secondary sensor (Figure 14). Of these, 12/17 were from event 146 at CSL 01, where the 
secondary oxygen sensor malfunctioned throughout the cast. These data points were 
removed before proceeding to the threshold step. 


The “threshold fields" were calculated for the both the primary (#3030) and secondary 
(#0042) sensors (Figures 15 and 16). Six data rows were removed during the primary 
threshold calculation (Figure 15) before proceeding to secondary threshold calculations 
(Figure 16), where no further threshold outliers were observed. 


Table 8 shows the previous and revised Soc values for the primary and secondary oxygen 
sensors (#3030 and #0042) between events 117 and 180. The ratio of the new and old 
Soc values was calculated for each sensor. 


The original outlier free primary and secondary sensor values were then multiplied by 
their new corresponding Soc ratios to produce corrected primary and secondary sensor 
values (Figures 17 and 18) that more closely matched corresponding Winkler values. 
The calculated Soc for the secondary sensor (#0042) for each event subset is very similar 
(Tables 7 and 8), which shows that this sensor was fairly constant throughout the mission. 
Figure 19 shows the difference between the primary and secondary sensors after 
correction for this range of events. 
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Figure 14. Twelve of the 17 outliers between the primary (#3030) and secondary 
(#0042) oxygen sensors occurred during event 146 at CSL OI (red circles). 
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Figure 15. Six outlier “threshold” field values were removed for the primary (#3030) 
sensor after other outlier data (1.e., Winkler and sensor differences) were removed. 
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Figure 16. No additional “threshold” field values for the secondary sensor were 
removed. 


Table 8. Previous and New Soc values for the primary and secondary SBE Oxygen 
sensors from events 117 to 180. 


|| | Old Sec Ratio (New:Old) 
Primary Sensor #3030  |4.5600e-I | 4.872490e-1 | 1.068529 


Secondary Sensor #0042 | 4.3400e-1 4.488096e-1 | 1.033409 
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Figure 17. Black dots — uncorrected outlier free primary sensor values (#3030) and Blue 
squares — Soc corrected primary sensor values. 
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Figure 18. Black dots — uncorrected outlier free secondary sensor values (#0042) and 
Blue squares — Soc corrected secondary sensor values. 
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Figure 19. Black dots — uncorrected difference between primary (#3030) and secondary 
sensors (#0042) between events 117 and 180 (mean = -2.3013e-001 ml/l). Blue squares — 
corrected difference between the primary and secondary sensors (mean = 1.2820e-002). 
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Salinity 
(With portions extracted from HUD2014017 Cruise Report) 


Conductivity Calibration 


The salinometer outputs the conductivity as a ratio with the standard; therefore, some 
conversions are done to get the conductivity of the bottle. The standard has a given K15 
value: 


K15 = conductivity of standard seawater at 15°C and 1 atm/conductivity of KCI solution 
(32.4356g/kg) at 15°C and 1 atm. 


Where K15 = 0.99984 for this particular standard and the conductivity of KCI standard = 
4.29140 S/m and can be found in the seawater Matlab package (gsw C3515 function). 
Knowing K15 and the conductivity of the KCI solution, the conductivity of the standard 
seawater can be determined. Then, by multiplying by the conductivity ratio from the 
salinometer, the conductivity of the sample can be determined. 


It should be noted that these samples were analyzed with a bath temperature of 24?C 
rather than the 15?C that the standard conductivity was defined. The salinometer program 
accounted for this temperature difference so that the output sample conductivity ratios 
with the standard are at 15?C. 


Now we have the conductivity of the sample at 15?C and at the pressure of the bath in the 
salinometer; however, this needs to be converted to conductivity at the temperature and 


pressure of the CTD. This can be done using some functions from the same Matlab 
package. 


First calculate the salinity of the bottle using the conductivity and pressure from the 
salinometer and a temperature of 15?C. 


Salinity bottle = gsw SP from C(Conductivity salinometer[mS/ecm]|, T[C],P bath) 


Then re-calculate the conductivity from this salinity value using temperature and pressure 
from the CTD. 


Conductivity bottle = gsw C from SP(Salinity bottle, T CTD,P CTD) %[mS/cm] 


This now gives conductivity values that can be compared to the CTD values. To correct 
the CTD conductivity a linear regression is done on this equation: 


Bottle conductivity = b1 + b2*CTD conductivity 


to find an intercept, bl, and slope, b2, that will make the CTD conductivity better match 
the bottle conductivity. 


First, a comparison of the primary (#3562, Calibrated December 7, 2017) and secondary 
(21076, Calibrated December 6, 2017) sensor data (P.S.U.) was performed to highlight 
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and remove any outliers beyond 1.5 * the inter-quartile range of the data (Figure 20). 
This revealed 55 outliers that were removed from the analysis. Next, the difference 
between the primary sensor and salinometer values was compared in a similar manner to 
identify outliers that should be removed from analysis (Figure 21) (n=29). The same 
process was completed for the secondary sensor, and a single outlier was identified and 
removed before proceeding (Figure 22). After outliers were removed, it was determined 
that the primary and secondary salinometers were on average -0.0012 and 0.0015 P.S.U 
different from their corresponding salinometer values throughout the mission (Figure 23). 


At this point the swCSTp function, which uses the Gibbs-Sea Water (gsw C from SP) 
formulation, from the R OCE package, is used to convert the salinity of the bottle sample 
to conductivity. The data were filtered and used to create a linear regression for both the 
primary and secondary CTD sensor conductivity cells. The intercept (b1) and slope (b2) 
values for both regressions were extracted from the linear regression summary and used 
to correct the sensor values. These terms were used to calibrate the sensor salinity values 
for CTD output files prior to data archiving (Table 9). Figure 24, shows the difference 
between the primary and secondary sensors both before and after calibration with the new 
coefficients applied. 


Outliers Outside 1.5"1Q9R 


Primary Salinity - Secondary Salinity (FP... 
40 


40.15 


Ü 200 400 500 800 1000 1200 
Ordered by Event and Increasing Sample ID 
Figure 20. The outlier sensor values (red dots) that were removed prior to further 


analysis (n=55). The average difference between the primary (#3562) and the secondary 
(#1076) was ~-0.0022 P.S.U. throughout the mission. 
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Figure 21. The outlier differences between the primary sensor and the salinometer 
values (red dots) are removed prior to further analysis (n=29). 
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Figure 22. The outlier differences between the secondary sensor and the salinometer 
values (red dots) were also removed (n=1). 
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With Outiier Salinometer Data Removed 
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Figure 23. Note that after the outliers have been removed, the differences between the 


primary (#3562) and secondary (#1076) sensors and corresponding salinometer values 
are -0.0012 (black line) and 0.0015 (blue line) respectively. 
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Figure 24. The mean difference between the primary (#3562) and secondary (#1076) 
conductivity values before correction was -0.00184 mS/cm (black) and after correction is 
0.00042 mS/cm (blue). 
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Table 9. The revised intercept (b1) and slope (b2) terms calculated for both the primary 
(#3562) and secondary (#1076) conductivity sensors. 


Conductivity Sensor bl b2 
Primary (#3562) 7.5762e-03 0.999802 
Secondary (#1076) 1.5085e-02 0.999504 


Chlorophyll a 


Throughout the mission, Chl a was measured in-situ via a SeaPoint fluorometer attached 
to the CTD rosette (Appendix la). Duplicate samples (n=668) were regularly taken for 
Chl a analysis with a Turner Fluormeter. A comparison of the replicates showed that 
while the mean difference between replicates was 3.1903e-003 ug/L, there were a total of 
94 out of 668 replicates are considered outliers (Figure 25). Outliers were selected via 
the 1.5 * interquartile range (1.5 IQR) method discussed in the previous oxygen and 
salinity sections of this report. These outliers were removed before making the 
comparison between the SeaPoint sensor values and the Turner sensor values. 


Similar outlier identification methodology was employed to remove data that showed 
larger than expected differences between the SeaPoint sensor and the Turner Fluorometer 
data (Figure 26). First, both the SeaPoint data and the Turner data were standardized by 
dividing both data sets by the SeaPoint data value. This made each SeaPoint data value 
for a bottle fire equal to 1, and the corresponding mean replicate Turner fluorometer 
value a percentage of the SeaPoint value. A value of 1.15 means that the Turner 
Fluorometer value was 15% greater than its corresponding SeaPoint value and a value of 
0.85 means that the Turner value was 15% less than the SeaPoint value. This was done, 
because calculating the straight difference between values was influenced greatly by their 
magnitude. The difference between 0.01 and 0.1 and the difference between 6.31 and 6.4 
are both 0.09, but the relative difference is ~90% and ~1.4 % respectively. Figure 26 
shows the outliers calculated in this way. Out of 574 comparisons between the CTD 
sensor and the mean of the Turner Fluorometer replicates, 50 outliers were identified and 
removed before proceeding. The blue line shows that on average, SeaPoint sensor 
concentration values are ~68.2 % less than corresponding Turner fluorometer values. 


Figure 27 shows the log/log relationship between the SeaPoint Fluorometer values and 
the Mean Turner ChlA values with the outliers from Figure 26 highlighted in red. The 
black line corresponds to the 1:1 line. When the outliers are removed and a linear 
regression is applied to the log10/log10 linear relationship between the CTD sensor and 
the mean replicates (Figure 28), the relationship is strong and significant (R-squared: 
0.8905, p«2.2e-16). 
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Figure 25. The outlier Turner replicates removed prior to determining the relationship 
between the Turner Fluorometer and SeaPoint sensor values (n=94). 
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Figure 26. The outliers identified from calculating the % difference between Turner 
Fluorometer values and the SeaPoint sensor values (n=50). 
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Figure 27. The log10 scale plot of SeaPoint Fluorometer values and the corresponding 


mean replicate Turner Fluorometer values. Note the highlighted 1.5 * IQR outliers from 
Figure 26 in red. 
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Figure 28. The log10/log10 plot of SeaPoint Fluorometer values and the corresponding 
mean replicate Turner Fluorometer values colour coded by depth, where red and dark red 
are shallow (closer to the surface) and purple and blue are deep (closer to 100 m). 
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Water Samples for Chemical Analyses 


Station specific rosette bottle firing depths and water collections for chemical analysis 
can be found by referring to the CTD deck sheet binder and/or water chemistry sampling 
document prepared upon the conclusion of the mission and provided to ODIS. Table 6 
highlights CTD casts where water collections were made. 


Photosynthetically Active Radiation Sensor (PAR) 


The CTD was outfitted with a Satlantic Cosine PAR (irradiance) sensor (#1043, 
calibrated December 1, 2015) enclosed in a 7000 m titanium housing. Unfortunately, it 
was determined that the sensor failed during event 98 at HL 07 and no spare was 
available, so there 1s no in water PAR data available after this. 


pH Sensor 


The pH sensor (#1129, calibrated February 12, 2018) was deployed on the rosette only 
when the maximum depth was less than or equal to ~1200 m. The CTD casts for which it 
was deployed are noted in Table 6. The sensor was included during the mission to 
support an ACCASP initiative investigating the delineation of ocean acidification and 
calcium carbonate saturation state of the Atlantic zone. The sensor will soon become part 
of the core AZMP suite of sensors. 
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Biological Program 


Narrative 


The “core” biological program conducted as part of cruise HUD2018004, with some 
modifications, was a continuation of studies began in pre-AZMP years to describe the 
large-scale (spatial and temporal) variability in plankton biomass, productivity and 
biogenic carbon inventories on the Scotian Shelf. 


The program currently consists of essentially 3 elements: 


1. mesozooplankton community structure, population growth and biomass, and 
dissolved organic carbon measurements; 

3. Pigment analysis and flow cytometry for assessment of phyto-, micro- and nano- 
plankton communities. 


Table 6 provides a review of the stations where water samples were taken from rosette 
bottles for element 2 above. The mesoplankton sampling program is described below in 
more detail. This is followed by descriptions of “non-core” or ancillary biological 
sampling that includes a description of water sampling efforts in support of projects 
investigating how organic and organometallic micronutrients influence primary 
productivity and phytoplankton community structure on the Scotian Shelf (Erin Bertrand 
— Dalhousie University), and an assessment microbial communities and their associated 
processes (Julie LaRoche — Dalhousie University). The Biological Program section is 
concluded with a summary of pelagic seabird and marine mammal observations provided 
by Carina Gjerdrum of the Canadian Wildlife Service. 


No integrated phytoplankton sampling took place during this mission. 
The ultimate aim of “core” studies is twofold: 


1. to provide a description of the inventories of biogenic carbon, their turnover rates and 
variability in space and time as part of Ocean Ecosystem Science Division’s (OESD) 
continuing climate studies, and 

2. to provide a description of plankton life-cycles and productivity on the Scotian Shelf 
and its influence or contribution to ecosystems in support of OESD’s ecosystem- 
related research. 
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Mesozooplankton Sampling 


Remarks/Comments 


In order to estimate the mesozooplankton community abundance and biomass, a conical 
ring net of 202 um mesh size with an aperture of 75 cm in diameter (filtering ratio 1:5) 
equipped with a KC Denmark flow-meter was towed vertically from the bottom to the 
surface at most stations (or from a maximum depth of 1000m — AZMP standard). In 
total, there were 85 vertical ring net tows during the mission (Table 10, Figure 29). Of 
these, 1 was a 76 um mesh tow (30 cm diameter and 1:5 filtering ratio) at HL, 02 (Events 
78), and 29 were 202 um mesh tows at core stations along core AZMP sections (HL, 
BBL and LL) (Table 10 - objective 1). The 76 um net tows at HL 02 serve the same 
purpose of quantifying the community but targets a smaller fraction of the 
mesozooplankton community (i.e. smaller developmental stages, eggs and nauplii). 
Regardless of the mesh size, contents of the cod end were preserved in 4% buffered 
formaldehyde. 


Table 10 provides a summary of the vertical net tow operations conducted throughout the 
mission and the objectives that they support. 
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Figure 29. Locations for vertical ring net tows. Each tow is labelled with the 
consecutive mission event. 
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Microbial Protein and Organic Micronutrient Sampling 


Principle Investigator: Dr. Erin Bertrand (Dalhousie University, Department of 
Biology) 
Sampling by: Magda Waclawik and : _ (Dalhousie University) 


Objective 


The objective was to collect rosette samples for protein and vitamin analyses in order to 
determine whether and how organic and organometallic micronutrients influence primary 
productivity and phytoplankton community structure on the Scotian Shelf. Sampling 
locations were coordinated with the LaRoche lab since our datatypes are synergistically 
informative. 


Microbial Protein Sampling 
Purpose 


Proteins are key to microbial activity: the type and amount of proteins present 
determines, in large part, the contributions microbes make to the ecosystems they occupy. 
Proteins can also be used as indices for nutritional status: elevated expression of specific 
proteins can be diagnostic for different nutritional states, such as nitrogen starvation, iron 
starvation, or vitamin starvation. Protein sequences also contain taxonomic information 
and can be used to assess contributions of different organisms to specific functions. 


samples were collected for targeted, mass spectrometry- based proteomic analyses of 
microbial communities in order to characterize the role of organic micronutrients in 
structuring phytoplankton communities on the Scotian Shelf. Primary objectives include 
measuring phytoplankton nutritional status indicator proteins (nitrogen, vitamin B12, 
vitamin B, starvation) and vitamin- production biomarker proteins. Development and 
application of peptides for primary producer community composition analyses is a 
secondary focus. 


sampling Methods 


10L samples: A total of 71 size-fractionated microbial protein samples (10L of water 
each) were taken from the CTD rosette at depths ranging from the surface to 250 m 
(Table 11) along the lines of the AZMP located in Canadian waters. Water was filtered 
sequentially through 3 and 0.2 um polycarbonate filters via peristaltic pumping. Filters 
were then frozen immediately at -80?C. 


Vitamin Sampling 
Purpose 


To determine the particulate and dissolved concentrations of organic and organometallic 
micronutrients on the Scotian Shelf. Organic and organometallic micronutrients are 
required by many phytoplankton groups and only produced by a select few microbes, 
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setting up a series of interactive dependencies between microbial groups. The importance 
of these dependencies are not well known, as they have not yet been studied on the 
Scotian Shelf. Measuring the concentrations of these micronutrients in the particulate 
and dissolved phases is one step towards understanding the role of microbial interactions 
in driving primary productivity and phytoplankton community structure. 


sampling Methods 


A total of 71 particulate and 31 dissolved vitamin samples (1L each) were taken from the 
CTD rosette at depths ranging from the surface to 250m depth along lines of the AZMP 
in Canadian waters (Table 11). Samples were protected from light and gently vacuum 
filtered through 0.2 um nylon filters. Filters were frozen at -80?C and dissolved samples 
were frozen in amber HDPE bottles at -20?C. 


Table 11. Protein and vitamin samples — Bertrand lab — AZMP Spring 2018 — 
HUD2018004. 


Depth Sample Protein Protein Particulate Dissolved 
Station Event # 


(m) ID 0.2um Sum Vitamin Vitamin 
I 445581 I I I - 
BBLI 5 
40 445574 I I I - 
I 445609 I I I : 
20 445603 I I I - 
BBL3 9 
40 445600 I I I - 
80 445595 I I I - 
I 445760 I I I - 
20 445755 I I I - 
BBL5 3l 
40 445751 I I I - 
80 445746 I I I - 
I 445800 I I I - 
20 445796 I I I - 
BBL7 35 
80 445789 I I I - 
250 445784 I I I : 
I 446032 I I I I 
20 446028 I I I I 
HL02 79 
40 446024 I I I I 
80 446018 I I I I 
I 446073 I I I I 
20 446068 I I I I 
HL04 86 
40 446064 I I I I 
60 446060 I I I I 
HL06 92 I 446117 I I I I 
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Station 


HL07 


GULD04 


LL7 


LL4 


LE1 


STABO1 


CSLOI 


CSL04 


STABOS 


LL9 


Event + 


105 


121 


124 


130 


136 


138 


146 


153 


159 


163 


Depth 
(m) 


20 


Sample 
ID 


446113 
446108 
446104 
446178 
446174 
446170 
446323 
446319 
446315 
446311 
446347 
446342 
446335 
446330 
446383 
446379 
446375 
4463069 
446423 
446419 
446415 
446411 
446435 
446429 
446485 
446479 
446473 
446469 
446534 
446528 
446522 
446514 
446583 
446579 
446571 
446565 
446625 
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Depth Sample Protein Protein Particulate Dissolved 
Station Event # 


(m) ID 0.2um Sum Vitamin Vitamin 
20 446621 I I I I 
80 446613 I I I I 
250 446609 I I I I 
I 446689 I I I - 
20 446683 I I I - 
BANQ B6 174 
= 80 446677 - I I - 
250 446671 - I I - 
I 446709 I I I - 
20 446703 I I I - 
BANQ B4 176 
E 40 446699 I I I - 
80 446695 I I I - 
I 446727 I I I - 
BANQ BI 180 
20 446721 I I I - 
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Microbial Community Analysis 


Principle Investigator: Dr. Julie LaRoche (Dalhousie University) 
sampling by: Magda Waclawik and (Dalhousie University) 


Purpose 


Microbial communities and their associated processes are the foundation of marine life. 
Of particular interest to our group 1s the marine nitrogen cycle, comprising complex 
microbially-driven reactions whereby atmospheric nitrogen is fixed into a biologically- 
available form and cycled through the ecosystem. Though nitrogen is an essential 
element for life, the availability of fixed nitrogen can be a limiting factor for primary 
production and thus diazotrophs — organisms capable of biological nitrogen fixation — can 
be key to the productivity of an ecosystem. 


samples were collected for genomic and fluorescence-based analyses of the microbial 
communities on the Scotian shelf. Community composition will be assessed via 16S 
(bacterial) & 18S (eukaryotic) tag sequencing, and the naturally-fluorescent population 
will be characterized via flow cytometry. The latter method can also be used to quantify 
the bacterial community via nucleic acid stain SYBR green. Community function will be 
assessed via metagenomic sequencing, and qPCR assays for selected functional genes. 
Further samples were taken for manipulation in the lab, including targeted metagenomics 
and single cell 1solation via fluorescence-associated cell sorting (FACS), and enrichment 
culturing of putative diazotrophs. 


sampling Methods 


Genomics: 

At 19 select stations along the lines of the AZMP located in Canadian waters, duplicate 
4L water samples were collected from the CTD rosette at each of 4 depths ranging from 
the surface to 300m (Table 12); at 3 additional stations, only 2 depths were collected, for 
a total of 81 samples (Table 12). Each water sample was sequentially filtered through 3 
and 0.2um polycarbonate filters by peristaltic pump until the water was depleted or the 
filters clogged. Filters were immediately frozen at -80?C. 


Flow Cytometry: 

At each station and depth where genomic samples were collected, duplicate 2mL water 
samples were fixed with 2% paraformaldehyde (PFA) for 10 minutes at room 
temperature, then frozen at -80°C for later enumeration of bacteria and characterization of 
the naturally fluorescent microbial community via the Accuri C6 flow cytometer. 


At select stations (Table 12), 45mL of water were mixed with 5mL of gly-TE buffer and 
frozen at -80°C for later cell sorting on the BD Influx FACS instrument. 
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Enrichment Cultures: 

At select stations (Table 12), 500mL water samples were collected. These samples were 
spiked with phosphate (200nM) and iron (2nM) and secured to the window of the lab to 
approximate natural light/dark cycles and ambient temperature until return to the lab. 


Table 12. Microbial community samples — LaRoche lab — AZMP Spring 2018 — 


HUD2018004. 
DNA 
Depth samples boc pone 500 mL 
Station Event # P Sample ID P Flow Flow 
(m) (size- "em TE MAH culture 
fractionated) d zi P dd 
I 445580 2 2 - - 
BBLI 5 
40 445573 2 2 - - 
I 445608 2 2 - - 
20 445604 2 2 - - 
BBL3 9 
40 445599 2 2 - - 
80 445594 2 2 = I 
I 445759 2 2 I I 
20 445754 2 2 - - 
BBL5 3l 
40 445750 2 2 = N 
80 445745 2 2 I I 
I 445799 2 2 I I 
20 445795 2 2 - - 
BBL7 35 
80 445788 2 2 I I 
250 445783 2 2 - : 
I 446033 2 2 Š E 
20 446027 2 2 I I 
HL02 79 
40 446023 2 2 Š I 
80 446019 2 2 - - 
I 446072 2 2 - - 
20 446067 2 2 - - 
HL04 86 
40 446063 2 2 - - 
60 446059 2 2 - - 
I 446118 2 2 I I 
20 446114 2 2 - - 
HL06 92 
50 446109 2 2 I I 
80 446105 2 2 - - 
I 446179 2 2 - - 
20 446173 2 2 - - 
HL07 105 
50 446169 2 2 - - 
80 446167 2 2 - - 
I 4462 2 2 I I 
HLO8 107 ae 
20 446200 2 2 - - 
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Station 


HL11 


GULD04 


LE 


LLA 


LLI 


STABOI 


CSLOI 


CSL04 


STABOS 


LL9 


Event + 


113 


121 


124 


130 


136 


138 


146 


153 


159 


163 


Depth 


Sample ID 


446193 
446191 
446274 
446272 
446266 
446265 
446234 
446318 
446312 
446306 
446346 
446341 
446334 
446329 
446384 
446378 
446374 
446370 
446424 
446418 
446414 
446410 
446346 
446248 
446484 
446478 
446474 
446468 
446533 
446529 
446523 
446515 
446584 
446578 
446572 
446566 
446624 
446620 
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DNA 
Standard Sorting 


Depth ] L 
Station Event + n Sample ID P i Flow Flow S00 m 
s. e cytomet Cytomet a 
fractionated) d y y ry 
80 446614 2 2 Í _ 
250 446610 2 2 š - 
I 446688 2 2 š " 
20 446684 2 2 I Ë 
BANQ B6 174 
7 80 446676 l 2 : _ 
250 446672 2 2 E - 
I 446708 2 2 2 I 
20 446704 2 2 P _ 
BANQ B4 176 
u 40 446700 2 2 - E 
80 446694 2 2 _ _ 
I 446726 2 2 š 3 
BANQ BI 180 
i 20 446722 2 2 - _ 
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Pelagic Seabird and Marine Mammal Observations 


Seabird Survey Report 

7-24 April 2018 

Canadian Wildlife Service, Environment and Climate Change Canada 
Prepared by: Carina Gjerdrum | 
Observer(s): Jeannine Winkel 


Background 


The east coast of Canada supports millions of breeding marine birds as well as migrants 
from the southern hemisphere and northeastern Atlantic. In 2005, the Canadian Wildlife 
Service (CWS) of Environment Canada initiated the Eastern Canada Seabirds at Sea 
(ECSAS) program with the goal of identifying and minimizing the impacts of human 
activities on birds in the marine environment. Since that time, a scientifically rigorous 
protocol for collecting data at sea and a sophisticated geodatabase have been developed, 
relationships with industry and DFO to support offshore seabird observers have been 
established, and over 100,000 km of ocean track have been surveyed by CWS-trained 
observers. These data are now being used to identify and address threats to birds in their 
marine environment. In addition, data are collected on marine mammals, sea turtles, 
sharks, and other marine organisms when they are encountered. 


Methods 


Seabird surveys were conducted from the port side of the bridge of the CCG Hudson 
during the Scotian Shelf AZMP from 7-24 April 2018. Surveys were conducted while the 
ship was moving at speeds greater than 4 knots, looking forward and scanning a 90? arc 
to one side of the ship. All birds observed on the water within a 300m-wide transect were 
recorded, and we used the snapshot approach for flying birds (intermittent sampling 
based on the speed of the ship) to avoid overestimating abundance of birds flying in and 
out of transect. Distance sampling methods were incorporated to address the variation in 
bird detectability. Marine mammal and other marine wildlife observations were also 
recorded, although surveys were not specifically designed to detect marine mammals. 
Details of the methods used can be found in the CWS standardized protocol for pelagic 
seabird surveys from moving platforms’. 


'Gjerdrum, C., D.A. Fifield, and S.I. Wilhelm. 2012. Eastern Canada Seabirds at Sea (ECSAS) 
standardized protocol for pelagic seabird surveys from moving and stationary platforms. Canadian 
Wildlife Service Technical Report Series. No. 515. Atlantic Region. vi 36 pp. 


Results 
Seabird Sightings 


We surveyed 1854 km of ocean from 7-24 April 2018. A total of 1785 birds were 
observed in transect (2731 birds in total) from 11 families (Table 13). Bird densities 
averaged 3.1 birds/km” (ranging from 0 — 183.8 birds/km^). The highest densities of birds 
(> 50 birds/km^) were observed in primarily in the eastern part of the survey area, in 
Sydney Harbour, Misaine Bank, the mouth of the Laurentian Channel, northern 
Banquereau and slope waters beyond the Haldimand Canyon, but also in far western 
survey area in the Gulf of Maine (Figure 30). 
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Dovekie comprised of 15% of the birds sighted (Table 13). The Scotian Shelf (and 
Grand Banks of NL) are important wintering grounds for Dovekie breeding in Greenland. 
The bulk of the Dovekies observed on this survey were in slope waters off Banquereau 
(Figure 31a), likely heading back to their breeding grounds. Common Murres were also 
encountered relatively frequently, accounting for 4% of the sightings (Table 13), and 
were observed in the highest densities on Misaine Bank (Figure 31b). Together, the 
alcids made up a total of 26% of the birds observed. 


Herring Gulls accounted for 9% of the sightings (Table 13), and were observed in the 
highest densities on Banquereau (Figure 31c). Great Black-backed Gulls were less 
numerous, but similarly distributed. Both these species breed in Nova Scotia. Together 
the gull species accounted for 15% of the birds sighted. 


Northern Gannets were sighted throughout the survey area, but in the highest densities 
around SW Nova Scotia in more inshore waters (Figure 31d). Northern Gannets breed in 
NL and the Gulf of St. Lawrence and are common on the Scotian Shelf this time of year 
as they make their way back their breeding colonies from the Gulf of Mexico and areas 
along the eastern seaboard of the US. 


Marine Mammal Sightings 


A total of 134 marine mammals were recorded during the surveys (Table 14), 46% of 
which were humpback whales, primarily observed in the eastern sections of the survey 
(Figure 32a). A single pod of 20 white-beaked dolphins were observed on Misaine Bank 
(Figure 32b). 


Gully MPA 
surveys were conducted within the Gully MPA on 17 April. A total of 46 birds were 


observed, most of which were Dovekie (85%), but no marine mammals were observed in 
this area (Table 15; Figure 32c). 


St. Anns Bank MPA 
Surveys were conducted within the St. Anns Bank MPA on 18 April and then again on 21 


April. Just 19 birds and 3 unidentified cetacean species were observed during the transits 
(Table 16; Figure 32d). 
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Table 13. List of bird species observed during surveys on the Scotian Shelf AZMP, from 


7-24 April, 2018. 


Family 


Gaviidae 


Procellariidae 
Hydrobatidae 


Phalacrocoracidae 
Sulidae 


Anatidae 


Laridae 


Alcidae 


Picidae 


Turtidae 


Emberizidae 


Total 


English 


Common Loon 

Northern Fulmar 

Sooty Shearwater 

Manx Shearwater 
Leach's Storm-Petrel 
Wilson's Storm Petrel 
Unidentified Storm-Petrels 
Double-crested Cormorant 
Northern Gannet 
Canada Goose 

Surf Scoter 

Black Scoter 
White-winged Scoter 
Unidentified Scoters 
Common Eider 
Long-tailed Duck 
Unidentified Ducks 
Pomarine Jaeger 
Unidentified Jaegers 
Black-legged Kittiwake 
Herring Gull 

Great Black-backed Gull 
Lesser Black-backed Gull 
Glaucous Gull 

Iceland Gull 
Unidentified Gull 
Dovekie 

Common Murre 
Thick-billed Murre 
Unidentified Murre 
Black Guillemot 
Razorbill 

Atlantic Puffin 
Unidentified Alcid 
Northern Flicker 
American Crow 

Song Sparrow 

Savannah Sparrow 

Fox Sparrow 

Dark-eyed Junco 
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Latin 


Gavia immer 

Fulmarus glacialis 
Ardenna griseus 
Puffinus puffinus 
Oceanodroma leucorhoa 
Oceanites oceanicus 
Hydrobatidae 
Phalacrocorax auritus 
Morus bassanus 
Branta canadensis 
Melanitta perspicillata 
Melanitta nigra 
Melanitta fusca 
Melanitta 

Somateria mollissima 
Clangula hyemalis 

All duck genera 
Stercorarius pomarinus 
Stercorarius Jaegers 
Rissa tridactyla 

Larus argentatus 

Larus marinus 

Larus fuscus 

Larus hyperboreus 
Larus glaucoides 
Larus 

Alle alle 

Uria algaa 

Uria lomvia 

Uria 

Cepphus grylle 

Alca torda 

Fratercula arctica 
Alcidae 

Colaptes auratus 
Corvus brachyrhynchos 
Melospiza melodia 
Passerculus sandwichensis 
Passerella iliaca 

Junco hyemalis 


Number 
observed 
in 
transect 
5 
73 
12 


50 


25 
249 


Total 
number 
observed 


5 
97 
19 
50 


33 
422 
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Table 14. List of marine mammals, fish and invertebrates observed during surveys on the 


Scotian Shelf AZMP, from 7-24 April, 2018. 


English 


Humpback Whale 
Long-finned Pilot Whale 
Fin Whale 

Minke Whale 
White-beaked Dolphin 
Common Dolphin 
Harbour Porpoise 
Unidentified Dolphins 
Unidentified Whales and Dolphins 
Gray Seal 

Harbour Seal 

Harp Seal 

Unidentified Seals 
Ocean Sunfish 
Portuguese Man-Of-War 


Total 


Latin 


Megaptera novaeangliae 
Globicephala melas 
Balaenoptera physalus 
Balaenoptera acutorostrata 
Lagenorhynchus albirostris 
Delphinus delphis 
Phocoena phocoena 
Delphinidae 

Cetacea 

Halichoerus grypus 

Phoca vitulina 

Pagophilus groenlandicus 
Phocidae 

Mola mola 

Physalia physalia 
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Table 15. List of species observed in the Gully Marine Protected Area during surveys on 
the Scotian Shelf AZMP on 17 April, 2018. 


English Total number observed 
Dovekie 39 

Herring Gull 2 

Northern Gannet l 

Northern Fulmar I 

Pomarine Jaeger l 

Atlantic Puffin I 

Leach's Storm-petrel l 

Total 46 
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Table 16. List of species observed in the St. Anns Bank Marine Protected Area during 
surveys on the Scotian Shelf AZMP, on 18 and 21 April, 2018. 


English Total number observed 
Northern Gannet 8 
Northern Fulmar 6 
Herring Gull 3 
Black-legged Kittiwake I 
Common Murre | 
Unidentified cetacean 3 
Total 22 
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TOO SUPE 


Figure 30. Density of birds observed during the seabird survey on the Scotian Shelf 
AZ MP, from 7-24 April 2018. 
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Figure 31. Density of a) Dovekies, b) Murres, c) Herring and Black-backed Gulls, and 
d) Northern Gannets observed during the seabird survey on the Scotian Shelf AZMP, 
from 7-24 April 2018. 
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Underway Sampling 


Contributions by: Adam Hartling 
Division: OESD 


Navigation and Bathymetry 


The science navigation system on the Hudson consists of both Novatel ProPak6 and 
Applanix PosMV systems. The Novatel is a two antenna receiver that provides both 
position and heading data. The PosMV includes an IMU and provides both position and 
attitude data. The NMEA output from both systems is logged by the GSC Navnet 
software. The output from one of the receivers can be selected to be added to the 
multiplexed serial data stream available throughout the ship. The output from this 
selected receiver is also logged by the Scientific Computing System along with the TSG 
data. The output rate for both receivers 1s 1 Hz. 


The dual frequency Knudsen sounder was used for determining the water depth. The 
3.5kHz channel was used exclusively for depth sounding during the cruise. During CTD 
stations the Raytheon PTR/LSR sounder located in the winch-room was used for depth 
sounding and monitoring the distance from bottom of a 12kHz pinger mounted on the 
CTD package. The offset for the ram 1s 6m when retracted and 8m when extended. Both 
12kHz transducers are mounted on the ram. One 12kHz transducer on the ram was used 
for tracking the pinger, the other has failed and must be replaced in the next dry-docking. 


Underway Seawater System — Thermosalinograph 


An underway system, also referred to as thermosalinograph (TSG), was placed in the 
forward lab and was connected to the pumped uncontaminated seawater plumbing. The 
configuration on HUD2018004 consisted of an SBE21 with conductivity (sn: 3396) and 
temperature (sn: 3396), an external temperature located at the ship's intake (SBE 38 sn: 
0766), WET Labs chlorophyll WETStar (sn: WSCHL-1468) with a scale factor of 15.5 
ug/l/V, Seapoint CDOM fluorometer (sn: 6211) with a 50x gain jumper and SBE pH 
sensor (sn: 1214). The sampling rate was 0.2 Hz. 


The pump was turned on in Bedford Basin at -0800 LT on April 6”, just prior to 
departure. The water pumped to the forward lab with exhaust routes (direct discharge 
over the side of the ship), through the TSG and from the de-bubbler. The flow rate 
remained fairly constant throughout the mission at an average flow rate of ~44 l/min. 
The system remained running until ~0900 on the 19% of April when we arrived in Sydney 
for repairs. The system was re-started on the 20% at ~0900 in Sydney Harbour as we 
resumed operations. The system was finally turned off on April 24" at ~0700 as we 
began our final approach to Mulgrave, N.S. for crew change. 


Over the duration of the cruise a single PCO2 and TIC sample, along with 2 Chla samples 
were acquired every day and provided with a unique sample ID. The scanned paper log 
for these samples should be located here: 
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\\densbionad la BIODataSveSre 20108201 V HUD2018004 SCANNED LOGS and the 
digital e-logs can be found here: 
\\densbionaO 1a BlODataSveSre 20108201 SHUD2O0180043ELOGALogbookslow- 

log. There were a total of 15 PCO2, 15 TIC and 30 Chla samples taken from 
April 6 — April 24. 


TSG underway data was managed the NOAA Scientific Computing Systems (SCS) 
software. These data are submitted to ODIS upon conclusion of the mission but Dr. Dave 


Hebert (Dave.Hebert@dtfo-mpo.gc.ca) is the point of contact for these data. 


Data Management 


Prepared by: Robert Benjamin 
Division: Program Coordination and Support 


Data Collection 


In addition to standard AZMP manual data collection methods (1.e., Bridge log, various 
equipment specific deck sheets) ELOG, an electronic logbook system for collecting event 
metadata including position and sounding, was again used during HUD2018004. This 
electronic logbook was accessible via computers connected to the science network on- 
board the vessel with one available at each major work area. Metadata related to each 
piece of equipment was collected in the electronic log including position/time deployed, 
on bottom and recovered. Additional logbooks were employed to act as an itinerary, a 
daily operational log and a logbook to monitor the flow through system in the forward 
lab. All digital logbooks were backed up hourly and at the end of the Mission and were 
sent to ODIS for storage. 


Nav-Net, an on board ship's data collection system was used to collect all streaming data 
available during the entire mission with the exception of the Flow-through. These data 
include GPS data, sounder data, gyro data, wind and motion data. These data will be 
located in the NavNet archive here: Sdensbiona0 1a BIO DataSveSreNavnet 


Data Input Template 


Reports were generated from shipboard input data in the AZMP Template Database to 
compare with corresponding CTD sensor data and conduct preliminary analyses included 
in this report. 


GIS 


Daily navigation and operations were maintained in a graphical information system (QGIS). 
Final line and point shapefile were generated from these data for the cruise report. 


Hardware 
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Regulus/Aldebaran computers (supplied by NRCAN) were placed in the Drawing room, 
the CTD computer room, the Forward lab and the general purpose lab (GP Lab) to 
provide positioning and station name information to operations in these locations. 


Systems Development Project 


There was some discussion on board the ship for streamlining the ship-board methods for 
data collection and from subsequent shore-side analyses. The emphasis over the past 4 
years has been to improve the quality of ship-board and shore-side data and meta-data 
management. Various systems and controls have been employed as immediate stop-gaps 
(ELOG; AZMP database template; standardized QC; improved data submission 
guidelines; improved archival folder structures; developing and improving feedback 
mechanisms; audits, repairs or existing databases; creation of new databases; etc...). 
These efforts are well underway and in some cases complete or near completion, and they 
have dramatically improved the quality, accuracy and comprehensiveness of the data 
submitted by AZMP. These systems have dramatically improved our efficiency while 
also decreasing errors so researchers can use these data with confidence. Nonetheless, 
more effort should be placed on improving efficiencies and modernizing outdated 
systems and data management approaches. The AZMP Steering Committee and ODIS 
will work together to create a plan to modernize our data collection and management over 
the coming years. 


Already systems are being explored to help with sample tracking, water budget planning 
and improved data entry for at sea laboratory analysis. The development of a data entry 
system for Chlorophyll that has been analysed at sea is well underway. The Chlorophyll 
Input System includes well known and accepted formulas used to convert meter readings 
into valid Chlorophyll A and Phaeophytin values eliminating the need to use and modify 
Excel formulas when system calibrations or new Blanks need to be incorporated. The 
new Lab system will also inventory and manage all fluorometers and there calibrations. 
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APPENDICES 


Appendix 1. CTD Configuration Files 


Appendix la. HUD2018004.xmlcon (Event 1) 


Date: 04/04/2018 


Instrument configuration file: 


Z:\Currently Deployed\Shipboard Environment CRUISE SETUPS\HUD2018004\CTD ` 


Acquisition 2018004HUD'ctd con HUD2018004.xmlcon 


Configuration report for SBE 911plus/917plus CTD 


Frequency channels suppressed : 0 
Voltage words suppressed :0 


Computer interface : RS-232C 

Deck unit : SBE11plus Firmware Version >= 5.0 
scans to average a 

NMEA position data adde : Yes 

NMEA depth data added : No 

NMEA time added : No 


NMEA device connected to : deck unit 
surface PAR voltage added  : Yes 
Scan time added : No 


1) Frequency 0, Temperature 


Serial number : 5083 
Calibrated on : 02-Dec-17 


A : 3.68121240e-003 
B : 5.97265115e-004 
e : 1.50057603e-005 
D : 1.97182834e-006 
FO : 2984.722 

Slope : 1.00000000 


Offset : 0.0000 
2) Frequency 1, Conductivity 


Serial number : 3562 
Calibrated on : 07-Dec-17 


G : -9.842 1 1899e+000 
H : 1.20149777e+000 
I : -1.12798269e-003 
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J : 1.33106553e-004 


CTcor : 3.2500e-006 
CPcor : -0.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
3) Frequency 2, Pressure, Digiquartz with TC 


Serial number : 1214 
Calibrated on : 30-Nov-16 


Cl : -4,470905e+004 
(2 : 3.840789e-001 

C3 : 1.367850e-002 

DI : 3.661600e-002 

D2 : 0.000000e+000 

TI : 3.015271e+001 

T2 : -1.367200e-004 

T3 : 3.926620e-006 

T4 : 3.761680e-009 

T5 : 0.000000e+000 
Slope : 0.99999865 
Offset : -1.26180 
AD590M : 1.280000e-002 
AD590B : -9,348400e+000 


4) Frequency 3, Temperature, 2 


Serial number : 1376 
Calibrated on : 08-Dec-17 


A : 3.68121183e-003 
B : 6.00661106e-004 
C : 1.51504431e-005 
D : 2.12620555e-006 
FO : 6469.435 

Slope : 1.00000000 


Offset : 0.0000 
5) Frequency 4, Conductivity, 2 


Serial number : 1076 
Calibrated on : 06-Dec-17 


G : -4,18969113e+000 
H : 5.65236010e-001 

I : 3.20574520e-004 

J : 1.37368013e-005 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 
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Offset : 0.00000 
6) A/D voltage 0, Altimeter 
Serial number : 59017 
Calibrated on : 01-Mar-2017 
Scale factor : 15.000 
Offset : 0.000 
7) A/D voltage 1, PAR/Logarithmic, Satlantic 


Serial number : 1043 
Calibrated on : 01-Dec-2015 


a0 : 1.03324700 
al : 0.80736900 
Im : 1.35890000 
Multiplier : 1.00000000 


8) A/D voltage 2, Oxygen, SBE 43 


Serial number : 0133 
Calibrated on : 01-Dec-2017 
Equation : Sea-Bird 


Soc : 4.04880e-001 
Offset : -6.32800e-001 
A : -3.54420e-003 
B : 1.53460e-004 
C : -2.54240e-006 
E : 3.60000e-002 
Tau20 : 1.01000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
H1 : -3.30000e-002 
H2 : 5.00000e+003 
H3 : 1.45000e+003 


9) A/D voltage 3, Oxygen, SBE 43, 2 


Serial number : 0042 
Calibrated on : 19-Dec-2017 
Equation : Sea-Bird 


Soc : 4.34300e-001 
Offset : -4.93300e-001 
A : -3.30140e-003 
B : 1.49820e-004 

C : -1.91610e-006 
E : 3.60000e-002 
Tau20 : 1.47000e+000 


66 


000151 


D1 : 1.92634e-004 


D2 : -4.64803e-002 
HI : -3.30000e-002 
H2 : 5.00000e+003 
H3 : 1.45000e+003 


10) A/D voltage 4, Fluorometer, Seapoint Ultraviolet 


Serial number : 3668 
Calibrated on : 1-Jan-2015 
Range : 50.000000 
Offset : 0.000000 


11) A/D voltage 5, Fluorometer, Seapoint 


Serial number : 6210 
Calibrated on : 1-Jan-2015 
Gain setting : 3 x, 0-50 ug/l 
Offset : 0.000 


12) A/D voltage 6, pH 

Serial number : 1129 

Calibrated on : 12-Feb-2018 

pH slope : 4.6510 

pH offset : 2.5163 

13) A/D voltage 7, OBS, WET Labs, ECO-BB 
serial number : 1490 

Calibrated on : 9-Aug-2016 

ScaleFactor : 0.002983 

Dark output : 0.048000 

14) SPAR voltage, Unavailable 

15) SPAR voltage, SPAR/Surface Irradiance 
Serial number  : 1069 

Calibrated on : 24-Jun-2016 
Conversion factor : 1.00000000 

Ratio multiplier : 1.00000000 


Scan length : 40 
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Appendix 1b. HUD2018004C.xmlcon (Events 2-35) 
Date: 04/30/2018 


Instrument configuration file: 

CAUsersVCogswellA Documents AZMP Missions 201832018 Spring at 

sea\2018 HUD2018004\CTD\CTD ACQUISITION'2018004HUD'ctd con HUD20180 
04C.xmlcon * 


Configuration report for SBE 911plus/917plus CTD 


Frequency channels suppressed : 0 
Voltage words suppressed :0 


Computer interface : RS-232C 

Deck unit : SBE11plus Firmware Version >= 5.0 
scans to average :1 

NMEA position data adde : Yes 

NMEA depth data added : No 

NMEA time added : No 


NMEA device connected to : deck unit 
Surface PAR voltage added  : Yes 
Scan time added : NO 


1) Frequency 0, Temperature 


Serial number : 5083 
Calibrated on : 02-Dec-17 


A : 3.68121240e-003 
B : 5.97265115e-004 
C : 1.50057603e-005 
D : 1.97182834e-006 
FO : 2984.722 

Slope : 1.00000000 


Offset : 0.0000 
2) Frequency 1, Conductivity 


Serial number : 3562 
Calibrated on : 07-Dec-17 


G : -9,84211899e+000 
H : 1.20149777e+000 

I : -1.12798269e-003 

J : 1.33106553e-004 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
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3) Frequency 2, Pressure, Digiquartz with TC 


Serial number : 1214 
Calibrated on : 30-Nov-16 


Cl : -4.470905e+004 
C2 : 3.840789e-001 

C3 : 1.367850e-002 

DI : 3.661600e-002 

D2 : 0.000000e+000 

TI : 3.015271e+001 

12 : -1.367200e-004 

T3 : 3.926620e-006 

T4 : 3.761680e-009 

T5 : 0.000000e+000 
Slope : 0.99999865 
Offset : -1.26180 
AD590M : 1.280000e-002 
AD590B : -9,348400e+000 


4) Frequency 3, Temperature, 2 


Serial number : 1376 
Calibrated on : 08-Dec-17 


A : 3.68121183e-003 
B : 6.00661106e-004 
C : 1.51504431e-005 
D : 2.12620555e-006 
FO : 6469.435 

Slope : 1.00000000 


Offset : 0.0000 
5) Frequency 4, Conductivity, 2 


Serial number : 1076 
Calibrated on : 06-Dec-17 


G : -4,18969113e+000 
H : 5.65236010e-001 

I : 3.20574520e-004 

J : 1.37368013e-005 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
6) A/D voltage 0, Altimeter 


Serial number : 49559 
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Calibrated on : 18-Feb-2010 
Scale factor : 15.000 
Offset : 0.000 


7) A/D voltage 1, PAR/Logarithmic, Satlantic 


Serial number : 1043 
Calibrated on : 01-Dec-2015 


a0 : 1.03324700 
al : 0.80736900 
Im : 1.35890000 


8) A/D voltage 2, Oxygen, SBE 43 


Serial number : 0133 
Calibrated on : 01-Dec-2017 
Equation : Sea-Bird 


Soc : 4.04880e-001 
Offset : -6.32800e-001 
A : -3.54420e-003 
B : 1.53460e-004 
C : -2.54240e-006 
E : 3.60000e-002 
Tau20 : 1.01000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
HI : -3.30000e-002 
H2 : 5.00000e--003 
H3 : 1.45000e+003 


9) A/D voltage 3, Oxygen, SBE 43, 2 


Serial number : 0042 
Calibrated on : 19-Dec-2017 
Equation : Sea-Bird 


Soc : 4.34300e-001 
Offset : -4.93300e-001 
A : -3.30140e-003 
B : 1.49820e-004 
C : -1.91610e-006 
E : 3.60000e-002 
Tau20 : 1.47000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
H1 : -3.30000e-002 
H2 : 5.00000e+003 
H3 : 1.45000e+003 
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10) A/D voltage 4, Fluorometer, Seapoint Ultraviolet 


Serial number : 3668 
Calibrated on : 1-Jan-2015 
Range : 50.000000 
Offset : 0.000000 


11) A/D voltage 5, Fluorometer, Seapoint 


Serial number : 6210 
Calibrated on : 1-Jan-2015 
Gain setting : 3 x, 0-50 ug/l 
Offset : 0.000 


12) A/D voltage 6, pH 


Serial number : 1129 
Calibrated on : 12-Feb-2018 
pH slope : 4.6510 
pH offset : 2.5163 


13) A/D voltage 7, OBS, WET Labs, ECO-BB 
serial number : 1490 
Calibrated on : 9-Aug-2016 
ScaleFactor : 0.002983 
Dark output : 0.048000 

14) SPAR voltage, Unavailable 

15) SPAR voltage, SPAR/Surface Irradiance 
Serial number : 1069 
Calibrated on : 24-Jun-2016 
Conversion factor : 1.00000000 
Ratio multiplier : 1.00000000 


* - The configuration was changed after the file was opened. 
Scan length : 40 
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Appendix 1c. HUD2018004D.xmlcon (Events 37-113) 
Date: 04/10/2018 


Instrument configuration file: 
CACTD ACQUISITION '2018004HUD'ctd con HUD2018004D.xmlcon * 


Configuration report for SBE 911plus/917plus CTD 


Frequency channels suppressed : 0 
Voltage words suppressed :0 


Computer interface : RS-232C 

Deck unit : SBE11plus Firmware Version >= 5.0 
scans to average :1 

NMEA position data adde : Yes 

NMEA depth data added : No 

NMEA time added : No 


NMEA device connected to : deck unit 
Surface PAR voltage added  : Yes 
Scan time added : No 


1) Frequency 0, Temperature 


Serial number : 5083 
Calibrated on : 02-Dec-17 


A : 3.68121240e-003 
B : 5.97265115e-004 
C : 1.50057603e-005 
D : 1.97182834e-006 
FO : 2984.722 

Slope : 1.00000000 


Offset : 0.0000 
2) Frequency 1, Conductivity 


Serial number : 3562 
Calibrated on : 07-Dec-17 


G : -9,84211899e+000 
H : 1.20149777e+000 

I : -1.12798269e-003 

J : 1.33106553e-004 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
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3) Frequency 2, Pressure, Digiquartz with TC 


Serial number : 1214 
Calibrated on : 30-Nov-16 


Cl : -4.470905e+004 
C2 : 3.840789e-001 

C3 : 1.367850e-002 

D1 : 3.661600e-002 

D2 : 0.000000e+000 

T1 :3.015271e+001 

T2 : -1.367200e-004 

T3 : 3.926620e-006 

T4 : 3.761680e-009 

T5 : 0.000000e+000 
Slope : 0.99999865 
Offset : -1.26180 
AD590M : 1.280000e-002 
AD590B : -9,348400e+000 


4) Frequency 3, Temperature, 2 


Serial number : 1376 
Calibrated on : 08-Dec-17 


A : 3.68121183e-003 
B : 6.00661106e-004 
C : 1.51504431e-005 
D : 2.12620555e-006 
FO : 6469.435 

Slope : 1.00000000 


Offset : 0.0000 
5) Frequency 4, Conductivity, 2 


Serial number : 1076 
Calibrated on : 06-Dec-17 


G : -4,18969113e+000 
H : 5.65236010e-001 

I : 3.20574520e-004 

J : 1.37368013e-005 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
6) A/D voltage 0, Altimeter 


Serial number : 49559 
Calibrated on : 18-Feb-2010 
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Scale factor : 15.000 
Offset : 0.000 


7) A/D voltage 1, PAR/Logarithmic, Satlantic 


Serial number : 1043 
Calibrated on : 01-Dec-2015 


a0 : 1.03324700 
al : 0.80736900 
Im : 1.35890000 
Multiplier : 1.00000000 


8) A/D voltage 2, Oxygen, SBE 43 


Serial number : 0133 
Calibrated on : 01-Dec-2017 
Equation : Sea-Bird 


Soc : 4.04880e-001 
Offset : -6.32800e-001 
A : -3.54420e-003 
B : 1.53460e-004 
C : -2.54240e-006 
E : 3.60000e-002 
Tau20 : 1.01000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
HI : -3.30000e-002 
H2 : 5.00000e--003 
H3 : 1.45000e+003 


9) A/D voltage 3, Oxygen, SBE 43, 2 


Serial number : 0042 
Calibrated on : 19-Dec-2017 
Equation  : Sea-Bird 


Soc : 4.34300e-001 
Offset : -4.93300e-001 
A : -3.30140e-003 
B : 1.49820e-004 
C : -1.91610e-006 
E : 3.60000e-002 
Tau20 : 1.47000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
H1 : -3.30000e-002 
H2 : 5.00000e+003 
H3 : 1.45000e+003 
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10) A/D voltage 4, Free 
11) A/D voltage 5, Fluorometer, Seapoint 


Serial number : 6210 
Calibrated on : 1-Jan-2015 
Gain setting : 3 x, 0-50 ug/l 
Offset : 0.000 


12) A/D voltage 6, pH 


Serial number : 1129 
Calibrated on : 12-Feb-2018 
pH slope : 4.6510 
pH offset : 2.5163 


13) A/D voltage 7, OBS, WET Labs, ECO-BB 
Serial number : 1490 
Calibrated on : 9-Aug-2016 
ScaleFactor : 0.002983 
Dark output : 0.048000 

14) SPAR voltage, Unavailable 

15) SPAR voltage, SPAR/Surface Irradiance 
Serial number  : 1069 
Calibrated on : 24-Jun-2016 
Conversion factor : 1.00000000 
Ratio multiplier : 1.00000000 


* - The configuration was changed after the file was opened. 
Scan length : 40 
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Appendix 1d. HUD2018004E.xmlcon (Events 116-159) 
Date: 04/16/2018 


Instrument configuration file: 
CACTD ACQUISITION\2018004HUD\ctd con HUD2018004E.xmlcon * 


Configuration report for SBE 911plus/917plus CTD 


Frequency channels suppressed : 0 
Voltage words suppressed :0 


Computer interface : RS-232C 

Deck unit : SBE11plus Firmware Version >= 5.0 
scans to average :1 

NMEA position data adde : Yes 

NMEA depth data added : No 

NMEA time added : No 


NMEA device connected to : deck unit 
Surface PAR voltage added  : Yes 
Scan time added : No 


1) Frequency 0, Temperature 


Serial number : 5083 
Calibrated on : 02-Dec-17 


A : 3.68121240e-003 
B : 5.97265115e-004 
C : 1.50057603e-005 
D : 1.97182834e-006 
FO : 2984.722 

Slope : 1.00000000 


Offset : 0.0000 
2) Frequency 1, Conductivity 


Serial number : 3562 
Calibrated on : 07-Dec-17 


G : -9,84211899e+000 
H : 1.20149777e+000 

I : -1.12798269e-003 

J : 1.33106553e-004 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
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3) Frequency 2, Pressure, Digiquartz with TC 


Serial number : 1214 
Calibrated on : 30-Nov-16 


Cl : -4.470905e+004 
C2 : 3.840789e-001 

C3 : 1.367850e-002 

D1 : 3.661600e-002 

D2 : 0.000000e+000 

T1 :3.015271e+001 

T2 : -1.367200e-004 

T3 : 3.926620e-006 

T4 : 3.761680e-009 

T5 : 0.000000e+000 
Slope : 0.99999865 
Offset : -1.26180 
AD590M : 1.280000e-002 
AD590B : -9,348400e+000 


4) Frequency 3, Temperature, 2 


Serial number : 1376 
Calibrated on : 08-Dec-17 


A : 3.68121183e-003 
B : 6.00661106e-004 
C : 1.51504431e-005 
D : 2.12620555e-006 
FO : 6469.435 

Slope : 1.00000000 


Offset : 0.0000 
5) Frequency 4, Conductivity, 2 


Serial number : 1076 
Calibrated on : 06-Dec-17 


G : -4,18969113e+000 
H : 5.65236010e-001 

I : 3.20574520e-004 

J : 1.37368013e-005 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
6) A/D voltage 0, Altimeter 


Serial number : 49559 
Calibrated on : 18-Feb-2010 
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Scale factor : 15.000 
Offset : 0.000 


7) A/D voltage 1, PAR/Logarithmic, Satlantic 


Serial number : 1043 
Calibrated on : 01-Dec-2015 


a0 : 1.03324700 
al : 0.80736900 
Im : 1.35890000 
Multiplier : 1.00000000 


8) A/D voltage 2, Oxygen, SBE 43 


Serial number : 3030 
Calibrated on : 21-Dec-2016 
Equation : Sea-Bird 


Soc : 4.56000e-001 
Offset : -5.24800e-001 
A : -3.41420e-003 
B : 1.83300e-004 
C : -2.88040e-006 
E : 3.60000e-002 
Tau20 : 1.11000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
HI : -3.30000e-002 
H2 : 5.00000e--003 
H3 : 1.45000e+003 


9) A/D voltage 3, Oxygen, SBE 43, 2 


Serial number : 0042 
Calibrated on : 19-Dec-2017 
Equation  : Sea-Bird 


Soc : 4.34300e-001 
Offset : -4.93300e-001 
A : -3.30140e-003 
B : 1.49820e-004 
C : -1.91610e-006 
E : 3.60000e-002 
Tau20 : 1.47000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
H1 : -3.30000e-002 
H2 : 5.00000e+003 
H3 : 1.45000e+003 
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10) A/D voltage 4, Free 
11) A/D voltage 5, Fluorometer, Seapoint 


Serial number : 6210 
Calibrated on : 1-Jan-2015 
Gain setting : 3 x, 0-50 ug/l 
Offset : 0.000 


12) A/D voltage 6, pH 


Serial number : 1129 
Calibrated on : 12-Feb-2018 
pH slope : 4.6510 
pH offset : 2.5163 


13) A/D voltage 7, OBS, WET Labs, ECO-BB 
Serial number : 1490 
Calibrated on : 9-Aug-2016 
ScaleFactor : 0.002983 
Dark output : 0.048000 

14) SPAR voltage, Unavailable 

15) SPAR voltage, SPAR/Surface Irradiance 
Serial number  : 1069 
Calibrated on : 24-Jun-2016 
Conversion factor : 1.00000000 
Ratio multiplier : 1.00000000 


* - The configuration was changed after the file was opened. 
Scan length : 40 
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Appendix le. HUD2018004F.xmlcon (Events 161—180) 
Date: 04/21/2018 


Instrument configuration file: 
CACTD ACQUISITION'2018004HUD'ctd con HUD2018004F.xmlcon * 


Configuration report for SBE 911plus/917plus CTD 


Frequency channels suppressed : 0 
Voltage words suppressed :0 


Computer interface : RS-232C 

Deck unit : SBE11plus Firmware Version >= 5.0 
scans to average : 1 

NMEA position data adde : Yes 

NMEA depth data added : No 

NMEA time added : No 


NMEA device connected to : deck unit 
Surface PAR voltage added  : Yes 
Scan time added : No 


1) Frequency 0, Temperature 


Serial number : 5083 
Calibrated on : 02-Dec-17 


A : 3.68121240e-003 
B : 5.97265115e-004 
C : 1.50057603e-005 
D : 1.97182834e-006 
FO : 2984.722 

Slope : 1.00000000 


Offset : 0.0000 
2) Frequency 1, Conductivity 


Serial number : 3562 
Calibrated on : 07-Dec-17 


G : -9,84211899e+000 
H : 1.20149777e+000 

I : -1.12798269e-003 

J : 1.33106553e-004 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
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3) Frequency 2, Pressure, Digiquartz with TC 


Serial number : 0370 
Calibrated on : 19-Jan-17 


Cl : -4,274542e+004 
C2 : 1.040996e-002 

C3 : 1.266000e-002 

D1 : 4.087300e-002 

D2 : 0.000000e+000 

T1 : 3.009606e+001 

T2 : -6.521164e-005 

T3 : 4.345040e-006 

T4 : 2.428800e-009 

T5 : 0.000000e+000 
Slope : 0.99940215 
Offset : -0.24955 
AD590M : 1.289670e-002 
AD590B : -8.390780e+000 


4) Frequency 3, Temperature, 2 


Serial number : 1376 
Calibrated on : 08-Dec-17 


A : 3.68121183e-003 
B : 6.00661106e-004 
C : 1.51504431e-005 
D : 2.12620555e-006 
FO : 6469.435 

Slope : 1.00000000 


Offset : 0.0000 
5) Frequency 4, Conductivity, 2 


Serial number : 1076 
Calibrated on : 06-Dec-17 


G : -4,18969113e+000 
H : 5.65236010e-001 

I : 3.20574520e-004 

J : 1.37368013e-005 
CTcor : 3.2500e-006 
CPcor : -9.57000000e-008 
Slope : 1.00000000 


Offset : 0.00000 
6) A/D voltage 0, Altimeter 


Serial number : 49559 
Calibrated on : 18-Feb-2010 
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Scale factor : 15.000 
Offset : 0.000 


7) A/D voltage 1, PAR/Logarithmic, Satlantic 


Serial number : 1043 
Calibrated on : 01-Dec-2015 


a0 : 1.03324700 
al : 0.80736900 
Im : 1.35890000 
Multiplier : 1.00000000 


8) A/D voltage 2, Oxygen, SBE 43 


Serial number : 3030 
Calibrated on : 21-Dec-2016 
Equation : Sea-Bird 


Soc : 4.56000e-001 
Offset : -5.24800e-001 
A : -3.41420e-003 
B : 1.83300e-004 
C : -2.88040e-006 
E : 3.60000e-002 
Tau20 : 1.11000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
HI : -3.30000e-002 
H2 : 5.00000e--003 
H3 : 1.45000e+003 


9) A/D voltage 3, Oxygen, SBE 43, 2 


Serial number : 0042 
Calibrated on : 19-Dec-2017 
Equation  : Sea-Bird 


Soc : 4.34300e-001 
Offset : -4.93300e-001 
A : -3.30140e-003 
B : 1.49820e-004 
C : -1.91610e-006 
E : 3.60000e-002 
Tau20 : 1.47000e+000 
D1 : 1.92634e-004 
D2 : -4.64803e-002 
H1 : -3.30000e-002 
H2 : 5.00000e+003 
H3 : 1.45000e+003 
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10) A/D voltage 4, Free 
11) A/D voltage 5, Fluorometer, Seapoint 


Serial number : 6210 
Calibrated on : 1-Jan-2015 
Gain setting : 3 x, 0-50 ug/l 
Offset : 0.000 


12) A/D voltage 6, pH 


Serial number : 1129 
Calibrated on : 12-Feb-2018 
pH slope : 4.6510 
pH offset : 2.5163 


13) A/D voltage 7, OBS, WET Labs, ECO-BB 
Serial number : 1490 
Calibrated on : 9-Aug-2016 
ScaleFactor : 0.002983 
Dark output : 0.048000 

14) SPAR voltage, Unavailable 

15) SPAR voltage, SPAR/Surface Irradiance 
Serial number  : 1069 
Calibrated on : 24-Jun-2016 
Conversion factor : 1.00000000 
Ratio multiplier : 1.00000000 


* - The configuration was changed after the file was opened. 
Scan length : 40 
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Appendix 3. Preliminary Section Plots and Anomalies (T/S/Sigma-T) 
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St. Anns Bank Line 
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Cabot Strait Line 
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ABSTRACT 


This analysis is intended to begin the dialogue required to develop a framework for monitoring 
and assessing spatially managed areas such as Marine Protected Areas (MPAs). In particular, 
we begin with the Oceans Act requirement to describe productivity, biodiversity, habitat and 
species of interest by identifying pragmatic metrics of these ecosystem-level attributes from pre- 
existing monitoring systems in the area of interest: St. Anns Bank in the Maritimes Region of 
Canada. We also identify a few human influences/pressures that are also readily quantified, 
namely fishing and vessel activity and some of the data gaps that are evident in the area. The 
nature of these data, and the steps required to use them properly, are made explicit in an open- 
sourced and revision-controlled environment (R and git) for the purpose of developing a 
transparent, vetted data and code system from which future assessment and modeling attempts 
can be staged. This should reduce the replication of effort in future. The rationale for the 
methods chosen to quantify the characteristic space and time scales of processes and features 
are identified and discussed. Ultimately, the intent is to develop approaches similar to the risk- 
based approach frequently encountered in fishery stock assessments, such that we can begin 
to express the “status” of an area or MPA. It is perhaps even possible to attempt to define 
reference points and obtain a better sense of the "health" of these areas, or at least assess the 
relative influence of area-based management closures. This first report will focus only upon the 
“Data” side of the issues and clarify the proposed “Methods”. 
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INTRODUCTION 
TERMS OF REFERENCE 


the conservation and restoration of lands and waters in Canada. In this context, the Science 
Branch of the Department of Fisheries and Oceans (DFO), has been tasked with developing a 
monitoring approach for Marine Protected Areas (MPAs) and, if possible, to assess their 
effectiveness in meeting their objectives. 


SCOPE OF THIS REPORT 


An Marine Protect Area (MPA) is defined in the Oceans Act 35(1) as, "An area of the sea that 
forms part of the internal waters of Canada, the territorial sea of Canada or the Exclusive 
Economic Zone of Canada and has been designated under this section for special protection for 
one or more of the following reasons: 


e conservation and protection of commercial and non-commercial fishery resources, including 
marine mammals, and their habitats; 


e conservation and protection of endangered or threatened marine species, and their habitats; 
e conservation and protection of unique habitats; 


e conservation and protection of marine areas of high biodiversity or biological productivity; 
and 


e conservation and protection of any other marine resource or habitat as is necessary to fulfill 
the mandate of the Minister." 


As a result, for the purposes of this report, we will likewise focus upon these key ecosystem 
attributes or components, namely: productivity, biodiversity, habitat, and species of 
interest. In reality, however, there are many other attributes or components of ecosystems that 
are also known to be important and relevant, including: ecological integrity and health, trophic 
structure and balance, ecosystem function, complexity, network structure, resilience, 
sustainability, as well as an open-ended number of species or life-history stages of the resident 
species. These other components will be touched upon where possible or necessary. 


ST. ANNS BANK 


The St. Anns Bank Marine Protected Area (herein, SAB) is an area of interest for eventual 
designation as an MPA. It is located east of Cape Breton Island, Nova Scotia, Canada 

(Figure 1). Previous advisory processes (DFO 2012; Kenchington 2013), identified the primary 
objectives of SAB as being to conserve, protect and, where appropriate, restore ecologically 
distinctive or significant areas and, overall, the ecosystem "health" of SAB. As in the Oceans 
Act, the focus was upon the above four ecosystem components: productivity, biodiversity, 
habitat and species of interest. 


Other MPA goals were also expressed in DFO (2012) and Kenchington (2013), but these were 
less emphatic: 


e contribute to the health, resilience and restoration of the Eastern Scotian Shelf ecosystem; 
e contribute to the recovery and sustainability of commercial fisheries; and 


e promote scientific research and monitoring to further understand and protect SAB. 
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Ford and Serdynska (2013) make more precise the ecological components that the SAB area of 
interest might help to protect and conserve, especially in the context of the definition of MPAs in 
the Oceans Act: 


e commercial and non-commercial fishery resources, including marine mammals, and their 
habitats (e.g., habitat for Atlantic Cod, redfish, American Plaice, sea urchins, White Hake, 
Witch Flounder, sea anemones, sponges, and sea pens); 


e Endangered or Threatened marine species, and their habitats (e.g., habitat for depleted 
species such as Atlantic Wolffish, Atlantic Cod, and Leatherback Turtles); 


e unique habitats (it is the only major bank on the Inner Scotian Shelf); and 


e marine areas of high biodiversity or biological productivity of invertebrates and fish. 


OBJECTIVES 
The primary objectives of this report are to: 
e Develop an Assessment Framework that can 


o monitor the status of an MPA; and 
o assess the effectiveness of an MPA in meeting its conservation objectives. 


e Identify data gaps and sources of uncertainty. 


The method by which this can be accomplished, however, is anything but straightforward. This 
is primarily due to the fact that the SAB is: 


e alarge ecosystem and as such complex, operating at various space, time and 
organizational scales; 


e connected in various ways to the surrounding environment and so cannot be treated as an 
isolated system; and 


e measures of system components of interest, namely, productivity, biodiversity, habitat and 
species of interest, are ambiguous and imperfect at best, and usually non-existent or poor in 
information quality and/or quantity. 


As such, we emphasize that this report is a simplistic first attempt at developing a general 
approach towards assessing MPAs. Indeed, given the above significant challenges, it is best 
viewed as a work in progress that will require further precision and improvement. To this end, 
we will, in the following, describe the data used for the assessment; outline the methods and 
assumptions associated with the modeling of this data; summarize the primary results of this 
analysis; provide some discussion of salient points; and conclude with general 
recommendations. The technical aspects of data Quality Assurance (QA)/Quality Control (QC), 
associated assumptions and methods are identified in the Appendices. 


NOTE: The primary purpose of this document is the Data and Methods. Results and 
Discussions will be the focus of a subsequent paper. 


DATA 


In this section, we focus upon a description of the data chosen for inclusion in this assessment. 
The purpose of the section is to clearly identify the data, sampling design, and the associated 
assumptions and methods required/used to filter and integrate them in an informative manner. 
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For the sake of transparency, all data assimilation and QA/QC methods have been encoded in 
an open-sourced analytical environment R (R Core Team 2015) and made publicly accessible 
at GitHub ( ) and to permit flexible and adaptive multiuser 
contributions through the Git revision control system and a uniform data interface system. This 
approach permits the development of a coherent and vetted approach that is completely open 
source in nature. In this way, we see this project as a real and flexible structural and 
collaborative framework in the sense of a real scaffolding, to build a monitoring and assessment 
system that can be easily transferred to other regions, domains and mandates, as well as, 
fostering collaborations and communication with universities and the general public. 


STUDY AREA 


Evaluating MPA status and effectiveness in meeting conservation objectives, requires explicit 
reference to changes both within and without the area of interest. Even in the most basic BACI- 
type design, this requirement is explicit (Green 1979; Underwood 1992). For this reason, and 
also to facilitate evaluations of other potential MPAs in the region, a much larger surrounding 
area was chosen for analytical purposes. This area is the continental shelf region of Nova 
Scotia (Figure 2), bounded by latitudes 37°N to 48°N and longitudes from 48°W to 71°W. [Note: 
It should be emphasized that this will not alleviate problems associated with pseudo-(spatial, 
temporal) replication (Hurlbert 1984), although an assumption of a Gaussian process (see 
Methods) may potentially alleviate this problem.] 


DATA SELECTION CRITERIA 


Exhaustive surveys of available data have been compiled by Ford and Serdynska (2013). Their 
conclusions were that most biological data and environmental conditions are poorly sampled in 
the SAB area. The decision criteria for inclusion of data in this study were as follows: 


e Part of an on-going sampling program. This is because the design principle behind this 
project is that the underlying assessment must smoothly transition into a routine monitoring 
approach into the future. 


e Sufficient and regular spatial coverage (> 100 sampling locations) inside MPA and 
throughout the study area. 


e Sufficient and regular temporal coverage (approximately annually, > 10 years) inside MPA 
and throughout the study area. 


e Informative — high data quality that is in some manner related to productivity, biodiversity, 
habitat and species of interest. 


The same decision criteria were applied to human usage data. The result was to include the 
following data streams for MPA characterisation: 


e Atlantic Zone Monitoring Program (AZMP)/chlorophyll-a and nutrients: BioChem bottle data 
(Devine et al. 2014). 


e AZMP/Zooplankton: BioChem database (Devine et al. 2014). 
e Remote Sensing Data: ocean colour and SST (Remote Sensing Group). 


e Groundfish: DFO's Groundfish Research Vessel surveys focus upon demersal fish species, 
since approximately 2000, upon invertebrates as well. 


e Snow Crab survey, focus upon benthic invertebrates. 
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e Clam survey data in Banquereau and Western banks (though it does not pass the temporal 
coverage conditions, it offers very high resolution multispecies data on the banks). 


e Temperature records: from various sources, especially, groundfish, Snow Crab and AZMP 
surveys. 


e Salinity (Groundfish surveys/AZMP, BioChem). 

e Oxygen and pH (once the data have been reloaded; Groundfish surveys/AZMP, BioChem). 
e Bathymetry (CHS, Groundfish survey, Snow Crab survey). 

To characterise human usage patterns, the following have been chosen for inclusion: 

e Logbook records of catch and effort (MARFIS/ZIFF). 

e AIS tracks — Radio-based Automatic Identification System. 


e VMS potentially — Satellite-based Vessel Monitoring System. 


DISCRETE BOTTLE DATA: CHLOROPHYLL-A AND NUTRIENTS 

e Relevance: productivity, biodiversity, habitat and species of interest (in relative order). 
e Sampling: AZMP surveys, Groundfish surveys, pelagic net tows and water profiles. 

e Spatial coverage: variable no. stations, 143,499 records, 829 missions. 

e Temporal coverage: 1955 to present, annual surveys. 

e Source code: htt El 


e Discrete bottle data consisting of chlorophyll-a and nutrient records (nitrate, phosphate and 
silicate) were obtained by laboratory analysis of water samples collected at discrete depths. 
For this study all available nutrient and chlorophyll-a discrete bottle data were extracted from 
DFO's BioChem database for the study area. This dataset contains more than 500,000 
records with the earliest records starting in 1955. After QA/QC, the discrete bottle data 
retained for analysis was comprised of 143,499 profiles, collected on 829 missions 
(Figure 3; Appendix 1). 


The number of profiles available in each year (Figure 4) shows that there were few profiles 
taken until the mid-1960s, and a relatively steady number of yearly profiles after the initiation of 
the AZMP in 1999. The peak sampling during the period 1976-1982 corresponds to DFO's 
ocotian Shelf Ichthyoplankton Program (SSIP) and foreign research vessels sampling programs 
that were obtained from the National Oceanic and Atmospheric Administration's (NOAA) 
National Oceanographic Data Center (Pierre Clement, personal communication). Monthly 
distribution of the profiles (Figure 5) demonstrates that most of the data was collected in July 
(mostly during DFO's groundfish surveys), followed by the months of September and April. Note 
that spatial distribution of the sampling varies among months with most data collected on the 
ocotian Shelf in July and the fewest data in January (Figure 3 and Figure 6). 


ZOOPLANKTON DATA 

e Relevance: productivity, biodiversity, species of interest, habitat (in relative order). 
e Sampling: AZMP surveys, groundfish surveys, pelagic net tows, 400 taxa. 

e Spatial coverage: 2367 net deployments, 126 missions. 


e Temporal coverage: 1999 to 2014, annual surveys. 
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e Source code: https://github.com/j¡ae0/aegis/ 


Number of net deployments for each month is shown in Figure 7 and the corresponding spatial 
distribution of the net deployments are shown in Figure 8. Note that most of the net data were 
collected in July during summer groundfish survey missions and in April and October on AZMP 
spring and fall missions, while winter months contain mostly fixed station data (Halifax 2 and 
Prince 5). Abundance patterns are found in Figure 9. The QA/QC issues are identified in 
Appendix 1. 


REMOTE SENSING DATA 


Ocean Colour 

e Relevance: productivity, habitat, biodiversity and species of interest (in relative order). 
e Sampling: MODIS. 

e Spatial coverage: 39°N to 62.5°N and 42°W to 71"W, resolution of 1.5 km. 


e Temporal coverage: August 2002 to March 2015, 610 quarter-monthly (8-day) composite 
images. 


e Source code: https://github.com/j¡ae0/aegis/ 


Ocean colour refers to the spectral distribution of light emerging from the ocean that carries 
information about water constituents, particularly about biologically useful chlorophyll 
concentration in the surface layer. When measured from satellites, it provides unique synoptic 
view of chlorophyll spatial distribution over large areas of the ocean on a daily time scale. 


The nominal uncertainty of chlorophyll products derived from ocean colour satellites is 35%, 
with better agreement with in-situ chlorophyll for the open ocean (Moore et al. 2008), while 
overestimation is often observed in the coastal ocean (Darecki and Stramski 2004). This is due 
to the inability of the algorithms to distinguish chlorophyll from suspended particulate matter and 
colored dissolved organic matter often present in the coastal waters as, for example, in the Bay 
of Fundy and Northumberland Strait. 


Ocean colour satellite data for this study was provided by the Remote Sensing Unit (RSU) from 
the Bedford Institute of Oceanography (BIO) as 8-day composite chlorophyll images, which are 
routinely produced by the unit for the AZMP. The dataset was created using the Moderate 
Resolution Imaging Spectroradiometer (MODIS-Aqua) data, where the chlorophyll-a values are 
based on the 2012 reprocessing carried by the National Aeronautics and Space Administration 
(NASA) using OC3 chlorophyll algorithm. Composite images were created from daily Level-2 
MODIS-Aqua files downloaded from NASA by averaging valid chlorophyll-a values for each 
pixel using all available daily images within that time period (Caverhill et al. 2015; Feldman and 
McClain 2015). The dataset is comprised of years 2002 to 2015 and area 39°N to 62.5°N and 
42°W to 71°W, with spatial resolution of about 1.5 km per pixel. 


Even though there is ocean color data available before the MODIS launch, it was decided to 
limit our dataset to a single sensor to avoid potential biases between the sensors. Due to the 
frequent cloud coverage of the Northwest Atlantic, it was decided to use 8-day composite 
images as daily images would not provide a sufficient number of valid pixels. 


The values of chlorophyll-a pixels within St. Anns Bank polygon were extracted from each 8-day 
composite image and average chlorophyll-a concentration was computed for the polygon. The 
time series for the polygon and the associated climatology that was computed from time series 
data show characteristic spring blooms in March and April, with varying intensity and timing 
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throughout the years (Figures 10 and 11). An example of MODIS semi-monthly chlorophyll-a 
products showing spring bloom progression in the St. Anns Bank area in 2012 is shown on 
Figure 12. 


Primary Production 

e Relevance: productivity, habitat, biodiversity and species of interest (in relative order). 
e Sampling: MODIS. 

e Spatial coverage: 39°N to 62.5%N and 42°W to 71°W, resolution of 1.5 km. 

e Temporal coverage: July 2002 to December 2014, 150 monthly images. 


e Source code: https://github.com/j¡ae0/aegis/ 


Marine primary production plays an important role in biogeochemical cycles, in food web 
dynamics, and in marine fisheries. It may be defined as the amount of organic material (or 
organic carbon) produced per unit area (or volume) per unit time by photosynthetic plants, 
predominately by phytoplankton. 


Primary production of the ocean on synoptic scale is estimated by models that use satellite data 
(Sea Surface Temperature (SST), ocean colour chlorophyll, and available light estimates at the 
surface), shipborne in-situ information on the vertical distribution of phytoplankton in the water 
column, and the phytoplanktons photosynthetic response to light. 


Monthly primary production data were provided by the Remote Sensing Unit at BIO that 
routinely generates production maps for the Northwest Atlantic as part of the AZMP. The 
general approach for the production computation is described in Platt et al. (2008) and employs 
chlorophyll and light estimates from MODIS, SST produced by the unit, and DFO's archive of 
shipborne measured parameters. This particular production algorithm has been validated with 
in-situ measured production (Platt and Sathyendranath 1988) and also has performed very well 
in global comparisons (Carr et al. 2006). 


Primary production for each pixel within St. Anns Bank polygon was extracted from the monthly 
images, and average production was computed for the polygon. The time series for the polygon 
and the associated monthly climatology are showing peaks in primary production in spring and 
summer, with varying intensity and timing throughout the years (Figure 13 and Figure 14). 


Sea Surface Temperature 

e Relevance: productivity, habitat, biodiversity and species of interest (in relative order). 
e Sampling: AVHRR. 

e Spatial coverage: 39°N to 62.5% N and 42"W to /1°W, resolution of 1.5 km. 

e Temporal coverage: December 1997 to March 2015, 845 8-day composite images. 


e Source code: https://github.com/j¡ae0/aegis/ 


Sea Surface Temperature (SST) from space was estimated using infrared channels of the 
Advanced Very High Resolution Radiometer (AVHRR) on board the polar-orbiting satellites. 


The SST data for this study was provided by the Remote Sensing Unit from the BIO that has 
been downlinking AVHRR data from NOAA satellites since 1997 on an L-band satellite receiver 
that resides on the roof of the BIO. They routinely process the received data and supplement it 
with data stream from the AVHRR onboard the European satellites. Composite SST images of 
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different periods are operationally produced as part of the AZMP. Here we used 8-day 
composite images with the same spatial coverage and spatial resolution as the ocean colour 
chlorophyll data. 


The SST pixels within St. Anns Bank polygon were extracted from each 8-day composite image 
and average SST was computed for the polygon. The time series for the polygon and the 
associated climatology that was computed from time series data are shown on Figures 15 

and 16. An example of semi-monthly SST products corresponding to the spring bloom 
progression in the St. Anns Bank area in 2012 is shown in Figure 17. 

BOTTOM TEMPERATURES 

e Relevance: productivity, habitat, biodiversity and species of interest. 

e Sampling: Groundfish survey, Snow Crab survey, AZMP profiles. 

e Spatial coverage: full extent, varied sampling. 

e Temporal coverage: 1950 - present (more historical data present but coverage is variable). 
e Source code: https://qithub.com/jae0/aegis/ 


Numerous data sources have been compiled by Ocean Sciences Division, DFO. The data were 
QA/QC controlled and then modeled in a two-part process, temporal (first order harmonic 
analysis) and then spatial interpolation as indicated in the Methods (Figure 18). 


DEMERSAL FISH AND MACRO-INVERTEBRATES 

e Relevance: productivity, habitat, biodiversity and species of interest. 
e Sampling: Groundfish survey, Snow Crab survey. 

e Spatial coverage. 


o Groundfish: full extent, random stratified, variable number of stations. 
o Snow Crab: Colder water environment, geostatistical grids of approximately 10 minutes 
and approximately 400 stations. 


e Temporal coverage. 


o Groundfish: 2000 - present (started in 1970, but consistent sampling since 2000). 
o Snow Crab: 2005 - present (started in 1996, but consistent sampling since 2005). 


e Source code: 


o Groundiish: https //aithub.com/iae0/aeais/ 
o snow Crab: https //athub. com/AaeO/bio snowcrab/ 


Groundfish Survey 


The Groundfish survey (Figure 19, left) utilizes a Western II-A Otter Trawl net with a wingspread 
that is assumed to be 12.5 m and a target distance of 1.75 nautical miles (3.24 km) and/or a 
20-30 min tow. This net was used from 1982 to the present. Prior to this period, a Yankee 36 ft 
trawl was used with unmeasured net configuration data. It operates night and day. Sampling 
occurrence as a function of time and season are shown in Figure 19 (right). The consistent 
identification of invertebrates in this survey began in approximately the year 2000. All species 
assemblage analyses will use data from the post 2000 period. 
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Net Mensuration 


Sensors measuring trawl net configuration and state are ubiquitous in modern surveys and 
commercial fishing practices. In the Snow Crab survey, net configuration has been monitored 
and used to determine swept area manually (since 1996) and also with an automated procedure 
(since 2004). Unfortunately, in the groundfish survey, this information is ignored, even though 
the survey series is a major source of information for many stock assessments. Net 
configuration data have been collected in the groundfish surveys sporadically since 2004 
(Figure 20). However, swept area estimates have never been directly computed. Instead, the 
catch data is used almost invariably with the assumption that each tow is equivalent in swept 
area (a “standard tow’). Alternatively, it has also been used by some assessments by 
"adjusting" catch data by the linear distance from “start” and “end” times and locations and so 
implicitly assuming the net configuration to be a constant (Don Clark, personal communication). 
Both approaches are problematic for the reasons identified below. 


To address these important and incorrect assumptions, algorithms were developed to 
automatically determine lift off and touch down times, locations and net width (Munden, J., and 
J.S. Choi. (in prep) Calculatina Swept Area for the Maritimes Region Trawl Survey). Based upon 
random visual inspection, the skill was determined to be reasonable, with >90% of the cases 
having estimates within 15 seconds of visual determination of touch down/lift off locations. 
Where bias was observed, this was mostly to underestimate total contact time due to an over- 


smoothing of the lift off or touch down profiles. 


° 


From this re-analysis, the actual time and distance of bottom contact was found to range from 
approximately 10 to 40 minutes and distances from 1.75 to 6 km. Historical studies assumed a 
"standard tow" of 3.24 km in length. Even when compared to the "length" of a tow determined 
from positions recorded of “start” and “stop”, times or locations were frequently in error relative 
to actual locations of net touch down and lift-off determined from the Global Positioning System 
(GPS) information (Figure 21). 


Net structure also varies along tows and between tows as it encounters differing substrate, 
bathymetric and hydrodynamic conditions and vessel speeds, currents, surface sea state, and 
net fullness/filtration efficiency due to contact with rocks, boulders and mud. Assuming net width 
to be a constant at 12.5 m is, therefore, problematic (Figures 22 and 23). 


When the variable nature of the length and width of survey tracks are both accounted for, the 
potential error in swept area estimates are evidently large (Figure 24). Indeed, the range in 
variation of swept area was approximately the same in magnitude as the swept area of a 
standard tow. This means that an unnecessary and large "observation error”, potentially as 
large in magnitude as the magnitude of the catch in a set, may be needlessly introduced to 
survey indices. Further, swept area estimates based upon standard tows were also biased 
down relative to estimates based upon net mensuration and actual tow tracks, meaning that 
there is a high probability of over-estimating catch densities. 


Unfortunately, net mensuration data is still not consistently recorded nor used in the groundfish 
surveys (Figure 20). As a result, it is not possible to satisfactorily estimate swept area for all 
historical sets. Statistical methods of recalibration were used to impute swept area in these 
latter sets (Figure 24), using GAM-based estimates based on relationships with location, depth 
and salinity (R* ranged from 40 to 6096, depending upon availability of covariates). However, 
due to issues with a significant bias detected in wingspread sensors since 2013 (Figure 23), this 
approach cannot be extended for the period post-2013. 


To conclude this section, the observation error (uncertainty) associated with this data series, 
used in so many research programs, is high and biased high. Removal of these errors by a 
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coherent representation of net configuration is simple to address and a major data gap that 
needs to be addressed before this data series can be used quantitatively. 


Snow Crab Survey 


The Snow Crab survey uses a Bigouden Nephrops trawl (Figure 19), a net originally designed to 
dig into soft sediments for the capture of lobsters in Europe (headline of 20 m, 27.3 m foot rope 
mounted with a 3.2 m long 8 mm chain, with a mesh size of 80 mm in the wings and 60 mm in 
the belly and 40 mm in the cod-end). Tows were conducted for approximately 5 minutes in 
duration with duration of bottom contact being monitored by Netmind and Seabird sensors. The 
width of the mouth of the net was also monitored with Netmind sensors. The ship speed was 
maintained at approximately 2 knots. The warp length was approximately 3 times the depth. 
Positional information, as well as, water temperature measurements, was collected using a 
global positioning system and Minilog and Seabird data recorders. The surface area swept by 
the net was calculated from swept distance and net width information. 


Supplemental SAB Stations 


The 2015 Snow Crab trawl survey increased sampling in the St. Anns Bank area to provide 
additional information about the marine macro-fauna. Fourteen stations (14) in, and adjacent to, 
the proposed MPA location were included in this additional sampling. These locations were 
close to previous Snow Crab survey stations and represent varied depths, bottom-types and 
habitats. The species composition of the catch at these stations varied greatly, as expected with 
differences in depth and bottom type. The sampling at these stations included: 


e all catch identified to species level; 

e all species counted and weighed to a tenth of a kilogram; 

e all finfish and crab species measured and weighed individually; and 
e Stomach samples taken from finfish for diet studies. 


An overview (in Google Earth format) of the latest catch and sampling at these stations can be 


found at: htip//vww.enssnowcrab.com/mpa/mpatows.kmz. 


FISHERIES ACTIVITY 

e Relevance: productivity, habitat, biodiversity and species of interest. 
e Sampling: MARFIS and ZIFF. 

e Spatial coverage: full extent. 

e Temporal coverage: 1999 — present. 

e Source code: 


o Net Vensuration: https://github.com/¡ae0/netmensuration 
o Data Wrangling: https://github.com/Maritines/Mar.datawrangling/ 


Commercial-fishing activities can modify the habitat and ecosystem and contribute to changes 
in the structure and functioning of exploited marine communities. Fishing impacts can be direct, 
such as the reduction of targeted and non-targeted species, as well as truncations in age and 
size distributions. Other direct effects due to fishing activities include habitat alterations and 
substrate modifications through interactions with fishing gear. Fishing can also cause indirect 
impacts via changes in foodweb structure within an ecosystem. 
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Direct impacts of commercial fishing can be measured using data from the Marine Fish 
(MARFIS) database that provides information on commercial-fishing activities. For most 
fisheries, the fishing position, gear type, catch per unit effort, and estimated weight of catch by 
species information is available from the database. The MARFIS database details information 
for all fishing trips where a landing is reported within the DFO Maritimes Region and includes 
data from 2002 through 2015. The exploitation of marine species and entanglement threat to 
cetaceans and sea turtles will be quantified from data derived from the MARFIS database. 


Trawling and dredging disturbances to the sea floor will require different estimation techniques. 
Vessel Monitoring System (VMS) point locations have been used to estimate fishing-effort 
distributions (e.g., Lee et al. 2010) and to estimate impacts on the seabed (Gerritsen et al. 
2013). Similar techniques can be developed within R to estimate the impacts to benthos from 
trawling and dredging activities within the St Anns Bank area. 


The St. Anns Bank Area of Interest has four zones within it with various levels of fishing 
restrictions (Figure 1). Commercial fishing is restricted in Zone 1, the largest area, with the 
exception of the seal harvest. The MPA requires monitoring to ensure that fishers are complying 
with these regulations. Monitoring could be done through a variety of techniques and data 
sources, including data reported in the MARFIS database and VMS data that allow the direct 
monitoring of fishing activities. Automatic Identification System (AIS) data may also be used to 
monitoring fishing activities if the vessel is large enough to require AIS (2 300 gross registered 
tonnes) or if a fishing vessel has voluntarily installed an AIS system onboard. 


MARFIS Data Extraction 


Prorated landings for all species from 2002 onwards were extracted from MARFIS. The 
proration process distributes the actual reported weights across the reported efforts as identified 
within the fishers' logs. 


In addition to the landings, we also included several other forms of catch so as to better reflect 
the removal of biomass. These catches were identified by their CATCH USAGECODE and 
include biomass used as bait or discarded (sometimes identified as dead). Live discards were 
not included. These catches were self-reported in a variety of units, so they have been 
converted to kilograms as necessary to match the logged landings. 


Most log records included the spatial location of the catch, and some catches have multiple sets 
of coordinates available within the table MARFISSCI, LOGEFRTSTDINFOID, ENTLATITUDE 
and ENTLONGITUDE are physically entered into the logbook by the fisher, and they were used 
preferentially over DETLATITUDE and DETLONGITUDE, which are determined from other 
sources (like Loran-C). 


Many logs have no usable coordinates since they are either left blank or are clearly incorrect 
(i.e. on land). Rather than discarding this data, we still extracted it and attempted to account for 
this biomass in the next section. 


MARFIS Data Quality Control and Aggregation 


Quality Control (QC) occurs in 2 stages and is quite simple. First, all records without coordinates 
are identified and retained, but they are removed from the main dataset. Next, all remaining 
data are compared against a high resolution coastline (the same as is shown in Figure 2), and 
those points falling on land are identified and retained but removed from the main dataset. 


The remaining data are considered good and are aggregated. The aggregation level is user- 
defined. A scale of 3 minutes is used by default since it is an even division of a degree, with no 
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potential for rounding errors. The aggregation process outputs a single record with a single 
position for all of the catch in the area. 


Following aggregation, the proportion of the total catch attributable to each cell is calculated. 
The data that failed the QC tests are then summed into a single value, representing the total 
catch that cannot be positioned. These data are then added to all of the cells in the same 
proportion as was calculated in the previous step. 


For example, one cell might have a total catch of 3269./ kg, and this cell represents 0.002167 
of the total catch. If there are 5954.1 kgs of un-positioned data, then the corrected value for the 
weight attributed to this cell would be calculated as: 3269.7 kg + (55954.1 kg * 0.002167) = 
3390.9 kg. 


An example of the annually aggregated (2010) reported commercial catches for sea scallops 
(Placopecten magellanicus) is presented in Figure 25. Nominal catches ranged to 14,073.8 kg 
in Bay of Fundy, Georges Bank, and the Scotian Shelf. 


oimilarly, an example of annually aggregated (2011) reported commercial catches for Atlantic 
Halibut (Hippoglossus hippoglossus) is shown in Figure 26. Nominal catches ranged to 
approximately 9626.9 kg in the same area. 


VESSEL ACTIVITY 
e Relevance: habitat, biodiversity and species of interest. 
e Sampling: AIS. 


e Spatial coverage: Global for satellite AIS, coastal (approximately 100 km) for Canadian 
Coast Guard terrestrial AIS network. 


e Temporal coverage: 2013 — present. 


e Source code: https://github.com/j¡ae0/aegis/ 


Commercial shipping can have various direct and indirect effects on an ecosystem. Direct 
effects including the contamination of the ecosystem from the discharge of contaminates, 
radiated underwater noise, the introduction of aquatic invasive species, and vessel-strike risk to 
marine mammals and sea turtles. Spatial-temporal data on vessel traffic is needed to determine 
the probability and/or magnitude of these effects on ecosystems. Such information is available 
through AIS data. 


The International Maritime Organization (IMO) requires AIS transponders on all international 
vessels z 300 gross tonnage and all passenger vessels. Many studies have used AIS data to 
examine risk of lethal vessel collisions to large whales (e.g., Vanderlaan and Taggart 2009; 
Wiley et al. 2011; Redfern et al. 2013; Guzman et al. 2013) or to assess and monitor ship noise 
and assess the impact on marine mammals (Hatch et al. 2008; McKenna et al. 2012; Hatch 

et al. 2012; Merchant et al. 2014). Similar exercises can be undertaken in the St. Anns Bank 
area with AIS data. 


Fisheries and Oceans Canada (DFO) has access to two different sources of AIS data. The first 
is from the Canadian Coast Guard (CCG) terrestrial system that was developed to track and 
monitoring coastal shipping and provides a real-time, continuous stream of AIS vessel positions. 
Archived historical data from this system is available for January 2012 through December 2015, 
and data from 2016 are currently streaming and archiving to a server within DFO. Decoding 
routines have been developed using native R methods for these data. Both sources of data 
provide dynamic and static data, where the dynamic data includes information on vessel 
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identity, speed, and location, and static data includes information vessel identity, name, size, 
and type. 


The CCG terrestrial system have 21 AIS coastal receiving stations in the Maritime region and 19 
AlS coastal receiving stations in Newfoundland and Labrador. These receiving stations have 
limited range for detecting vessels (Figure 27) as AIS transmission detectability are primarily a 
function of the receiving tower height above sea level and the height of the AIS antenna on the 
transmitting vessel. AIS data is transmitted via Very High Frequency (VHF) marine radio on two 
channels (161.9765 and 162.025 MHz). Based on the height of the associated towers and a 
vessel with an AIS antenna 100 m high, line of sight calculations for VHF provide a reception 
range of 57 to 113 m in and around St. Anns Bank (Figure 28). However, there are several other 
factors that will contribute to transmission range, including weather conditions. Simard et al. 
(2014) estimated that coastal antennas within this network generally provide a reception range 
of 100 km (Figure 29). In either case, the terrestrial network is insufficient to monitor vessels 
across the entire AO! and just north of the AOI. These data can be combined with satellite AIS 
data. 


The satellite AIS data are available globally for the years 2013, 2014, and 2015 with on-going 
data collection for 2016. Although satellite AIS data coverage is global, data are limited 
temporally as large sections of vessel transits are unavailable due to a limited number of 
satellites (n = 8) and their orbital paths (see Figure 30). Spatial interpolation must be completed 
to fill in missing data. Spatial interpolation is achieved using an A* function (Hart et al. 1968) that 
estimates the minimum cost to get from one point to another based on a cost map. Using 
seasonally aggregated annual density distributions of satellite AIS data for the years 2013 
through 2015 (Figure 31), cost maps were estimated (Figure 32). Grid resolution for this 
analysis was initially set to 0.01 degrees and, within each grid square, the number of unique 
vessels identified by the vessel's Maritime Mobile Service Identity (MMSI) number was counted 
daily and summed through time. Cost maps were estimated quarterly. Two different cost maps 
were developed for the St. Anns Bank area to interpolate vessel transits north of Cape Breton 
into and out of the Gulf of St. Lawrence. Interpolation was heavily influenced by the ferries 
transiting between Cape Breton and Newfoundland and, therefore, a cost map was developed 
without the data derived from these ferries. A bathymetric restriction can also be built into the 
cost maps. 


DATA GAPS 


Feeding Relationships — Stomach Database 


Bottom trawl surveys have been conducted by DFO on the Scotian Shelf annuallt since 1970 
using a stratified random design. Sampling protocols changed in the late 1990s with the focus 
changing from commercially important finfish species to more comprehensive ecosystem 
monitoring that included the sampling of macro invertebrates (Tremblay et al. 2007). Stomach 
contents samples were collected from finfish using a length-stratified sampling protocol. Prey 
were quantified by weight and number, and they are often identified to the genus or family level, 
or to the species level when possible. 


For the purpose of determining any change in the diet of finfish, pre- and post-implementation of 
the MPA on St. Anns Bank, or if there were differences within the area compared to other areas 
on the Scotian Shelf, we explored the stomach database to determine if these dietary 
differences could be described and detected. Within the database, 54% of the prey number 
observations were missing. Due to the large interannual variation in prey weights, estimating the 
prey number consistently from non-missing data where both the prey number and weight was 
available was determined to be impracticable and unreliable. Due to the sampling stratification 


12 


000201 


both by depth and length classes, it was further determined that there would be insufficient 
samples available to reach the asymptote of prey species accumulation curves (Cook and 
Bundy 2010) and total diet composition could not be detected. Prey species identified in the 
stomach samples could be used for quantifying biodiversity and species richness (Cook and 
Bundy 2012) for the proposed MPA and comparing it to similar ecosystem or pre- and post- 
implementation. 


Other Ecosystem Metrics 


On the side of ecosystem characterization, data gaps in the following are evident: pelagic fish 
(small and large bodied) and invertebrates (e.g., squid, jellyfish), substrate characterisation via 
multibeam surveys, marine mammals, reptiles, birds and genetic diversity. They are gaps in that 
they are expensive and/or difficult to monitor and/or contain information that is not readily 
available at present. 


Other Human Usage Metrics 


A large number of variables and ecosystem descriptors are being ignored. In particular, these 
include human influences such as seismic activity, pollution, ballast water, etc. They may be 
addressed once the basic biological features have been fully addressed. 


METHODS 


All methods have been implemented in R, an open-sourced programming environment. The 
methods are shared via GitHub using the git revisioning system. These architectural choices 
were adopted to enhance the transparency and ease of sharing and collaborating with all 
interested parties. It is structured such that any additional data series can be easily added to the 
system to permit adaptive change. In this way, the approaches developed represent a true 
structural framework in which to further develop methods and approaches. 


The main methods used/developed in this report will be described in this section. 


BIODIVERSITY AND TAXONOMIC RICHNESS 


Biodiversity is seemingly a simple concept that is, in fact, fundamentally complex. This is 
because it ranges in focus from genetic and phenotypic variations within a local population, to 
breeding populations, to biome or even larger scaled genetic and phylogenetic and community 
variability, both in terms of their number and relative dominance. Any and all of these aspects of 
biodiversity can be estimated. However, it is the number of unique kinds of organisms found in a 
given location (commonly called taxonomic richness) that is most readily quantified and 
monitored. 


Taxonomic richness is known to increase asymptotically with sampling intensity. As such, a 
statistical correction (“rarefaction”) for spatial and temporal sampling intensity must be applied to 
be meaningfully comparable across locations and time. Specifically, a simple regression model 
was used to predict an expected richness R at a standard surface area and time depth: 


R ~ Lognormal((f + a,log(SA) + a;log(TS)), T) 


The other terms are f a constant, SA surface area (ranging from 1 to 50 km, radial length 
scale), TS, the number of years entering into the count (ranging from 0 to 5 years) and t is a 
lognormal error. The autocorrelation in SA and TS are ignored for the present but will eventually 
be modeled as well via a Poisson process. 
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The intent is to model the spatial/temporal patterns and then integrate them in a risk-based, 
probability model to permit formal statements of risk and probability of exceeding thresholds. 


PRODUCTIVITY 


Total system standing biomass is generally used as a proxy for productivity. They are not the 
same; however, they will be used to describe aggregate abundance of: various categories of 
organisms such as total bottom biomass, macroinvertebrates, zoopankton, phytoplanton, 
chlorophyll-a, etc. 


To estimate true production, a modeled approach is necessary. These indices can then be 
coupled with spatially explicit total landings to estimate the biomass and secondary production 
associated with the biota and fishery exploitation/footprint. 


The intent is, therefore, to model the spatial/temporal patterns and then integrate them in a risk- 
based approach to permit formal statements of risk and probability of exceeding thresholds and 
perhaps even estimate net production. 


HABITAT 


The basic Hutchinsonian notion of "ecological niche" is closely tied to our current understanding 
of “habitat”. It is a multi-dimensional concept in that it incorporates an undefined set of 
environmental variables and the associated biological constraints/specialisations/requirements 
(e.g., nutrients, thermal, oxygen, pH, etc.) pertinent to an organism of interest. 


Two notions of habitat can be discriminated, depending upon outlook: 


e Functional (H;) — make increasing more precise habitat definitions by adding more 
environmental and biological factors for increasingly more precise categories of organisms. 


e Integrative (H;) — the biota living in a given time and location represent a full integration of all 
relevant environmental and biological factors at their proper space time scales and so, in 
effect characterises the full system-level concept of "habitat space". 


Functional-habitat Modeling 


A utilitarian way of describing the Functional-habitat (H;) space of an organism is to examine 
the presence-absence or relative abundance of organisms as a function of environmental 
gradients and biological/life history constraints. From such information, the likelihood of a given 
location to be potential habitat for an organism of interest can be derived. The most robust 
method is to develop a probability model under the assumption that the presence or absence of 
an organism is a Bernoulli process. This can be readily parameterised using standard 
Generalized Linear Models. However, as environmental constraints are almost always modal in 
influence given a wide enough environmental gradient, a nonlinear model is more useful. 
Generalized Additive Models and Stochastic Partial Differential Equation Models are two 
methods that can deal with these environmental constraints in a simple and efficient manner 
(Choi 2010; Appendix 3). 


The utility of such an approach is most relevant for organisms with highly specific habitat 
requirements. They are used in this framework. The intent is to develop such Functional-habitat 
models for key species of interest: wolffish, Cod, etc., to assess changes in their available 
“habitat”. Examples are provided from those derived in the Snow Crab assessment (see 
Results). Though this Functional-habitat concept itself does not exclude species interactions, in 
actual practice, they are generally ignored as they make the statistical estimation impossibly 
over-parameterized (but see Choi et al. 2012, for one possible solution). Further, as there will 
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always be factors that are either poorly known, poorly sampled, or poorly parameterized, 
(e.g., dissolved oxygen, pH, redox, bacteria, jellyfish, squid, pollution, substrate type, etc.), 
these models will always necessarily be incomplete. 


It is relatively straightforward to model the spatial/temporal patterns and then "integrate" them in 
a risk-based approach to permit more formal statements of risk and probability of exceeding 
thresholds. 


Integral Habitat - Whole System Level 


While the Functional-habitat concept is interesting and pragmatic, it is decidedly a reductionistic 
perspective. The monitoring and assessment requirements for an MPA also demands a whole 
system (phenomenological) perspective. The assumption we make is that the relative 
abundance of organisms found in a given location and time defines and mirrors the kind 
of habitat in which they live. Sessile organisms that require high flow environments and 
associated biota tend to exist and flourish with a given group of other organisms similarly 
adapted, and they are different from those that require cold waters and minimal water flow, etc. 


Thus, if we can quantify the observed species composition in a given location and time, we 
would, in effect, be describing the habitat. We will define this as Integral-habitat (H;): species 
assemblage information that directly reflects all biological and environmental interactions 
simultaneously, both measured and unmeasured, and some too fleeting to be measurable. 


Fortunately, these associations are readily quantified using multivariate methods of data 
ordination. Here we focus upon Principal Components Analysis, which focuses upon an 
eigenanalysis of correlational structure of the species assemblages. An example derived from 
the Snow Crab assessment is provided in the Results. 


The intent is to model the spatial/temporal patterns and then integrate them in a risk-based 
approach to permit formal statements of risk and probability of exceeding thresholds. 


CONNECTIVITY: SPACE AND TIME SCALES 


Spatial Scale 


Marine Protected Areas exist in a spatial context. The characteristic spatial scale of productivity, 
diversity and habitat found in the MPA will determine which processes will be relevant to these 
aspects of an MPA. If the spatial variations in the productivity of a species of interest is small 
relative to the size of an MPA, the chances of the MPA having an influence upon the species is 
enhanced. This is usually the case when short-range processes dominate (e.g., less mobile 
species, weakly dispersing, low currents, habitat heterogeneity at small scales). If, however, the 
spatial scale is larger than the MPA, then it would mean that broader/larger processes were 
influencing the productivity of the species (e.g., higher mobility or dispersal processes/current, 
and stronger spatial connectivity, habitat heterogeneity at larger scales) — resulting in a lower 
likelihood of the MPA having an influence upon the species or components of interest. 


A second important factor is the relationship of the characteristic spatial scale to monitoring. As 
organisms exist at a given spatial scale in a given area, a sampling/monitoring protocol must 
reference/address these spatial scales. For example, when a spatial feature (e.g., biodiversity) 
demonstrates short characteristic spatial scales (i.e., a lot of spatial structure at smaller scales), 
any sampling approach must respect this and similarly operate at such shorter scales or even 
smaller, if one is to be able to resolve the patterns and describe properly the subject of interest. 
oimilarly, if a feature (e.g., biodiversity) is long-ranged and one wishes to resolve the patterns 
properly, then a sampling protocol must be similarly long-ranged to resolve the pattern. A 
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sampling program much smaller than the characteristic spatial scale would be beneficial, but the 
accrued benefits relative to cost of sampling would be rapidly diminishing. In that time, effort and 
resources requirements generally increase more rapidly than any benefit, (e.g., in the simplest 
case, if one is looking only naively at standard error as a measure of benefit, then it would 
increase asymptotically with increased effort with a power of —1/2. 


For these fundamental reasons, defining the spatial scale of a given observation or process is 
imperative for the development of any assessment or monitoring of MPAs. To this end, we 
represent any spatially explicit observation as Y (s), which are measured in a coordinate space 
(s € $ € Rİ} and domain S of dimensionality d. In this framework, we will mainly focus upon the 
case of d — 2 spatial dimensions (e.g., longitude and latitude or northing and easting). The 
observations Y (s) are assumed to be realizations of a spatial stochastic process, y, that is, 
some latent, unobservable but real stochastic generative function. These are also known as 
(spatial) "Random Fields” in the literature. The manner in which the variability of y changes as a 
function of distance, h =|| s — s' ||, is known as the spatial autocorrelation function. This spatial 
dependence is highly informative in that it defines how the similarity of observations changes 
with distance and so ultimately defines spatial scale. 


The spatial model is succinctly expressed as a regression model of a stochastic process 
(Banerjee et al. 2014): 


Y(s) = a(s) + els) 


where, the observations Y (s) are a function of some mean process u(s), and a residual error 
process e(s). The latter are further defined as: 


u(s = x(s)'B 
e(s) = w(s)+e(s) 


where x(s) are spatially referenced predictors with associated parameters 5; and the residual 
error process is decomposed into spatial c (s) and nonspatial e(s) components. The latter is 
also known as “nugget” error in geostatistics and represents error associated with measurement 
and/or microscale variability/processes. 


The error structures are usually assumed to be the following: 


e(s) ~ N(0,t7) 
w(s) ~ GP(0,C(s,s';0)) 
Y ~ MVN(H Y). 


The nonspatial error is assumed normal with mean O and standard deviation t. The spatial error 
is assumed to follow a Gaussian Process with mean 0 and covariance C(s,s’; 0), that is, a 
spatial covariance function with parameters 0. A multivariate normal likelihood is usually 
assumed for the observations Y = (Y (s1), ..., Y (s,))!, with mean u = XB and covariance X = 


[C(S; Sy; 0). + 1?1,. The matrix of regressors is X = [x(s;)1]|#_, and I, is an identity matrix of 
size n. 

The spatial covariance function C(h) = C(s, s^) expresses the tendency of observations closer 
together to be more similar to each other than those further away; h =|| s — s' || is the distance 


separating observations. Historically, a number of different forms have been used. The most 
frequently used forms include (for h > 0): 
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3 1 
o (1—-75h/8 +> (/9)); 0<h<=¢ 


C(A) spherical 
0; isg 
CCh)exponential = o? (exp(-h/d)) 
CHGausstan = G^ (expC- (h/0)^) 
C(h)powered exponential = o° (exp(—|h/|”)) 
Cuan = 0 VEO KVEG): 


At zero distance, C(0) = Cov(Y (s), Y (s)) = Var(Y(s)) = t? + o* (i.e., global variance), where + 
is the nonspatial error, o is the spatial error, and 0 = ($,v, p, ... } are function-specific parameters 
including @ the range parameter. I'(-) is the Gamma function and K,,(-) is the Bessel function of 
the second kind with smoothness v. The latter, Matérn covariance function is frequently used as 
the shape of this function is more flexible (Figure 33), albeit at the cost of an additional 
parameter (Appendix 2). 


The semivariance, more commonly used in geostatistics, is the covariance reflected on the 
horizontal axis of the global variance: y(h) = C(0) — C(h) = : Var[Y (s) — Y (s^)]. 


The spatial autocorrelation function is defined as the covariance function scaled by the global 
variance: p(h) = C(h)/C(0). For the purposes of this framework, we will define the spatial 
scale to be the distance at which the spatial autocorrelation decreases asymptotically to p(x) ^ 
0.05 ("practical range). This will be determined in modeled form where possible and, when it is 
not possible, an empirical estimate will be made derived from the empirical semivariogram. We 
focus upon spatial patterns larger than 1 km in resolution and smaller than the size of study 
domain, the "Scotian Shelf Ecosystem" (SSE). 


Temporal Scale 


Marine Protected Areas also exist in a dynamic frame. As such, similar to the above spatial 
considerations, there also exists some characteristic temporal scale upon which an MPA and its 
subcomponents operate. If the overall temporal variations of the biota and environment is small 
relative to the overall life of an MPA, the chances of the MPA having an influence upon the 
species is enhanced. Presumably, the longevity of MPAs will be long lasting and so will 
guarantee some influence, however small that may be; an effect that would be enhanced if the 
subject is shorter-ranged in spatial scales. 


Again, similar to the spatial scale case, this also has a simple and obvious implication in terms 
of monitoring and assessment. Short-range variations require higher sampling effort to 
resolve/understand the issues and vice-versa. 


As the temporal scale is an informative metric for monitoring and assessment of an MPA, we 
must be precise in its definition. Similar to the spatial case, we focus upon how the correlation 
and variability of some quantity changes with greater difference in time. The analogue to the 
semivariogram in a time series context is known as a cumulative periodogram. 


A periodogram expresses the amount of variance found at different wavelengths (w). It is a 
discrete sample estimate of the continuous concept of spectral density, y(t): 


1/2 
y(t) = | ¿riot F (a0)doo. 
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It is easily obtained from a Fast Fourier Transform of any arbitrary time-series and so the 
cumulative distribution permits a rapid identification of the time scale at which correlation drops 
to some arbitrary level. To be approximately comparable to the spatial scale, we define the 
temporal scale as the time difference by which the Cumulative Power Spectral Density 
increases to 9596 of the total variance. 


If the goal is to resolve short-term processes, then sampling must also, of necessity, be more 
frequent. However, similar to spatial scale issues, there is a point where there will be 
diminishing returns for any increase in the resolution of a temporal signal. It is the intent of this 
framework to operate upon timescale of 1 year or greater. Sub-annual signals, where they are 
available, would be used to decompose the seasonal signals from the inter-annual signals to 
avoid bias due to discretization errors. 


Space-time Models 


In reality, spatial and temporal patterns coexist and evolve. They are correlated processes and, 
as such, a challenge to model properly ('nonseparable" — see below). This renders the 
independent treatment and estimation of autocorrelation in time and space problematic. 
Nonetheless, new developments in computational methods are bringing such models within 
range of use for a framework such as this. This is primarily due to efficient methods associated 
with numerical modeling of Stochastic Partial Differential Equations (SPDEs), especially in 
spectral (Fourier) space. 


Again, we follow Banerjee et al.'s (2014) development of spatio-temporal models as a simple 
extension of the spatial regression model. The observations, Y (s, t) are measured in a 
coordinate space (s,t) € $ € RI x R in the domain S of dimensionality d + 1. The space-time 
regression model can then be specified as: 


Y(s,t) — u(s,t) + e(s,t), 


where, p(s, t) = x(s,t)! B(s,t) is the mean process and e(s, t) the residual error. The 
parameters $ of the spatially referenced predictors x can have variable forms: 


e p(s,t) — varying in both time and space. 

e f(s) — spatially varying (fixed in time). 

e f(t) — temporally varying (and fixed in space). 

e Pf — completely fixed (no variation in time and space). 
e f.(s) + B,(t) — complex, hierarchical (mixed). 


The error process can be separated into a spatiotemporally structured component w and an 
unstructured component es: 


e(s,t) = w(s,t) + e(s, t). 


The unstructured error is usually assumed to be a white error process: e(s,t) ~ N(0, 02). 
However, the manner in which the spatiotemporally structured error is handled is not straight- 
forward: 


e w,(t) — temporal effects nested in sites (temporal evolution at each site, space is not 
modelled). 


e  (w¿(s)- spatial effects nested in time (spatial evolution at each time, time is not modelled). 
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e w(s)u(t) or w(s) + v(t) — separable (structure in space and structure in time are 
independent). 


e w(s,t) — non-separable (both time and space structure evolve in a nonsimple manner). 


If spatial and temporal autocorrelation act independently (“separably”) then: 
v(t) ~ GP(0,C(t,t';0,)) 
w(s) ~ GP(0,C(s,s';0,)). 


The spatial and temporal errors are usually assumed to follow Gaussian Processes with mean 0 
and covariance C(:,:: 8). The spatial covariance can be modelled with any spatial form such as: 
C(s, S5 0,) = C(Mmatérn- Similarly, the empora covariance can be formulated as any 
autocorrelation model such as: C(t, t’; m = ofexp(—¢,|t — t'|), or if discrete in time: 


AR(1): v(t + 1) = p,v(t) + n(t) with nt) — ? N(0, of), etc. 


While computationally appealing, even in this simple case of “separable” models, an evaluation 
of the likelihood requires the inverse X1, , which happens to scale with O(n") operations. This 
has been a strong bottleneck to further development of these covariance-based methods in 
large-scaled problems of space and space-time. Approximations have been suggested to 
overcome this crippling numerical problem: modeling the spatial process w(s) with a lower 
dimensional process via kernel convolutions, moving averages, low rank splines/basis functions 
and predictive processes (projection of spatial process onto a smaller subset; Selna and Switzer 
1996, Wikle and Cressie 1999, Hung et al. 2004, Xu et al. 2005, Banerjee et al. 2014); 
approximating the spatial process as a Markov random field with Laplace and SPDE 
Approximations (Lindgren and Rue 2015); and approximating the likelihood of the spatial- 
temporal SPDE process with a spectral domain process (Sigrist et al. 2015). 


In this framework, we focus upon Sigrist et al.'s (2015) approach due to the more realistic 
assumption of non-separability of space and time processes and the high computational 
performance of operating in spectral space. Specifically, the space-time error w(s, t) is 
formulated as an advection-diffusion process (SPDE): 


us t) = —u’Va(s,t) + V-XVo(s,t) — Co(s,t) + e(s,t). 


Here, s = (x, y)! € &^: u = (uy, uy)! parameterizes the drift velocity (advection); V = (E 
OF: y 


T | 
ox? 57 is 


the gradient sans V-F= S is the divergence operator for a given vector field F — 


pum COSQ Sind T COSQ sina 
Ge E) e EP TNR sina y: Cosa A sina y-:cosc 
diffusion via (y > 0, a € [0,7/2]) with D, > 0 parameterizing the diffusion range; € > 0 
parameterizing local damping; and e(s, t) parameterizing a Gaussian random field that accounts 
for source-sink processes with white noise in time and Matérn spatial covariance (aka, 
"innovation"). 


) parameterizes the anisotropy in 


The full advection-diffusion model specification in state space form is (Sigrist et al. 2015): 
Y(st) = x(st)g-o(s.t))-e(st) 


w(s,t) = @a(s,t) {discretized advection — diffusion process} 
a(s,t) = Ga(s,t—1)+é€(s,t) {transition model} 

e(s,t) ~ N(0,7%1) {unstructured error} 

€(st) ~N (0, Q) {innovation} 
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where, a(s, t) are Fourier coefficients; b = [@(s,), ..., o (sy)]/ is a matrix of spatial basis 


~ 


functions; G is the transition matrix; and Q is the innovation covariance matrix (residual errors). 
See Appendix 3 for more details. 


Once the form of the space-time model is formulated, it can be used as a predictive/interpolating 
method. Such interpolation is required to describe and understand the linkages between the key 
ecosystem attributes of productivity, biodiversity, habitat, and species of interest, as well as, to 
estimate the spatial and temporal scales. 


We use an R-package implementation of the above space and space-time (etmy, 
nttps //github.com/JaeU/stmv). The stmv library additionally permits localised, hierarchical, 
model-based interpolations, for a given (local) subdomain, indexed by 2. Note that y indicates a 
global regressor in a linear, nonlinear or binomial functional form f, and u, represents a local 
regressor specific to the subdomain 1. The model specifications for each data series are as 


follows: 


e Bathymetry (pure space model): 


(s) = (s)- €(S) 
wi(s) ~ GP(0,C(s,s';B,)) 
&(s) ~ N(0,t,(s)*) 
e Substrate grainsize (space with covariates): 
(Ss) = u(s) + @,(s) + £(s) 
u(s) = f( depth, slope, curvature ) 
wi(s) ~ GP(0,C(s,s’;8,)) 
&(s) ~ N(0,t(s)*) 
e Bottom temperature (temporal effects nested in sites): 
(St) = u(s) + O@A(S) + &(S) 
a(s) = f (depth ) 
(st) = a(s, t) 
a. (s, t) = G10 (S, t — 1) + €, (s, t) 
&(st) ~ N(0,:21) 
(st) ~ N(0,Q,) 
e Ecosystem indicators (including integrative habitat H;; temporal effects nested in sites): 
(st) = (s) + u1(s,t) + ox(s) + &(s) 
u(s) = f( depth, slope, curvature, substrate grainsize ) 
pa(s,t) = f,( temperature ) 
w,(s,t) = ),a,(s,¢) 
a. (S, t) = G,a,(S, t — 1) + €, (s, t) 
&(st) ~ N(0,171) 
(st) ~ N(0,Q,) 


e Abundance and functional habitat (Hs; temporal effects nested in sites): 
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(st) = p(s)+ 15, t) + 0,(s) + &(s) 


u(s) = f(depth,slope, curvature, substrate grainsize ) 
ua(s,t) = f,( temperature, ecosystem indicators ) 
w,(s, t) = a; (Ss, t) 
a, (S, t) = G,a,(s, C= 1) F €, (s, t) 
&(st) ~ N(0,171) 
(st) ~ N(0,Q,) 


Tagging and Mark-recapture 


A tangible way of quantifying time and space scales (connectivity) is to demonstrate movement 
and genetic similarity. In the latter, no effort has been made. In the former, due to synergies with 
the fishing industry, increased tagging efforts have been made in the vicinity of SAB. Most of the 
effort has been driven by industry, Ocean Tracking Network (OTN), and Emera interest in Snow 
Crab movement near SAB. However, acoustic tagging of other species of interest have been 
completed as well. In total, 80 tags were deployed: 58 Cod, 14 Atlantic Striped Wolffish, 

1 Shorthorn Sculpin and 7 Snow Crab. These tags will allow us to track these species over the 
coming years within the MPA (through the two existing receiver lines), as well as outside the 
MPA with other OTN receiver lines and potentially with an OTN glider. They will help define the 
spatial connectivity and range of the species of interest. The intent is to develop movement 
models where possible and estimate spatial range directly. 


Mark-recapture information for sea turtles, seals, sharks and various other species exist in the 
area. These data have not been examined nor are they always available and represent a data 
gap at present. 


RISK MODELING 


Risk means many things to many people. We use the term specifically in the sense of having a 
believable error distribution for some quantity of interest such that probabilistic inferences can 
be made. Once the error distributions are determined, it is simple to make probability 
statements related to how likely the current state is different from some previous state or 
arbitrary threshold. 


The estimation of such error distributions requires a reliable method to propagate the errors 
from observations to predictions. One venerable (deterministic) method is to build a 
mechanistic model and then using approximations or simulations to determine the error 
distributions of interest. A second (phenomenological) method empirically quantifies the errors 
and propagates them via statistical/correlational methods. The latter is chosen in this framework 
as it is very general and simple to implement. The former is not chosen as no operating model 
of an ecosystem nor any subcomponent is known to the authors that can be said to be able to 
perform with sufficient skill to be able to propagate errors, let alone, magnitudes. 


The proposed approach is to use the discrete form of the logistic equation for this purpose. 
Verhulst, Pearl and Lotka in the late 1800s and early 1900s popularized this equation, using it to 
describe patterns of asymptotic increase (population growth, economic growth, etc.). The model 
is sufficiently general that it can be readily applied to most quantities of interest that show some 
dynamic behaviour, such as, for example: aggregate estimates of biomass (productivity), 
biodiversity (biodiversity change), and habitat (habitat change). Similar techniques can be 
applied to the other indicators of interest and may be applied depending upon availability of 
time. The justification of why it can be so generally applied is developed in Appendix 4. 
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The discrete form the basic logistic equation is, after normalization to K, the "carrying capacity" 
(see Appendix 4): 


Ve = rye (1 — yc). 


Perhaps the most powerful and flexible method of estimation of the parameters, 0 = fr, K}, is a 
state space representation where an additional observation model is added to connect observed 
indices O to the unobserved and (usually unobservable) real system state y. The simplest such 
model is a linear scaling factor, q, though again others are possible, at the cost of more 
complexity: 


Ot = qyt 
where, 0 = ír, K,q]. 


The use of a Bayesian approach to solve the above nonlinear state space problem is frequently 
used as it provides an opportunity to: have greater numerical stability; incorporate prior scientific 
knowledge in a formal manner; realistically propagate credible errors; estimate unobserved 
(“true”) states; and simultaneously estimate model “process” errors and data "observation" 
errors. [Note that process errors (o) are the uncertainties that feeds back into future states via 


error propagation — e.g., via the recursive form of the logistic equation (i.e., errors in y,,4 in the 
state space of y, vs y,,4). They are important if predictive risk is being assessed. Observation 
errors (o2) refer to the uncertainties associated with measurement and observation 

(i.e., measurement/data-related errors of both variables in the state space of y, vs y,..).] 


This latter ability is particularly important as parameter estimates and forecasts based on 
observation-only errors provide unrealistically optimistic (small and constant) error bounds. 
Parameter estimates and forecasts based on process-only errors expand rapidly into the future, 
resulting in potentially unrealistically pessimistic (large and usually growing) error bounds. 


The posterior distribution of the parameters of interest, 9, conditional upon the data are 
estimated via MCMC (Gibbs) sampling using the JAGS platform (Plummer 2003). The JAGS 


ANTHROPOGENIC THREATS 


Productivity, habitat and biodiversity monitoring is essential if one is to assess if MPAs are 
achieving their primary objectives. It is also essential to examine the anthropogenic threats in 
the St. Anns Bank area and examine the cumulative impacts on productivity, habitat, 
biodiversity, and endangered or threatened species. From the data available, we can examine 
trawling and dredging disturbances, exploitation of marine resources by fisheries, fishing-gear 
entanglement threats to marine mammals and sea turtles, and vessel collision threats due to 
marine traffic, and vessel-noise disturbances. Each threat can be normalised on a zero-one 
scale to compare the intensity of threats across the region and then weighted and combined to 
examine cumulative anthropogenic threats as in Coll et al. (2012). 


RESULTS 


Most analyses have not been completed. The effort so far has been to assimilate and develop 
the scaffolding to support the analyses and future monitoring and assessment. These headings 
are here mostly as place holders. However, some preliminary results can be reported upon to 
show direction. We highlight a few of these results but emphasize their very preliminary nature. 
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BIODIVERSITY 


No analytical results are available for presentation at this time. 


PRODUCTIVITY 


At present, only the models for Snow Crab have been completed using this method (Figure 34; 
see Choi et al. 2012 for more details). 


HABITAT 


Functional Habitat 


At present, only the models for Snow Crab have been completed using this method (Figure 35 
and 36; see Choi et al. 2012 for more details). 


Integral Habitat 


An example of Integral habitat (H;) as expressed through an ordination of taxa found together 
in various bottom trawls (Figure 37). 


CONNECTIVITY 


Spatial Scale 


This is a first attempt at describing the spatial scale of SAB and outlying areas. Essentially, S. 
represents the distance one must walk before one loses memory of where one started 
(arbitrarily defined at the 596 level of autocorrelation). Using this level of similarity as the 
standard, depth variations are demonstrably patchy/clustered (Figure 38). In the SAB area, 
there is a mixture of large autocorrelation scales (exp(6)=400 km) in the north-east areas, while 
the south-west areas have spatial scales of approximately 20 km or less. 


Recall that if S. is small, then short-range processes dominate (e.g., less mobile species, weakly 
dispersing, low currents, habitat heterogeneity at small scales). If S. is large, then broader/larger 
processes were influencing the productivity of the species (e.g., higher mobility or dispersal 
processes/current, and stronger spatial connectivity, habitat heterogeneity at larger scales). 


As an another example, Snow Crab densities suggest different spatial autocorrelation patterns 
(Figure 39), with the majority of core fishing grounds being in the range of 60-90 km in length 
scale. In the vicinity of SAB, these spatial scales decline to about 10-60 km. 


Temporal Scale 


Similar to spatial scale, we have defined temporal scale (S,) as the time required to reach a 
state where autocorrelation drops to an insignificant level (5%). Recall that if S, is small, short- 
range variations dominate (higher sampling effort to resolve/describe). If S, is large, long-range 
variations dominate (lower sampling effort required to resolve/describe). No analytical results 
are available for presentation at this time. 


Tagging 
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DISCUSSION 


No discussions are forthcoming until more analysis can be completed. 


CONCLUSIONS AND RECOMMENDATIONS 


This report was an exercise in demonstrating what is available and possible. All methods 
discussed are viable and informative and so should serve as a strong basis for any future 
monitoring framework in SAB and MPAs in general. 
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Figure 1. Bathymetric chart of the St. Anns Bank (SAB) area with the proposed Marine Protected Area 
(MPA, thick line) and limited fishing zones (hatched areas). See Figure 2 for geographic location in a 
larger map. 
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Figure 2. Map of the data extraction area 37°N to 48°N and from 48°W to 71°W and the relative location 
of the St. Anns Bank MPA. 
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Discrete bottle data, 13855-2014 
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Figure 3. Monthly spatial distribution of discrete bottle data for the time period 1955-2014. 


Figure 4. Number of chlorophyll and nutrient profiles extracted from the BioChem database for each year 
since 1955. 
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Figure 5. Number of chlorophyll and nutrient profiles extracted from the BioChem database for the time 
period 1955-2014, grouped monthly. 
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Figure 6. Depth profiles of chlorophyll-a and nutrients; all data for the time period 1955-2014. 
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Total number of net deployments for each month during the time period 1999-2014. 
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Total number of net deployments for each month during the time period 1999-2014. 
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Figure 9. Monthly averages of all data from 1999 to 2014: total abundance (top left), total biomass 
computed from wet weight (top right), ratio of total biomass computed from wet weight to total abundance 
(bottom left) as a potential measure of the average weight of the individual organism, and abundance of 
Calanus finmarchicus, C. hyperboreus, and C. glacialis (botton right). 
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Figure 10. Chlorophyll-a concentration extracted from MODIS 8-day composite images for St. Anns Bank 
polygon for the time period 2002-201 5. 
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Figure 11. Average Chorophyll-a concentration computed from 8-day composite images for St. Anns 
Bank polygon for the time period 2002-2015. 
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Figure 12. MODIS semi-monthly chlorophyll-a concentration showing spring bloom progression in the 
Northwest Atlantic in 2012. Note the intense bloom at St. Anns Bank during the last two weeks in March. 
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Stánns MPA Polygon, Primary Production 2002-2014 
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Figure 13. Annual monthly Primary Production (PP) computed from PP composite images for St. Anns 
Bank polygon for the time period 2002-2014. 
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Figure 14. Average Primary Production (PP) computed from monthly composite images for St. Anns Bank 
polygon for the time period 2002-2014. 
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Figure 15. Sea Surface Temperature (SST) extracted from 8-day AVHRR composite images for St. Anns 
Bank polygon for the time period 1997-2015. 
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Figure 16. Average Sea Surface Temperature (SST) computed from 8-day AVHRR composite images for 
St. Anns Bank polygon for the time period 1997-2015. 
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Figure 17. Bi-weekly composites from AVHRR showing SST in the Northwest Atlantic in the spring of 
2012, corresponding to the intense spring bloom at St. Anns bank shown in Figure 2. 


38 
000227 


Agen 


EPN 


BHW GE GFW GEN GO SB^ WW SE 


Fractional year 


0,3 


06 


100 km 


Temperature 


1970 1920 1990 2000 2010 
Year 


Figure 19. Left: Survey locations in the Groundfish survey (orange) and Snow Crab survey (green). 
Right: Timing of surveys in the Groundfish survey (orange) and Snow Crab survey (green). 
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Figure 20. Number of sets in the groundfish surveys (blue) and the number of sets with usable net 
configuration data (orange). 
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Figure 21. Towed distance comparisons in the groundfish survey. 
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Figure 22. Net spread variations by year. Note in 2011, the doorspread sensors seem to have failed 
completely. Note also that wingspread has been significantly larger from 2013 to 2015. 
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Figure 23. Net spread variations: doorspread vs wingspread. Note also that wingspread has been 
significantly larger from 2013 to 2015 but not doorspread. The cause of this divergence is not known. 
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Figure 24. Left: Surface area estimates based on GSINF logged start-end positions vs computed surface 
area estimated from tow track and net configuration. Right: Surface area estimates based on GSINF 
logged start-end positions vs computed surface area estimated from tow track and net configuration, as 
well as, modeled solutions. 
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Figure 25. Commercial catch weights of Sea Scallops (Placopecten magellanicus) on Georges Bank, the 
Scotian Shelf, and in the Bay of Fundy. 
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Figure 26. Commercial catch weights of Atlantic Halibut (Hippoglossus hippoglossus) on Georges Bank, 
the Scotian Shelf, and in the Bay of Fundy. 
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Figure 27. Automatic Identification System (AIS) data collected from the Canadian Coast Guard terrestrial 
network of AIS receiving stations on December 8, 2015. A total of 127 vessels were detected in the area 
with each colour representing a unique vessel. 
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Figure 28. Bathymetic (100 m resolution) chart of the St. Anns Bank area with line of sight detection (red 
circles) for the terrestrial AIS receiving stations (red dots) around St. Anns Bank Area of Interest. 
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Figure 29. Bathymetic (100 m resolution) chart of the St. Anns Bank area with Simard et al. (2014) 
estimated vessel detection distances (blue circles) for the terrestrial AIS receiving stations (blue dots) 
around St. Anns Bank Area of Interest. 
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Figure 30. Hypothetical vessel positions (large filled circles) and interpolated vessel positions (lines) 
based on the A* algorithm for three vessels transiting through the St. Anns Bank area. Each colour 
represents a unique vessel. 
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Figure 31. Vessel density maps for the first quarter of a year based on satellite AIS data from 2013-2015 
for all vessels (left panel) and all vessels except of the Newfoundland ferries (right panel). 
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Figure 32. Cost maps developed for the A* function to interpolate undetected vessel positions as vessels 
transit in and out of the Gulf of St. Lawrence. 
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Figure 33. Matérn autocorrelation function, p(h) = C(h)/C(0), the covariance function C(h) scaled by the 
total variance C(0), for two values of v (dark lines). At nu = 100, it approaches the Gaussian curve (upper 
dark curve on the left side) while at v = 0.5 the curve is exponential (lower dark curve on the left side). 
The associated semi-variograms (scaled to unit variance) y (h) are shown in light stippled lines. Spatial 
scale is defined, in this framework, as the distance h at which the autocorrelation falls to 0.05% (dashed 
horizontal line) — in this example between 2.5 and 3 units, depending upon value of v. 
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RO.mass 2010 


Figure 34. Predicted biomass density of Snow Crab in Maritimes Region based upon a combination of a 
Functional-habitat method and simple spatial interpolation. 


Figure 35. Functional habitat, (Hs), the predicted probability of observing Snow Crab. 


47 
000236 
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Figure 36. Surface area of potential Functional habitat (H;) of Snow Crab in Maritimes Region. Note the 
large interannual variability and a decadal scale decline in the southern areas (lower two panels). 
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Figure 37. Integral habitat (H;) based upon species associations in Maritimes Region. Note the 
first axis is primarily a temperature gradient and the second associated with depth. 
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Figure 38. First estimate of log(spatial range; km) based upon depth variations. 


S ABUNDANCE RANGE 


Figure 39. A first estimate of log. (spatial range; km) of Snow Crab abundance variations. 
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APPENDICES 
APPENDIX 1. DATA QUALITY CONTROL OF AZMP DATA 


All is extraction, See Conte and processing methods are documented in R scripts found in 


DISCRETE BOTTLE DATA 


BioChem data are generally not subject to any quality control (QC). As such, substantial QC 
was required. The QC protocol used was based on procedures designed at DFO's Institut 
Maurice-Lamontagne (IML). These were in turn based on procedures developed by NOAA's 
National Oceanographic Data Center/World Ocean Database, as well as many of the tests 
proposed in the GTSPP Real-Time Quality Control Manual. The quality control procedure was 
as follows: 


Step 1: Impossible Dates 


Due to known issues with dates, database query for nutrients and chlorophyll were designed to 
extract records for which the sampling date is within start and end dates of the mission. Another 
check includes comparing HEADER START and EVENT START dates that should be the 
same. It is often found that month and day were reversed in EVENT START field. Those 
records were retained and HEADER. START was used as a more reliable date. 


Step 2: Quality Control Flags 


A small number of records in BioChem was subject to quality control and include flags for 
position (POSITION QC CODE) and data (DATA QC CODE). The meaning of the codes are 
as follows: 


e O- No quality control performed 

e 1 = Value appears correct 

e 2- Value appears inconsistent 

e 3- Value appears doubtful 

e 4- Value appears erroneous 

e 5- Value changed as result of quality control 


The QC flags were checked for parameter and inconsistent (2), doubtful (3) and erroneous (4) 
values were removed from the dataset. 


Step 3: Depth Check 


For bottle data, start and end depth at which the water samples were collected are verified to be 
the same. Records with different start and end depths were removed from the dataset. 


Step 4: Duplicated Records 


BioChem often contains duplicated records as the same data was sometimes loaded into a 
database twice and treated as different records. Duplicated records were removed, and the first 
record of each duplicate is kept in the dataset. 


Step 5: Suspect Missions 


Missions with suspect data were identified and removed. Those missions often show unusual 
data values (for example, integer numbers without decimal places with only some values out of 
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range), suggesting that the data for the whole mission might be compromised. The suspect 
mission for each parameter are following: 


e  Chlorophyll-a OC7908, 32G879008 
e Phosphate 18HU88026 
e Silicate 180167005, 311R26870 


Chlorophyll and all nutrients values were examined if they fell within expected limits for the 
Northwest Atlantic, which are adopted from IML quality control procedure (IML Test 2.1). The 
expected range of values are: 


e Chlorophyll-a: 0-50 mg/m? 
e Nitrate: 0-515 mmol/m? 
e Phosphate: 0-4.5 mmol/m? 
e Silicate: 0-250 mmol/m* 


Any values outside of expected range were removed for open ocean data only. Coastal data (up 
to 5 km from the coast) were not filtered using expected ranges as in coastal water chlorophyll-a 
and nutrients concentrations can be higher. 


Step 6: Profile Envelope 


Data for each parameter are checked if they fall within the expected limits by depth interval, as 
shown in Table A1.1 (IML Test 2.4). This test does not allow zero values for silicate and 
phosphate in the deep water. Again, only open ocean data were subject to the profile envelope 
test. 


Table A1.1. Expected ranges of parameters for the profile envelope test (IML Test 2.4) 


Parameter Depth Interval Expected Range 
Chlorophyll-a 0-1500 m 0-50 mg/m? 
Silicate 0-150 m 0-250 mmol/m? 
Silicate 150-900 m 0.01-250 mmol/m? 
Phosphate 0-500 m 0-4.5 mmol/m? 
Phosphate 150-1500 m 0.01-4.5 mmol/m? 
Nitrate 0-1500 m 0-515 mmol/m? 


Step 7: Impossible Profiles 


This check was not implemented in the code and impossible profiles were identified by 
investigating unusual outliers. 


Additional steps from IML QC procedure, such as checks for constant profile, excessive 
gradient and inversions, were not implemented in this quality control procedure. However, due 
to eutrophication from terrestrial sources, phosphate levels in coastal regions often exceeded 
the upper limits for globally and locally possible values, with the phosphate concentrations 
sometimes six times higher than the upper limits for the Northwest Atlantic in offshore waters. 
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Therefore, coastal and open ocean data were examined separately; non-coastal data were 
filtered using the expected limits for the Northwest Atlantic. Coastal data were defined as the 
ones collected less than 5 km away from the coast, where a 5 km limit was chosen as an 
optimal distance at which all coastal inlets are included. Buffer polygons along the coastline 
were created (Figure A1.1) and used for flagging the data as open ocean records (Flag 1) and 
coastal records (Flag 2). 


FEN 
= ios 
— 
- 
Polygons: 
# 1 Mainiand 
a 2: Anticosti Island 


3: Magdalen Island 
* 4 Sable Island 
* $ Newfoundiand 


Longitude 


Figure A1.1. Polygons used for separation of coastal and ocean data, arbitrarily assumed to be a 5 km 
buffer from the coastline. 


CHLOROPHYLL-A 


Chlorophyll-a data are derived from four methods. The methods are listed and described in 
Table A1.2 and the aggregate time series associated with each method is shown in Figure A1.2. 
For most of the chlorophyll data, the method is not specified (unknown); Holm-Hansen 
fluorometric method is the standard AZMP method and is the second most frequent; 
Welschmeyer fluorometric method is used least frequently, often by the Quebec and 
Newfoundland regions. 


In a number of cases, the same water sample was processed using two different methods, 
resulting in two sets of chlorophyll estimates for the same samples. Comparisons between 
these two sets of values are shown in Figure A1.3. In both cases, chlorophyll-a estimated by the 
Welschmeyer method are lower than the ones using the Holm-Hansen method or “unknown” 
method. Since there is more data mapped to the Holm-Hansen method than to the 
Welschmeyer method, only data derived from Chl a and Chl-a Holm Hansen. sF methods 
were retained. No corrections were applied to correct for differences in methodology. 
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Table A1.2. Methods associated with chlorophyll-a records in BioChem. 


Method Description 

Chl a Unknown method 

Chl a Holm-Hansen F Holm-Hansen method; Prefiltered; Frozen before analysis (-20 ) 
Chl a Holm-Hansen sF Holm-Hansen method; Super Frozen before analysis (-196 ) 
Chl a Welschmeyer sF Welschmeyer method; Super Frozen before analysis (-196 ) 
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Figure A1.2. Time series of chlorophyll-a data from BioChem grouped by methods. 
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Number of overlapping samples: 1115 Number of overlapping samples: 1204 
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Figure A1.3: Comparison of the chlorophyll-a values collected using different methods is shown on the 
scatter plots on the top panel and the geographical locations of those samples is shown on the maps in 
the bottom panel. 
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NITRATE 


Nitrate estimates are derived from 10 methods. The methods are listed and described in 

Table A1.3 and the time series of data associated with each method is shown in Figure A1.4. 
Most of the methods measure nitrate and nitrite together. We also included data for nitrate only 
since in most seawater the concentration of nitrite is small compared to that of nitrate. 


Table A1.3. Methods associated with nitrate records in BioChem. 


Method Description 
NO2NO3 0 Nitrate+Nitrite / Unknown method 
NO2NO3_Alp_F Nitrate+Nitrite / Alochem / Frozen 


NO2NO3_Alp_SF 
NO2NO3_S8P1968 
NO2NO3_Tech_F 
NO2NO3_Tech_Fsh 
NO2NO3 Tech. SF 
NO2NO3 Tech2 F 
NO3 Tech F 


NO3 Tech SF 


Nitrate+Nitrite / Alochem / SuperFrozen 
Nitrate+Nitrite/ S&P(1968) / filtered and frozen 
Nitrate+Nitrite / Technicon / Frozen 

Nitrate + Nitrite / Technicon / Fresh / Strain / Unfiltered 
Nitrate+Nitrite / Technicon / SuperFrozen 
Nitrate+Nitrite / Technicon2 / Frozen 

Nitrate / Technicon / Frozen, corrected for NO2 


Nitrate / Technicon / SuperFrozen 
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Figure A1.4: Time series of nitrate data from BioChem grouped by methods. 


PHOSPHATE 


Phosphate data available in BioChem are mapped to 7 methods. The methods are listed and 
described in Table A1.4 and the time series of data associated with each method is shown in 
Figure A1.5. 


56 


000245 


Table A1.4. Methods associated with phosphate records in BioChem. 


Method Description 

PO4 0 Phosphate / Unknown method 

PO4 Alp SF Phosphate / Alpchem / SuperFrozen / Filtered 
PO4 Tech 2 Phosphate / Murphy and Riley / filtered and frozen 
POA Tech F Phosphate / Technicon / Frozen / Unfiltered 

PO4 Tech Fsh Phosphate / Technicon / Fresh / Strain / Unfiltered 
PO4 Tech SF Phosphate / Technicon / SuperFrozen / Filtered 
PO4 Tech2 F Phosphate / Technicon2 / Frozen / Unfiltered 


Phosphate [rmmalim3] 


Time 


Figure A1.5: Time series of phosphate data from BioChem grouped by methods. 
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SILICATE 


Silicate data available in BioChem are mapped to 8 methods (Table A1.5). Their time series of 


are shown on Figure A1.6. 


Table A1.5. Methods associated with silicate records in BioChem. 


Method Description 

SiO4_0 Silicate, Unknown methods and handling 
SiO4_1 Silicate / Mullin and Riley / filtered and frozen 
oiO4 Alp F oilicate / Alpochem / Frozen / Unfiltered 

SiO4 Alp SF Silicate / Alpchem / SuperFrozen / Filtered 
SiO4 Tech F oilicate / Technicon / Frozen / Strain / Unfiltered 


SIO4 Tech Fsh 
SiO4 Tech SF 


SiO4 Tech2 F 


Silicate / Technicon /Fresh /Strain / Unfiltered 
oilicate / Technicon / SuperFrozen / Filtered 


Silicate / Technicon2 / Frozen 
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Silicate [mrnolférmn] 


Time 


Figure A1.6. Time series of silicate data from BioChem grouped by methods. 


ZOOPLANKTON 


Zooplankton data was extracted from DFO's BioChem database (Devine et al. 2014) from 1914 
to 2014. They were comprised of 687 missions and 53,787 samples, using 13 different kinds of 
nets, 35 different mesh sizes ranging from 20 microns to 4.23 mm, and various net deployment 
and sample processing protocols. To ensure data consistency and comparability, only samples 
collected and analyzed using Atlantic Zone Monitoring Program (AZMP) protocol (Mitchell et al. 
2002) were retained for study, reducing the time scope to 1999-2014. 


As AZMP samples are not always properly flagged in the BioChem database, a list of missions 
that followed the AZMP protocol were provided by Ocean and Ecosystem Science Division 
(OESD) and Ocean Data and Information Services (ODIS). The relevant missions include 
AZMP spring and fall cruises, summer and winter groundfish survey missions, bi-weekly 
sampling at fixed stations (Halifax Station 2 and Prince 5) and samples collected on the Scotian 
Shelf during Labrador Sea missions. 


AZMP protocol samples zooplankton with Ring nets (0.75 m diameter, mesh size of 202 um) 
deployed as vertical tows from either near bottom or 1000 m (whichever is shallower) to the 
surface. The sample analysis includes estimation of abundance, species composition and 
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biomass in terms of wet and dry weight in two size fractions; one for organisms ranging from 
0.2 mm to 10 mm in size and the other for all organisms larger than 10 mm. The protocol is as 
follows: 


Organisms larger than 10 mm are manually separated from the sample, identified and 
counted. Wet weight is determined and reported for each individual species. In addition total 
wet weight for all large organisms is reported as the sum of all individual wet weights. 


Captured organisms smaller than 10 mm (0.2 - 10 mm fraction) are identified and counted. 
Dry and wet weight is determined and reported for the whole sample containing all 
organisms in that size fraction. 


Total wet weight is reported for all captured organisms as the sum of wet weights of large 
and small organisms. 


Developmental stages are identified for Calanus finmarchicus, C. glacialis and 
C. hyperboreus. 


oince the data hosted in BioChem is generally not subject to quality control and can contain 
erroneous records, substantial quality control was conducted to ensure correct representation of 
actual measurements. The quality control included verification of the following fields: 


time stamps: comparing mission dates with header dates and event dates start and end 
depths of the nets that cannot be equal or close together; 


volume of the samples: all records with volumes O, or NA were removed; 
split fraction of the sample: all the records with split fraction NA, O or » 1 were removed; 
minimum and maximum sieve for dry weight records: records NA sieve were removed; and 


repeated records. 


Finally, the numerical and biomass density for each species per unit surface area of a tow was 
computed as follows: 


abs(depth — depth 
abundance = counts x abs(deptisrart — depthena) [individuals/m7] 
split fraction x volume 


abs(depthstart — deptheng) 
split fraction X volume 


[g/m], 


biomass = weight x 


where counts refer to number of organisms encountered in the sample, depth star: and depth, 
to the start and end depth of the net deployment, split fraction to the fraction of the sample 
analyzed and volume to the sample volume. The final filtered dataset includes 126 missions in 
the time period 1999 to 2014, with 2,367 net deployments and more than 400 taxonomic 
species. 
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APPENDIX 2. MATÉRN FUNCTION 


The Matérn correlation function (Figure 33) is parameterised in a number of ways, and 
nomenclature of variables have also been inconsistent and can potentially cause confusion. The 
geoR library (Diggle and Ribeiro 2007; Ribeiro and Diggle 2001) defines it as: 


= area e) «l(5) 
p(x) zd 2K-1]['(x) d K d J 
where, œ > 0 is the “range parameter”; k > 0 is the shape (smoothness) parameter; K,.(-) is the 


modified Bessel function of the third kind of order x; and T(>) is the Gamma function. It is also 
related to the fractal dimension of the surface complexity (Constantine and Hall 1994). 


spBayes's (Finley et al. 2007) parameterization just a little different as well as is the 
nomenclature: 


p(x) — (px)" Ky, (px). 


1 
2V71 l'(v) 
INLA's (Lindgren and Rue 2015) parameterization is the same as spBayes’, but nomenclature is 
slightly different: 

1 a 
p(x) = Dt Fay 0%) Kj (kx). 
Thus, the following identities are important to interpret the discussions and outputs from these 
authors: 


KgeoR VspBayes = AINLA 


1/Pgeor 
Additionally, INLA defines range, 39, = V82/x (see rationale below). 


DspBaves = KINLA: 


The Matern covariance function, y (x) is the correlation function p(x) scaled by the variance. An 
offset of 7^ the variability at local scales associated with sampling error smaller than the unit of 
measurement (sometimes called “nugget” variance) is also modeled: 


y(x) = 1^ + o*(1 — p(x)). 
When Vsppayes = 1/2, this become the exponential model. When v;,payes > ©, the model 
becomes the Gaussian model. The curves become only incrementally different once Vspgayes > 
2 and so Finley et al. (2007) suggests limiting the prior to the interval (0,2) as a pragmatic 
solution to speeding up MCMC convergence. We have used the interval Vspgayes = (0,5), Just to 
be sure. Further, the effective range of spatial dependence (distance at which the correlation 
drops to 0.05 is given by —In(0.05)/Psprayes (Finley, Banerjee, and Carlin 2007). As such, they 
recommend a uniform prior in the interval of (Pspgayes: —In(0.05)). 


The spatial scale we define as exactly this practical or effective spatial range: the distance at 
which spatial dependence drops asymptotically to p(x) — 0.05. 


One of the reasons for the popularity of the Matérn covariance is that it is a stationary solution to 
the SPDE: 

(x? — 4)*"? (s£(s)) = W(s), 
where W (s) is Gaussian spatial white noise and t controls the variance and x is the scale 


parameter. INLA defines the effective range at V8/PspBayes: the approximate distance where 
correlation falls to 0.13 (vs 0.05): 
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A =a-d/2 
(A) 


2 


ü rT'(a)(4m)84/2<k22+2 


which for d=2 dimensions: 


A =a-1 
I'(A) 


» l'(aà)(41)k?^?' 


and a = 2 by default and so A = 1 such that: 
range =v8/x 
1 


"E 


Ank?t1?' 


range = V8/k. 
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APPENDIX 3. ADVECTION-DIFFUSION STOCHASTIC PARTIAL DIFFERENTIAL 
EQUATION (SPDE) 


Following the development of Sigrist et al. (2015), we begin with the general regression model 
Y(s,t) = u(s,t) + els, t), 


where, p(s, t) = x(s,t)! f(s, t) is the mean process and e(s, t) the residual error. This error 
process can be separated into a spatiotemporally structured component w and an unstructured 
component e: e(s, t) = w(s,t) + e(s,t). The unstructured error is usually assumed to be a white 
error process: <(s,t) ~ N(0, o). 


The structured error w(s, t) can be defined in terms of the following advection-diffusion SPDE: 


Se (s.t) = -u Va(s t) + - 270(5,t) — (S t) + e(s,D, 


where, s = (x, y)! € «^: u = (uy, uy)” parameterizes the drift velocity (advection); V = (E 22 rok 
is the gradient — V.F- E, y is the divergence operator for a vector field F = 

ot Ls cosa sind `r, cosa sina 
(Bo Ey) ; Fr pz a sing y: cosa? a: sina y: cosg 
diffusion via (y > 0, æ € [0,7/2]) with D, > 0 parameterizing the diffusion range; ¢ > 0 
parameterizing local damping; and e(s, t) parameterizing a Gaussian random field that accounts 
for source-sink processes with white noise in time and Matérn spatial covariance (aka, 
"innovation"). 


) parameterizes the anisotropy in 


If e(s, t) is stationary and has the form of a Whittle spatial covariance function, then it has the 
following spectrum: 


CE ME 


where k are the spatial wave numbers. The spectrum of the w(s,t) process is then: 


f (k) = 


f(w,k) = f (K) 5 ((K7 Ek + 7)? + (w + u! k)2) 2, 


B m) 
where w is temporal frequency. The covariance function can be recovered as: 
C(s, t) $ f(w, k) exp(i - tw) exp(i - s'k) dkdw 

~  exp(-i-u'kt — (k' Xk + Q)|t|) 
POUR AR 


exp(i - s'k) dk. 


2(kTXk + O 
Upon discretization: 
o(st) = o(s)"a(t) 
K 
= > 4040 
j=1 
4 K/2+2 
= De (SDAP + DOGO PERO, 
j=1 J=5 
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where k = (kyj, ky. ¡) is the spatial wavenumber of the j component of K Fourier components 
and $;(s) = exp(t: kj s) = (cos(k; s), sin(k; s)} is the spatial function. And, for each time t and 
location s;, | = 1, ..., n? in the SPDE: 


WG; (s) = -uko (s) 
wd (s) = -ukh (sp, 
and, 
D-20p (5) = -kf Ekg; (s) 
V.XV$P (s) = -k7zkof™ (sp. 
The full non-separable, advection-diffusion model specification is thus: 
Y(st) = x(st)g-o(st)-s(st) 
w(s,t) = da(s t) {advection — diffusion model} 
a(s,t) = Ga(s,t—- 1)+€(s,t) {transition model} 
e(s,t) ~ N(0,7%1) {unstructured error} 
€(st) ~ N(0,Q) {innovation}, 


where, a(s, t) are Fourier coefficients; b = [@(s,), ..., o (sy)]/ is a matrix of spatial basis 


functions; G is the transition (propagator) matrix; and Q is the innovation covariance matrix 
(residual errors). The Fourier functions are: 


$(s;) = (costs). ..., COS(K4S)), cos(K¿s; ), sin(K£s;) ..., cos CN "sin CN | 


2 
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APPENDIX 4. LOGISTIC MODEL 


The derivation of the logistic model presented below is more phenomenological than 
mechanistic and attributed to Lotka (1925). It is our opinion that this phenomenological 
interpretation renders it useful as a general characterisation of system state. The argument is 
that the rate of change in time t of any state Y (in our context: abundance, species richness, 
“habitat”, size) can be expected to be in some way, a function of itself: 


dY /dt = g(Y). 


It is expected that when Y = 0, dY /dt will also be zero and so represents an algebraic root of g. 
A Taylor series expansion of g near this root Y — 0 gives: 


dY/dt =  g'(0)Y + g"(0)Y?/2 + higher order terms ...; 
Y[g (0) + g(0)Y/2]. 


With the identities g'(0) = r and g"(0) = —2r/K, the standard form of the logistic equation is 
obtained: 


à 


dY /dt = rY(1 — Y/K). 


The intrinsic rate of increase, r, is, therefore, some abstract and aggregate function that 
describes the net increase or decrease of the system state Y when Y is small. In biological 
systems, this means a maximum rate of growth, recruitment, mortality, movement, climatic 
change, extinction, speciation, etc. as the rate is expected to decline as it approaches some 
upper limit K of the magnitude of the system state Y. 


With normalization by K such that y = Y/K: 


dy/dt = ry(1— y). 


Many variations of this basic model are known, mostly different ways of adjusting the shape of 
the curve and/or the location of the inflection point (i.e., K/2). A flexible family is the gamma- 
logistic: 


dy/dt = ry*(4-— y)P + f(...). 


The additional parameters exponentiate different components which ultimately amounts to 
adding higher order and even fractional order terms in the Taylor series) and also adding 
additional terms f(... ) that are external to the dynamic that r and K govern, such as fishing, 
advection, diffusion, noise. For the purposes of this discussion, we will use only the basic 
model and estimate the parameters of interest, 0 = fr, K}, but other parameterizations are 
possible, such as 0 = {r,K,a, p). The intent is to explore the utility of this more flexible 
formulation in conjunction with the basic model. 


In discrete form, where At — 1 year and after a Euler discretization, the basic model becomes: 


Yt = Tyra (1 — Yt-1). 
This, we call the "basic" form of the discrete logistic model. 
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